Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1710TA20.4% 202.7 / 52.8 93intergenic (+43/‑187)tufB/secEtranslation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (30/44);  new base A (9/10);  total (39/54)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.76e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 102 total aligned reads displayed.

CCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATA  >  NC_000913/4177074‑4177269
                                                                                                 |                                                                                                  
ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAttt                                                                                                  <  2:2280548/100‑1 (MQ=255)
   tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                               <  1:1925553/100‑1 (MQ=255)
    gAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                              <  1:1958568/100‑1 (MQ=255)
    gAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                              <  1:546955/100‑1 (MQ=255)
     aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGc                                                                                             <  1:1464509/100‑1 (MQ=255)
      aGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGcc                                                                                            >  1:2049535/1‑100 (MQ=255)
       ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                           >  1:1369870/1‑100 (MQ=255)
       ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                           >  2:21478/1‑100 (MQ=255)
       ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                           >  2:752346/1‑100 (MQ=255)
       ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                           >  1:1127068/1‑100 (MQ=255)
         cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtt                                                                                         <  1:1120074/100‑1 (MQ=255)
         cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtt                                                                                         <  2:1126586/100‑1 (MQ=255)
         cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtt                                                                                         <  1:1549590/100‑1 (MQ=255)
          ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttt                                                                                        <  1:319135/100‑1 (MQ=255)
          ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                        <  1:1137732/100‑1 (MQ=255)
           gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                       <  1:1267612/100‑1 (MQ=255)
           gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                       <  1:2373398/100‑1 (MQ=255)
            ccGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                      <  2:737147/100‑1 (MQ=255)
              gTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGc                                                                                    >  1:761088/1‑100 (MQ=255)
              gTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagt                                                                                    <  1:587729/100‑4 (MQ=38)
               tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTTTGAGCTAATTGGCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtg                                                                                   <  2:580295/100‑5 (MQ=38)
               tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                   <  1:2187949/100‑1 (MQ=255)
               tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtg                                                                                   <  1:1919961/100‑5 (MQ=38)
               tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtg                                                                                   <  1:1218648/100‑5 (MQ=38)
               tACCGGTGGCGCGGGAGTTGTGGCAAAAGTTCTGAGCTAAGTGCCGATAACATTTGGCGCAATGCGCACTAAAAGGGCATGAAATGATGCACttttaggg                                                                                   <  2:690861/100‑5 (MQ=255)
                   ggtGGCGCGGGCGTTGTAGCAAACGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCGCTTTTTGCACGCt                                                                               >  2:1866293/3‑100 (MQ=255)
                       gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                           <  2:1575791/100‑1 (MQ=255)
                       gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                           <  1:1252651/100‑1 (MQ=255)
                       gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                           <  2:1574707/100‑1 (MQ=255)
                       gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                           >  2:321921/1‑100 (MQ=255)
                       gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                           <  1:2386472/100‑1 (MQ=255)
                        cgcgGGCGTTGTAGCAAAAGTTCTGAGCTAATGGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGt                                                                          >  2:456563/1‑100 (MQ=255)
                        cgcgGGCGTTGTAGCAAAAGTTCTGAGATAATTGCCGATAACGTTTGACGCAATGCGCAGTAAAAGGGCTTCATTTGATGCCCTATTTGCACGCTTTCGt                                                                          <  1:2296552/100‑1 (MQ=255)
                         gcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTa                                                                         <  1:2458302/100‑1 (MQ=255)
                          cgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAc                                                                        <  2:595909/100‑1 (MQ=255)
                          cgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAc                                                                        >  2:268102/1‑100 (MQ=255)
                               gTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGaa                                                                   >  1:1582824/1‑100 (MQ=255)
                                 tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                 >  2:167050/1‑100 (MQ=255)
                                  gTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCt                                                                <  2:174491/100‑1 (MQ=255)
                                    aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                              >  2:1265275/1‑100 (MQ=255)
                                    aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                              >  1:579110/1‑100 (MQ=255)
                                    aGCAAAAGGTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCCTTCGTACCAGAACCttg                                                              >  2:1382393/1‑98 (MQ=255)
                                     gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                             <  1:45852/100‑1 (MQ=255)
                                      cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                            <  2:719038/100‑1 (MQ=255)
                                      cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                            <  1:84106/100‑1 (MQ=255)
                                       aaaaGTTCTGAGCTAATTGCCGATACCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCTGAACCTGGCtc                                                           >  2:2069403/1‑100 (MQ=255)
                                       aaaaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtc                                                           >  1:1497003/1‑100 (MQ=255)
                                        aaaGTTCTTAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtca                                                          <  1:524081/100‑1 (MQ=255)
                                          aGTTCTGAGCTAATTGCCGATCACACTTGACGCAATCCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatc                                                        <  2:2453352/100‑1 (MQ=255)
                                          aGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatc                                                        <  2:2244966/100‑1 (MQ=255)
                                                gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                  >  2:1917339/1‑100 (MQ=255)
                                                    tAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTc                                              <  2:1700510/100‑1 (MQ=255)
                                                    tAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTc                                              <  1:507529/100‑1 (MQ=255)
                                                     aaTTGCCGATAAGATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                             <  2:2277190/100‑1 (MQ=255)
                                                     aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCGCGCTTTCGTACCCGAACCTGGCTCATCAGCGATTTTCt                                             >  1:2521723/1‑100 (MQ=255)
                                                      aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                            <  2:2290640/100‑1 (MQ=255)
                                                        tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                          <  2:2576469/100‑1 (MQ=255)
                                                        tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                          <  1:1232947/100‑1 (MQ=255)
                                                         gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                         <  2:1007159/100‑1 (MQ=255)
                                                          ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTc                                        <  1:1760465/100‑1 (MQ=255)
                                                           cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                       <  2:680409/100‑1 (MQ=255)
                                                            gATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCAt                                      <  1:586109/100‑1 (MQ=255)
                                                                ataTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCCTTCGTACCAGAACCTGGCTCATCTGTGGTTTTCTTTGTCATAATc                                  >  2:513470/3‑100 (MQ=255)
                                                                aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                  <  1:249322/100‑1 (MQ=255)
                                                                aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                  >  1:1725079/1‑100 (MQ=255)
                                                                 cATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCa                                 >  2:1278792/1‑100 (MQ=255)
                                                                  aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                                >  2:467113/1‑100 (MQ=255)
                                                                  aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                                <  1:1689581/100‑1 (MQ=255)
                                                                   tttGACGCAATGCGCACTAAAAGGGCATCATTTGCTGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAtt                               <  2:2082810/100‑1 (MQ=255)
                                                                      gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCt                            >  2:239973/1‑100 (MQ=255)
                                                                       aCGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                           <  2:945442/100‑1 (MQ=255)
                                                                       aCGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                           <  1:2385912/100‑1 (MQ=255)
                                                                       aCGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                           <  1:315639/100‑1 (MQ=255)
                                                                        cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                          >  1:1159509/1‑100 (MQ=255)
                                                                         gCAATGCGCACTAAAAGGGCATCATTTTATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATGAGTGATTTTCTTTGTCATAATCATTGCTgag                         <  2:143758/100‑1 (MQ=255)
                                                                         gCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgag                         >  2:512712/1‑100 (MQ=255)
                                                                           aaTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGAc                       >  2:128346/1‑100 (MQ=255)
                                                                           aaTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGAc                       <  2:348077/100‑1 (MQ=255)
                                                                                gtgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                  >  1:1323227/6‑100 (MQ=38)
                                                                                 tgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                 >  1:1906538/5‑100 (MQ=38)
                                                                                  gcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctct                >  1:647017/4‑100 (MQ=39)
                                                                                  gcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCGGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctct                >  2:2353354/4‑100 (MQ=38)
                                                                                   cTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCGCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg               >  1:2127367/1‑100 (MQ=255)
                                                                                    aaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACAAGAACCTGGCGCATCAGTGATTTTCTTTGTCATACTCATTGCTGAGACAGGCTCTGt              >  2:1170408/2‑100 (MQ=38)
                                                                                     aaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGtt             >  1:395548/1‑100 (MQ=255)
                                                                                      aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg            <  1:489604/100‑1 (MQ=255)
                                                                                      aaaGGGCATCATTTGATGCCCGTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg            <  1:1351657/100‑1 (MQ=255)
                                                                                       aaGGGCATCAAATGATGCCCTTTTTTTTTTTTTTTTTACCAGAACCTGGCTCATCAGTTATTTTCTTAGTCATAATCATTGCTGTGACAGGCCCTGTTGa           >  2:2597486/1‑100 (MQ=255)
                                                                                       aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa           >  1:1964608/1‑100 (MQ=255)
                                                                                        aGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg          <  1:1591340/100‑1 (MQ=255)
                                                                                        aGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg          >  2:1625319/1‑100 (MQ=38)
                                                                                         gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg         >  1:860108/1‑100 (MQ=255)
                                                                                            catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg      <  1:1667575/100‑1 (MQ=255)
                                                                                            catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg      <  2:1487280/100‑1 (MQ=255)
                                                                                              tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGta    <  1:362613/100‑1 (MQ=255)
                                                                                              tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGta    <  1:1840635/100‑1 (MQ=255)
                                                                                              tcatTTGATGCCCTTTATGCACGCTTTTGTACCAGAACGTGGCTCAACAATGATGTTCTTTGTCATAATCATTGCTTAGACGGGCTCTGTTGAGGGCGta    >  2:1155400/1‑100 (MQ=255)
                                                                                               catTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCGCATCAGTGATTTTCTTTGTCATAATCCTTGCTGAGACAGGCTGTGTTGAGGGCGtat   >  1:795694/1‑100 (MQ=255)
                                                                                                atTTGATGCCCTTTTTTTTTTTTTTTTTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGATACACTCTCTGTTGAGGGCGtata  >  2:1987993/1‑100 (MQ=255)
                                                                                                atTTGATGCCCTTTTTGCAAGCTTTCGTACCAGAACCTGGCTCATCAGTGCTTTTCTTTGTCATAATCATTGCTGAGCCAGGCTCTGTTGAGGGCGtata  >  2:697106/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATA  >  NC_000913/4177074‑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: ‑