Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,577,1230AC26.9% 192.4 / 10.4 78V166G (GTA→GGA) yddBputative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base A (42/15);  new base C (0/21);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.02e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAG‑TACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATA  >  NC_000913/1576923‑1577356
                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
cGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAG‑TACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaa                                                                                                                                                                                                      >  2:1118173/1‑238 (MQ=255)
 agtttctagATTAAAGGA‑AAGGTTAAAAAGTTTGTTTTAAAAAGCGGGGAGGTATAAGTTATTTCCGAAGGATCTTGGTTTAAAGGCGCTTTTGTTATTCC‑CACCAATTGGC‑GA‑TGTTACCCCGCCCAACCGCGTG‑GTTCGATAA‑CCAAATTCCCCCTTCCTATCATCGGCGTTGAAGG‑TTTGGATCTTTGAAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa                                                                                                                                                                                                     <  2:201446/230‑1 (MQ=255)
    ttCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAG‑TACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaa                                                                                                                                                                                                      >  1:332368/1‑234 (MQ=255)
          aTTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAG‑TACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGTCACAAAACTGACATaa                                                                                                                                                                                            >  2:62834/1‑238 (MQ=255)
              aaCGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAG‑TACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCg                                                                                                                                                                                        >  1:713558/1‑238 (MQ=255)
                  aCAAGATATAAATGTTTTTTTTAAAATCTGGGGAGTAATAAG‑TCCTTCCAGAAGAGCTTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATTTGC‑GA‑TGTTACCCAGTCAGAACGCGTA‑GTGCGATAC‑CCCAATTCCCCCTTGCTACCATCAGCGTTGACGC‑TTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca                                                                                                                                                                                   <  2:292242/239‑1 (MQ=255)
                                ttCTTTTTAAAATCGGGGGTATAATAAG‑TATTTCCTAAAGGCCCTTTATTAAATGCTCCTTTTTTTTTCC‑CATCGATTTGC‑GA‑TTTTAACCACCCCGGACGCGGATGTG‑AATAA‑CCCAATTTCCCCTTGCTACCATCAGCGTTGAACC‑GTTTGATCTTTCCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAg                                                                                                                                                                       <  2:603489/237‑1 (MQ=255)
                                  cTTTTTAAAATCTGGGGAGTAATAAG‑TACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTaa                                                                                                                                                                        >  1:966079/1‑234 (MQ=255)
                                             cgtgggAGTAATAAG‑TCCTCCCTGAGGTGCCTTAATTAAATGCGTTTTTGTTATTCC‑CTTCGTTTGGG‑GT‑TGTTACCCCGGCC‑AACCTGTA‑ATTGGAAAACCCCATTTTCCCCTTCCTATCATCAGCGTGGAAGC‑TTTTGATCTTGGCAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa                                                                                                                                                        <  2:374100/236‑1 (MQ=255)
                                                 ggtgTATTAAG‑TACTTCCTGAAGAGCTTTGATTAATTGCGCTTTTGTTATTCT‑CATCGATATGC‑GA‑TGTAACCCAGTCCGACCGCTTG‑GTGCGATAA‑CCCATTTCCCCCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAc                                                                                                                                                     <  2:1179654/235‑1 (MQ=255)
                                                             tccTTCGTGAAGAGCCTGGATAAAATCCCTTCTTTTTATTCT‑CATCGAT‑TGCGGT‑TGTTACCCCGTCCGACCGCGTG‑GTGCGAAAC‑CCAAATTCCCCCTTCCTATCATAAGCGTTGAAGC‑GTTTGATTTTTGCAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATc                                                                                                                                         <  2:966079/237‑1 (MQ=255)
                                                                ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATccc                                                                                                                                       >  1:1138731/1‑238 (MQ=255)
                                                                ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATccc                                                                                                                                       >  1:1206655/1‑238 (MQ=255)
                                                                ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGACATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAAAAc                                                                                                                                                     >  1:1283308/1‑224 (MQ=255)
                                                                ttCCTAAAGAGCTTTGATTAAACGCCCTCTTGTTTTTCC‑CATCG‑TATGG‑GATTTTTAACCAGCCCGAACGGGTA‑GTGCGATAA‑CCCAATTACCCCTTCCTATCATCAGCGTAGAAGC‑TTTTATTTTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCg                                                                                                                                      <  2:332368/239‑1 (MQ=255)
                                                                          tccTTGATAAAAGGCTCTCTTGTTATTCT‑CATCGATTGGC‑GA‑TGTTAACCAGCCCGAACGCGTA‑GTGCTATCA‑CCCAATTTCCCCTTGCTACCATCAGCGTTGAAGC‑TTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAg                                                                                                                                                                       <  2:897502/195‑1 (MQ=255)
                                                                           tcTTGATTAAATGCCTTTTTTTTTTTTC‑CACCGATTTGC‑GA‑TTTAAACCAGCCCGAACGCGTA‑TTGCGAAAC‑CCCAATTCCCCCTTCCTATCATCAGCGTTGAGGC‑TTTTGATCTTTGCAACATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca                                                                                                                                                                                   <  2:274462/182‑1 (MQ=255)
                                                                           ccTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGt                                                                                                                           >  1:59927/1‑239 (MQ=255)
                                                                           ccTTGATAAAATCCCTTCTTGTTATTCT‑CAGCGATAGGC‑GG‑TGTAAACCGGCCCGAACGCGTG‑GTGCAATAC‑CCCAATTTCACCTTGTTATCATCGGCGTTGAAGC‑GTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGa                                                                                                                                                           <  2:719550/207‑1 (MQ=255)
                                                                             ttgtatAAATGCGTTCTTTGTATGCT‑CATCGATTTGC‑GA‑TTTTAACCATTCCGAACGGGTG‑TGGCAATAA‑CCCAATTTCCCCTTGCTAACATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTg                                                                                                                         <  2:59927/234‑1 (MQ=255)
                                                                                aTTAAATGCGTTCTTTTTATTCT‑CATCGATAGGC‑GA‑TGTTACCCAGTCCGAACGCTGA‑GGGCGATAA‑CCCAATTTCCCCATCCTACCACAAGCGTTAAAGG‑TTTTTATCTTGCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             <  2:1420500/200‑1 (MQ=255)
                                                                                   aattGCGCT‑TTTTTATTCTCCACCGATTGGC‑GT‑TGTTACCCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTCCTATCATCAGCGTTGAAGC‑GTTTGATCTTGCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactac                                                                                                                   <  2:234206/236‑1 (MQ=255)
                                                                                   aatggcgcTCTTGTTCTGCT‑CATGGATTTGC‑GA‑TGTTACCAAGTCCGAGCGCGTA‑GTGCGATAA‑CCCATTTCCCCCTGCCTACAACCAGGGTTAAACC‑TTTTGATTTCTACATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactac                                                                                                                   <  2:1126309/235‑1 (MQ=255)
                                                                                    aaTGCGCTTTGATTATTTC‑AATCAATTTGC‑GT‑TGTCAACCAGTCCAAACGCGTA‑GTGCAAAAC‑CCAAATTCCCCCTTCCTATCTGCAGCGTTGAAGC‑TTTTGATCTTTCCAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactact                                                                                                                  <  2:1206655/239‑1 (MQ=255)
                                                                                    aaTGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             >  1:897502/1‑196 (MQ=255)
                                                                                     aTGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTa                                                                                                                                                                         >  1:274462/1‑183 (MQ=255)
                                                                                     aTGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCt                                                                                                                                                 >  1:719550/1‑207 (MQ=255)
                                                                                     aTGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCact                                                                                                                     >  1:457347/1‑235 (MQ=255)
                                                                                          ctctTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAccc                                                                                                                                                   >  1:1420500/1‑200 (MQ=255)
                                                                                          ctctTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCg                                                                                                                                      >  2:326264/1‑213 (MQ=255)
                                                                                          ctctTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCg                                                                                                                                      <  1:326264/213‑1 (MQ=255)
                                                                                                        cATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTc                                                                                                  >  2:456631/1‑236 (MQ=255)
                                                                                                        aatcaatgggC‑GA‑TGTTAACCAGTCCGAACGCGTG‑GTGCGATAC‑CCCATTTACCCCTGGCTACCATCGGC‑TGGAAGCTGTTTGACTTTTCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCt                                                                                               <  2:756580/230‑1 (MQ=255)
                                                                                                            cgattggagA‑TGTCAACCAGCCCGAACGGGTG‑GTGCAATAA‑CCCACTTTCCCCTTCCTACCACCGCCTTTAACGC‑TTTTTATCTTTGCACCATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTgg                                                                                                                            <  2:793142/199‑1 (MQ=255)
                                                                                                                    gA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCAtt                                                                                                                                                                                                                                 >  1:734496/1‑98 (MQ=255)
                                                                                                                       tGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactact                                                                                                                  >  1:793142/1‑207 (MQ=255)
                                                                                                                         ttACCCAGTCGGAGCGCATA‑GTGCGAAAC‑CCCAATTTCCCCTCGCTACCATCAGCGTTGAAGC‑GTTTGATTTTAGCATCAATTACCCCGCCATTGAAGCGACCAAATACAAc                                                                                                                                                                                                               <  1:813097/112‑1 (MQ=255)
                                                                                                                          tAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAg                                                                               >  1:125797/1‑239 (MQ=255)
                                                                                                                           accattctcaaaCTCGTA‑GGGCGAACC‑CCCAATTCCAACTTGTTTACAACGGCGTCGACGT‑TTTAGATCTTTACATAAATAACCCCGCCATTGAAGCGACCAAATTCAACCGGAACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGt                                                                              <  2:640481/229‑1 (MQ=255)
                                                                                                                                 tCCGAACGGGTA‑GGGCAATAC‑CCCAATTTCCCCTTTCTATCATAAGCGTGGAAGC‑GTTTGATCTTTACATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCa                                                                                                                       <  2:498145/192‑1 (MQ=255)
                                                                                                                                 tCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACt                                                                        >  1:803442/1‑239 (MQ=255)
                                                                                                                                   cGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa                                                                                                                                                                                                     >  2:813097/1‑112 (MQ=255)
                                                                                                                                           tA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc                                                                                                                                                                                                               >  2:1075543/1‑94 (MQ=255)
                                                                                                                                           tA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc                                                                                                                                                                                                               <  1:1075543/94‑1 (MQ=255)
                                                                                                                                           tA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCAt                                                                                                             >  1:498145/1‑192 (MQ=255)
                                                                                                                                           tA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATAAAGATAATACCCCTGTGCCATCCAGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGtca                                                              >  2:397114/1‑239 (MQ=255)
                                                                                                                                                 cGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATAGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAACACCGtccatcaa                                                         >  2:1050295/1‑239 (MQ=255)
                                                                                                                                                  gATAA‑CCCATTTCCCCCTTCCTACCATCGGCGTTGAAGC‑GTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa                                                                                                                                                        <  2:1208604/143‑1 (MQ=255)
                                                                                                                                                   aTAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAAt                                                       >  1:363127/1‑239 (MQ=255)
                                                                                                                                                    tAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATa                                                      <  1:1118173/239‑1 (MQ=255)
                                                                                                                                                        cccAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTGGCACTACTGGCATCGGACTCATTCGCTgg                                                                                             >  1:1228939/1‑197 (MQ=255)
                                                                                                                                                         ccAAGTACACCTTGCTACCATCAGCGTTGAAGC‑GTTTGATTTTGTCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATGCCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTGATTAGTTGCACTAATACCGTCAATCAAATAGGCa                                                  <  2:521779/239‑1 (MQ=255)
                                                                                                                                                         caaaTTTCACCTGCCTATCATCAGCGGTGAAGC‑GTTTGATCTTGACATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCa                                                  <  2:125797/237‑1 (MQ=255)
                                                                                                                                                          cAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGttatta                                                                                >  1:10287/1‑208 (MQ=255)
                                                                                                                                                             tttCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCtgtg                                                                                                                                              >  1:1208604/1‑143 (MQ=255)
                                                                                                                                                                       cTATCATCAGCGTTGAAGC‑TTTTGATCTTTACATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACAt                                                                                                                                          <  2:639174/137‑1 (MQ=255)
                                                                                                                                                                            accaGCGTGAAAGC‑TTTTGATCTGTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCAc                               <  2:1138731/237‑1 (MQ=255)
                                                                                                                                                                                 cGTTGAGGC‑TTTTATTCTTAGCATAATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCa                                                  <  2:154544/215‑1 (MQ=255)
                                                                                                                                                                                 cGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGc                                                                                                                                                                      >  2:620485/1‑99 (MQ=255)
                                                                                                                                                                                 cGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGc                                                                                                                                                                      <  1:620485/99‑1 (MQ=255)
                                                                                                                                                                                 cGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAata                                                                                                                                >  1:639174/1‑137 (MQ=255)
                                                                                                                                                                                 cGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGa                          >  2:330818/1‑239 (MQ=255)
                                                                                                                                                                                    tGAAGC‑GTTTAATTTTTGACTCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATa                       <  2:1356212/239‑1 (MQ=255)
                                                                                                                                                                                      aaGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGca                                                                                                                                                                                                          <  2:967850/58‑1 (MQ=255)
                                                                                                                                                                                      aaGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGca                                                                                                                                                                                                          >  1:967850/1‑58 (MQ=255)
                                                                                                                                                                                      aaGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             >  2:787598/1‑103 (MQ=255)
                                                                                                                                                                                      aaGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             <  1:787598/103‑1 (MQ=255)
                                                                                                                                                                                      aaGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATaga                     >  1:225638/1‑239 (MQ=255)
                                                                                                                                                                                            tttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc                                                                                                                                                                                                               <  1:336157/48‑1 (MQ=255)
                                                                                                                                                                                            tttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc                                                                                                                                                                                                               >  2:336157/1‑48 (MQ=255)
                                                                                                                                                                                            tttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTTGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGAttt                >  2:270348/1‑239 (MQ=255)
                                                                                                                                                                                            tttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAgg                                        >  1:154544/1‑215 (MQ=255)
                                                                                                                                                                                            tttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGAttt                >  1:535841/1‑239 (MQ=255)
                                                                                                                                                                                               gATCTTTGCATCTATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGt                                                                >  1:662389/1‑188 (MQ=255)
                                                                                                                                                                                                atttttCCATCAATCCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTtctc            <  2:225638/236‑1 (MQ=255)
                                                                                                                                                                                                   tttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATa                                                                                                                                                                                             >  1:501426/1‑59 (MQ=255)
                                                                                                                                                                                                   tttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATa                                                                                                                                                                                             <  2:501426/59‑1 (MQ=255)
                                                                                                                                                                                                           aaTTACCCCGCCATTGAAGCGCCCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCCCATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGaaca  >  1:1048864/1‑236 (MQ=255)
                                                                                                                                                                                                               aCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             >  2:193424/1‑79 (MQ=255)
                                                                                                                                                                                                               aCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             <  1:193424/79‑1 (MQ=255)
                                                                                                                                                                                                               aCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa                                                                                                                                                        >  2:1090285/1‑84 (MQ=255)
                                                                                                                                                                                                               aCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa                                                                                                                                                        <  1:1090285/84‑1 (MQ=255)
                                                                                                                                                                                                               aCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCa                                                  <  2:1360374/186‑1 (MQ=255)
                                                                                                                                                                                                               aCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCa                                                  >  1:1360374/1‑186 (MQ=255)
                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
CGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAG‑TACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCT‑CATCGATATGC‑GA‑TGTTAACCAGTCCGAACGCGTA‑GTGCGATAA‑CCCAATTTCACCTTGCTATCATCAGCGTTGAAGC‑GTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATA  >  NC_000913/1576923‑1577356

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑