Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,171 T→A 15.7% intergenic (+43/‑187) tufB → / → secE translation elongation factor EF‑Tu 2/preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1710TA15.7% 266.9 / 39.5 108intergenic (+43/‑187)tufB/secEtranslation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (42/49);  new base A (9/8);  total (51/57)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 119 total aligned reads displayed.

TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATA  >  NC_000913/4177073‑4177269
                                                                                                  |                                                                                                  
tCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAtt                                                                                                   <  1:1615367/100‑1 (MQ=255)
 ccGTGAAGGCGGCCGTGCAGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATGGCCGATAACCTTTGACGCGATGCGGGCTAAAAGGGCATCAtt                                                                                                   >  1:1343333/1‑99 (MQ=25)
 ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAtt                                                                                                   >  2:881845/1‑99 (MQ=255)
  cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATCCGCACTAAAAGGGCCTCATTTg                                                                                                 >  1:2858219/1‑100 (MQ=255)
      aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGc                                                                                             <  2:585768/100‑1 (MQ=255)
        ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                           >  2:679476/1‑100 (MQ=255)
         gcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCt                                                                                          <  2:1644190/100‑1 (MQ=255)
          cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtt                                                                                         <  2:111061/100‑1 (MQ=255)
           ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                        <  1:1576828/100‑1 (MQ=255)
           ggcCGGACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                        <  1:1988795/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                       <  2:1023706/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                       <  2:2558436/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                       <  2:800973/100‑1 (MQ=255)
             ccGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                      >  2:2382859/1‑100 (MQ=255)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttag                                                                                     <  1:1390342/100‑3 (MQ=255)
               gTACCGTTGGGGCGGGCGTTGTAGGAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                       <  1:2722127/97‑1 (MQ=39)
               gTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACGTTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagt                                                                                    <  2:1919198/100‑4 (MQ=38)
               gTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGc                                                                                    <  2:1839663/100‑1 (MQ=255)
                tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACAGTTGACGCAAAGCGCACTAAAAGGGAATCATTTGATGCCCTTTTTGCa                                                                                   >  1:1961141/1‑100 (MQ=255)
                    gTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCt                                                                               <  1:2444542/100‑1 (MQ=255)
                    gTTGGCGCGGGCGTTGTAGCAAAAGTGCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGCGCATCATTTGATGCCCTTTTTGCACGCt                                                                               >  1:1006889/1‑100 (MQ=255)
                       ggCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                            >  2:392617/1‑100 (MQ=255)
                        gcgcGGGCGTTGTAGCAAAAGTTCTGCGCTAATTGCCGATAACATTTTACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                           <  1:1692353/100‑1 (MQ=255)
                        gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                           <  2:1701650/100‑1 (MQ=255)
                         cgcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGt                                                                          >  2:195704/1‑100 (MQ=255)
                          gggggCGTTGTAGCAAAAGTTGTGAGCTAATTGTCGTTACGATTTGACGCTATGCGCCCTAAAAGGGCATCATTTGATGTGCTTGTTGCACGCTTTCg                                                                           <  1:2426915/96‑1 (MQ=255)
                          gcgGGCGTTGTAGCAAAAGTTCTGTGCTAATTGCCGATAACATTTGACGCAATGCGGACTAAATGGGCATCATTTGATGCCCTTGTTGCACGCTTTCGTa                                                                         <  1:891196/100‑1 (MQ=255)
                           cgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAc                                                                        <  1:788782/100‑1 (MQ=255)
                           cgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAc                                                                        <  2:1642723/100‑1 (MQ=255)
                           cgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAc                                                                        >  2:621417/1‑100 (MQ=255)
                             ggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                      <  2:707118/100‑1 (MQ=255)
                             ggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                      >  1:1167601/1‑100 (MQ=255)
                               cGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGa                                                                    >  1:895120/1‑100 (MQ=255)
                                  tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                 <  2:2094336/100‑1 (MQ=255)
                                    tAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                               >  1:2346512/1‑100 (MQ=255)
                                     aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                              >  1:2251325/1‑100 (MQ=255)
                                     aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                              <  2:672532/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                             >  2:1496130/1‑100 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                             >  1:2772793/1‑100 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                             <  2:2833935/100‑1 (MQ=255)
                                       caacggctCTGAGCTAATTGCCGATAACATTTGAAGTAATGCGCACTAAAAGCGCATCATTTGATGCCCTTTTGGCACGCTTGCGTACCAGAACCt                                                                >  1:2881732/8‑96 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGAAATGCGCACTAAAATGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                            >  1:1365067/1‑100 (MQ=255)
                                          aaGTTCTGAGCTAATTGCCGATAACATTTGACACAATGCTCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcat                                                         <  1:2041210/100‑1 (MQ=255)
                                            gTTCTGAGCTAATTGCCGATAACATTTGACGAAATGCGCAATAAAAGAGCATCATTTAATGCCGTTTTTGAAAGCTTTCGTAACAGAACCTGGCtc                                                           >  1:711045/1‑96 (MQ=255)
                                               cTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTg                                                    >  2:2039327/1‑100 (MQ=255)
                                                tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGa                                                   >  1:2013950/1‑100 (MQ=255)
                                                 gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                  <  2:883850/100‑1 (MQ=255)
                                                 gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                  >  1:1247270/1‑100 (MQ=255)
                                                   gCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAttt                                                >  1:758145/1‑100 (MQ=255)
                                                    cTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAtttt                                               <  1:921980/100‑1 (MQ=255)
                                                     tAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTc                                              >  1:414934/1‑100 (MQ=255)
                                                      aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                             <  2:2346523/100‑1 (MQ=255)
                                                       aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                            <  2:2402258/100‑1 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                          <  1:1063646/100‑1 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                          <  2:187495/100‑1 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                          <  1:2176865/100‑1 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                          >  2:2560160/1‑100 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                          <  2:31040/100‑1 (MQ=255)
                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                         <  2:838120/100‑1 (MQ=255)
                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                         >  2:1982360/1‑100 (MQ=255)
                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGACGCCCTTTTTGCACGCCTCTGTACCCGAACCTGGCTCACCAGGCATTCTGTTTg                                          >  2:2637606/1‑99 (MQ=255)
                                                           ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTc                                        <  2:396659/100‑1 (MQ=255)
                                                            cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                       >  2:1076805/1‑100 (MQ=255)
                                                            cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                       >  1:1091789/1‑100 (MQ=255)
                                                               tAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATaa                                    >  1:1597009/1‑100 (MQ=255)
                                                                aaCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAAt                                   <  1:2526951/100‑1 (MQ=255)
                                                                 aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                  >  2:162307/1‑100 (MQ=255)
                                                                 aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                  >  2:2392741/1‑100 (MQ=255)
                                                                      tGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                             <  2:2110247/100‑1 (MQ=255)
                                                                      tGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                             >  2:1819459/1‑100 (MQ=255)
                                                                      tGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                             <  1:392615/100‑1 (MQ=255)
                                                                       gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCt                            <  1:724988/100‑1 (MQ=255)
                                                                         cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                          >  1:2671231/1‑100 (MQ=255)
                                                                         cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                          >  1:937320/1‑100 (MQ=255)
                                                                           cAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgaga                        >  1:711767/1‑100 (MQ=255)
                                                                            aaTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGAc                       <  1:2022158/100‑1 (MQ=255)
                                                                             aTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACa                      <  2:420249/100‑1 (MQ=255)
                                                                                cgtgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                   >  2:2572378/7‑100 (MQ=38)
                                                                                 gtgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                  >  1:2138202/6‑100 (MQ=38)
                                                                                  tgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                 >  1:645859/5‑100 (MQ=38)
                                                                                  cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                 <  2:818167/100‑1 (MQ=255)
                                                                                   gcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctct                >  2:1839865/4‑100 (MQ=39)
                                                                                   aCTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctct                >  1:1532514/1‑100 (MQ=255)
                                                                                    caaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg               >  2:2490692/3‑100 (MQ=255)
                                                                                    cTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctct                >  1:1090485/1‑99 (MQ=255)
                                                                                      aaaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACAAGAACCTGGATCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGATCTGtt             >  1:2458508/1‑100 (MQ=255)
                                                                                      aaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGtt             >  2:1656129/1‑100 (MQ=255)
                                                                                       aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg            <  1:383129/100‑1 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa           >  1:874006/1‑100 (MQ=255)
                                                                                         aGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg          <  1:2245797/100‑1 (MQ=255)
                                                                                         aGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg          <  2:293595/100‑1 (MQ=255)
                                                                                          gggCATCATTTGATGCCCTTTTTGCTCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg         <  1:51451/100‑1 (MQ=255)
                                                                                           ggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAggg        <  2:1021668/100‑1 (MQ=255)
                                                                                             catcatTTGATGCCCTTTTTTTTTTTTTTTGTACCAGAACCTGGCTCATCAGTGCTTTTCTTTGTCATAATCATTGCTTGGACAGTCTCTTTTGacctc       >  1:1119593/1‑95 (MQ=255)
                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg      <  1:162306/100‑1 (MQ=255)
                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg      <  2:1382831/100‑1 (MQ=255)
                                                                                               tcatTTGATGCCCTTTTTGCACGCTTTCTGACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGta    <  2:2426903/100‑1 (MQ=255)
                                                                                               tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGta    <  1:586261/100‑1 (MQ=255)
                                                                                               tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGta    >  1:2556293/1‑100 (MQ=255)
                                                                                                 atTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGtata  >  1:1311386/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATA  >  NC_000913/4177073‑4177269

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: