Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,775,4670AC27.4% 130.7 / 12.5 62intergenic (+23/‑120)aroD/ydiF3‑dehydroquinate dehydratase/putative acetate‑CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (38/7);  new base C (0/17);  total (38/24)
Fisher's exact test for biased strand distribution p-value = 4.68e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAAC  >  NC_000913/1775333‑1775602
                                                                                                                                      |                                                                                                                                       
tGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTc                                                                                                                                    >  1:396964/1‑140 (MQ=255)
     tATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCccgc                                                                                                                               >  1:220121/1‑140 (MQ=255)
       tttGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCccgccg                                                                                                                             >  1:458606/1‑140 (MQ=255)
               ggcggcAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaa                                                                                                                       >  1:536269/1‑138 (MQ=255)
               ggcggcAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaa                                                                                                                       <  2:536269/138‑1 (MQ=255)
                 cggcAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                                                                                                                   >  2:1013339/1‑140 (MQ=255)
                 cggcAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                                                                                                                   >  1:1004413/1‑140 (MQ=255)
                 cggcAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                                                                                                                   >  2:106873/1‑140 (MQ=255)
                 cggcAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                                                                                                                   >  1:531407/1‑140 (MQ=255)
                     aaCTTTTGGTGCGGTAAAAAAAGCGTTTGCGCCAGGGAAAATCTCGGTAAATGATTTGCGCCCGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAatata                                                                                                               <  2:1308230/140‑1 (MQ=255)
                          ttGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCaa                                                                                                          <  2:1170432/140‑1 (MQ=255)
                          ttGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCaa                                                                                                          <  2:695948/140‑1 (MQ=255)
                                        aaaGCTTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATaaa                                                                                            <  2:76581/140‑1 (MQ=255)
                                             gTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAatat                                                                                       >  1:1306460/1‑140 (MQ=255)
                                             gTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAatat                                                                                       >  1:414997/1‑140 (MQ=255)
                                             gTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAatat                                                                                       >  2:698854/1‑140 (MQ=255)
                                              tCTGCGCAGGGGAAAACCTCGGTAAATGATTTGCGCACGGTATAACCTATTTCACACCAGCCATAAGAAATAATATTCCGCGGGAAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATAtt                                                                                      <  2:403522/140‑1 (MQ=255)
                                                tGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACCCCAGGCATAAGCAATAATATTTCGGCGGAAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTc                                                                                    <  1:558467/140‑1 (MQ=255)
                                                           aaaTCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCtg                                                                         <  1:1068264/140‑1 (MQ=255)
                                                           aaaTCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCtg                                                                         <  2:1004413/140‑1 (MQ=255)
                                                                  gaaaaaTGATTTGCGCTCGTTATAAACAATTATCCAACAGATATAAGTATTAATATTTCGGCGAAATCACCAACCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAATCTCTGTGCGCTa                                                                  <  2:1306460/137‑1 (MQ=255)
                                                                     aagtattttGCGCAGGGTATTATTTCTTTTACCCCAGGCATAAGCAATACTTTTTTGGCGGTATCCCCCTCCCCGCCGATCTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAAtt                                                               <  2:685210/134‑1 (MQ=255)
                                                                      aaTGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTGGGCGGGACCCCCCCCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATtc                                                              <  2:523624/140‑1 (MQ=255)
                                                                       aTGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATtct                                                             <  1:2557/140‑1 (MQ=255)
                                                                       aTGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGCATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATtct                                                             <  2:407412/140‑1 (MQ=255)
                                                                       aTGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGGACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTTGTGAGTCTCCGTGCGCTAATtct                                                             >  1:1081986/1‑140 (MQ=255)
                                                                         gATTTGCGCACGGTATTAACTATTTTACACCAGGCTTAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTcc                                                           >  2:204793/1‑140 (MQ=255)
                                                                            ttGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCAtt                                                        <  1:653605/140‑1 (MQ=255)
                                                                              gcgcACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTg                                                      <  2:531407/140‑1 (MQ=255)
                                                                                 cACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCg                                                   <  1:793498/140‑1 (MQ=255)
                                                                                 cACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCg                                                   >  2:779944/1‑140 (MQ=255)
                                                                                      tATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAggg                                              <  2:935391/140‑1 (MQ=255)
                                                                                      tATTAACTATTTTACACCAGGAATAAGCAATAATATTTCGGGGGGAACCCCCCCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAggg                                              <  2:48850/140‑1 (MQ=255)
                                                                                         tAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGaaa                                           >  2:1446996/1‑140 (MQ=255)
                                                                                            cTATTTTACACCAGGCTTAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATc                                        <  1:204793/140‑1 (MQ=255)
                                                                                              aTTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATcac                                      >  1:1206571/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAAt                                         >  1:876994/1‑135 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAAt                                         <  2:876994/135‑1 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  2:750861/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  2:815750/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  2:1060466/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  2:1253298/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  2:1222816/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  2:117718/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  2:1175529/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  1:899761/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  1:272508/1‑140 (MQ=255)
                                                                                                tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACAt                                    >  1:493007/1‑140 (MQ=255)
                                                                                                  tACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCt                                  >  1:706339/1‑140 (MQ=255)
                                                                                                    cacCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGa                                >  2:1100853/1‑140 (MQ=255)
                                                                                                      ccAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAAt                              <  1:1222816/140‑1 (MQ=255)
                                                                                                        aGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCa                            >  1:1373033/1‑140 (MQ=255)
                                                                                                         ggcgtAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAg                           >  2:1133194/5‑140 (MQ=255)
                                                                                                         ggCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAg                           >  1:428133/1‑140 (MQ=255)
                                                                                                         ggCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATAAATTTAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAg                           >  2:572323/1‑140 (MQ=255)
                                                                                                         ggCATAAGCAATAATATTTCGGCGGGAACACCCCCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAAACACATCTGAATCAg                           >  2:1482894/1‑140 (MQ=255)
                                                                                                               aGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTa                     <  2:1206571/140‑1 (MQ=255)
                                                                                                               aGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTa                     <  2:706339/140‑1 (MQ=255)
                                                                                                               aGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTa                     <  1:539739/140‑1 (MQ=255)
                                                                                                                    taataTTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAAt                >  1:362346/1‑140 (MQ=255)
                                                                                                                    taataTTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAAt                >  1:121633/1‑140 (MQ=255)
                                                                                                                     aataTTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATg               >  1:590037/1‑140 (MQ=255)
                                                                                                                       tatTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGaa             <  2:705502/140‑1 (MQ=255)
                                                                                                                       tatTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGaa             <  1:1175529/140‑1 (MQ=255)
                                                                                                                          ttCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGAAAcc          <  1:750861/140‑1 (MQ=255)
                                                                                                                          ttCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGAAAcc          >  2:1201262/1‑140 (MQ=255)
                                                                                                                            cggcggGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGAAACCTg        <  1:1253298/140‑1 (MQ=255)
                                                                                                                                  gAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAac  <  1:1013339/140‑1 (MQ=255)
                                                                                                                                  gAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAac  <  2:458606/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAAC  >  NC_000913/1775333‑1775602

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: