Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,221,886 A→T 6.4% Y154F (TAT→TTT)  yebT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,221,8860AT 6.4% 274.9 / 9.1 94Y154F (TAT→TTT) yebTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (50/38);  new base T (2/4);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCA  >  Exported/1221818‑1221956
                                                                    |                                                                      
tatCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAatat                                                                       >  1:2924238/1‑70 (MQ=255)
tatCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATc                                                                      <  1:2798219/71‑1 (MQ=255)
tatCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATc                                                                      >  1:1528017/1‑71 (MQ=255)
  tCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCgg                                                                    >  1:2140211/1‑71 (MQ=255)
     gCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                                 >  1:70409/1‑71 (MQ=255)
       atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                               <  1:250044/71‑2 (MQ=38)
       atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                               <  1:2652691/71‑2 (MQ=38)
       atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                               <  1:3955505/71‑2 (MQ=38)
       atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                               <  1:2802799/71‑2 (MQ=38)
        tgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                               >  1:2012182/1‑70 (MQ=255)
         gatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTgggta                                                             >  1:3587952/1‑68 (MQ=37)
         gatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTgggta                                                             >  1:906849/1‑68 (MQ=37)
          atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                            >  1:1029857/1‑71 (MQ=255)
          atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                            >  1:802471/1‑71 (MQ=255)
          atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                            >  1:2725512/1‑71 (MQ=255)
          atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                            >  1:3852980/1‑71 (MQ=255)
                 gtaaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                                     >  1:4177585/1‑71 (MQ=255)
                  taaaggtaaagAGCAGTATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                                       <  1:4320184/68‑1 (MQ=255)
                  taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                                       <  1:3412172/68‑1 (MQ=255)
                  taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                                       <  1:3559872/68‑1 (MQ=255)
                  taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATc                                                    >  1:2849333/1‑71 (MQ=255)
                    aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                                  >  1:3292777/1‑55 (MQ=255)
                    aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                                  >  1:3794254/1‑55 (MQ=255)
                    aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                                  >  1:4073706/1‑55 (MQ=255)
                    aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCt                                                   <  1:4350618/70‑1 (MQ=255)
                     aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                                 >  1:3089539/1‑71 (MQ=255)
                     aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                                 >  1:1707162/1‑71 (MQ=255)
                     aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                                 >  1:3518189/1‑71 (MQ=255)
                     aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                                 >  1:3657989/1‑71 (MQ=255)
                      ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                               >  1:1809656/1‑56 (MQ=255)
                      ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                               >  1:2567144/1‑56 (MQ=255)
                      ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                                <  1:3418502/71‑1 (MQ=255)
                       gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                                 <  1:859358/53‑1 (MQ=255)
                       gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                                 <  1:2656646/53‑1 (MQ=255)
                       gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                                 <  1:4046307/53‑1 (MQ=255)
                       gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                            >  1:2414956/1‑58 (MQ=255)
                       gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgatg                                               <  1:452570/71‑1 (MQ=255)
                       gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgatg                                               <  1:1226216/71‑1 (MQ=255)
                       gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                                >  1:4003162/1‑70 (MQ=255)
                        taaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGAGATCtgatg                                               >  1:2503251/1‑70 (MQ=255)
                         aaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                                     >  1:3519554/1‑63 (MQ=255)
                          aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                                     >  1:3959734/1‑62 (MQ=255)
                          aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                                     >  1:3108247/1‑62 (MQ=255)
                          aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATc                                            >  1:3429980/1‑71 (MQ=255)
                           agagCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCg                                                                     <  1:1423413/45‑1 (MQ=255)
                           agagCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCg                                                                     <  1:2750500/45‑1 (MQ=255)
                           agagCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCg                                                                     <  1:868985/45‑1 (MQ=255)
                              gCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCAc                                         <  1:2187188/70‑1 (MQ=255)
                               cAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCt                                       <  1:1546748/71‑1 (MQ=255)
                               cAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCt                                       <  1:1151494/71‑1 (MQ=255)
                                aGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTg                                      >  1:3115701/1‑71 (MQ=255)
                                 ggATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGc                                     <  1:3742953/71‑1 (MQ=255)
                                   aTCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGc                                     <  1:2464709/69‑1 (MQ=255)
                                        tttGTCGCACTCGATACCCAACCGAAATATCGGCTGGAc                                                              <  1:806404/39‑1 (MQ=255)
                                        tttGTCGCACTCGATACCCAACCGAAATATCGGCTGGAc                                                              <  1:1780079/39‑1 (MQ=255)
                                          tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                            >  1:1275829/1‑39 (MQ=255)
                                          tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                            >  1:3476280/1‑39 (MQ=255)
                                          tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                                     >  1:1969534/1‑46 (MQ=255)
                                          tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGccc                               <  1:1301739/68‑1 (MQ=255)
                                              gCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtct                        <  1:1955496/71‑1 (MQ=255)
                                               cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc                       >  1:3299012/1‑71 (MQ=255)
                                               cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc                       <  1:1100457/71‑1 (MQ=255)
                                               cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc                       >  1:191557/1‑71 (MQ=255)
                                                 cTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCgg                     >  1:2097516/1‑71 (MQ=255)
                                                 cTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCgg                     <  1:2811734/71‑1 (MQ=255)
                                                 cTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCgg                     >  1:77109/1‑71 (MQ=255)
                                                  tCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                                <  1:3488793/43‑1 (MQ=255)
                                                  tCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                                <  1:3575275/43‑1 (MQ=255)
                                                   cGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCAc                                         >  1:2631649/1‑49 (MQ=255)
                                                   cGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGtt                   >  1:2608283/1‑71 (MQ=255)
                                                    gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAg                                  >  1:2740496/1‑55 (MQ=255)
                                                    gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc                  >  1:3863571/1‑71 (MQ=255)
                                                    gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc                  >  1:4136526/1‑71 (MQ=255)
                                                    gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc                  >  1:1348553/1‑71 (MQ=255)
                                                    gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc                  >  1:1940407/1‑71 (MQ=255)
                                                    gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc                  <  1:2952791/71‑1 (MQ=255)
                                                      tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                             <  1:1265361/42‑1 (MQ=255)
                                                      tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                             <  1:3089195/42‑1 (MQ=255)
                                                      tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                             <  1:2289012/42‑1 (MQ=255)
                                                      tACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCa                                    >  1:741907/1‑51 (MQ=255)
                                                      tACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCa                                    >  1:1649786/1‑51 (MQ=255)
                                                        cccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCt               >  1:614977/1‑70 (MQ=255)
                                                        cccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTg              <  1:979830/71‑1 (MQ=255)
                                                        cccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTg              <  1:4404829/71‑1 (MQ=255)
                                                         ccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtct                        <  1:2667694/60‑1 (MQ=255)
                                                         ccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtct                        <  1:1911608/60‑1 (MQ=255)
                                                          cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa            >  1:3159878/1‑71 (MQ=255)
                                                          cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa            >  1:2990037/1‑71 (MQ=255)
                                                          cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa            >  1:2272879/1‑71 (MQ=255)
                                                          cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa            >  1:1774182/1‑71 (MQ=255)
                                                           aaCCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAAc           <  1:2911884/71‑1 (MQ=255)
                                                            aCCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa            >  1:3417214/1‑69 (MQ=255)
                                                            aCCGAAATATAGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa            >  1:1017712/1‑69 (MQ=255)
                                                             ccGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAg         <  1:814123/71‑1 (MQ=255)
                                                              cGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGc        <  1:495221/71‑1 (MQ=255)
                                                                   tatCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGTTTCGCTGAACAGCGGTTc   <  1:2176364/71‑1 (MQ=255)
                                                                   tatCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGtt    >  1:3391205/1‑70 (MQ=255)
                                                                    atCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCa  >  1:810067/1‑71 (MQ=255)
                                                                    |                                                                      
TATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCA  >  Exported/1221818‑1221956

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: