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.

CGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATA  >  NC_000913/1576923‑1577356
                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
cGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaa                                                                                                                                                                                                      >  2:1118173/1‑238 (MQ=255)
 agtttctagATTAAAGGAAAGGTTAAAAAGTTTGTTTTAAAAAGCGGGGAGGTATAAGTTATTTCCGAAGGATCTTGGTTTAAAGGCGCTTTTGTTATTCCCACCAATTGGCGATGTTACCCCGCCCAACCGCGTGGTTCGATAACCAAATTCCCCCTTCCTATCATCGGCGTTGAAGGTTTGGATCTTTGAAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa                                                                                                                                                                                                     <  2:201446/230‑1 (MQ=255)
    ttCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaa                                                                                                                                                                                                      >  1:332368/1‑234 (MQ=255)
          aTTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGTCACAAAACTGACATaa                                                                                                                                                                                            >  2:62834/1‑238 (MQ=255)
              aaCGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCg                                                                                                                                                                                        >  1:713558/1‑238 (MQ=255)
                  aCAAGATATAAATGTTTTTTTTAAAATCTGGGGAGTAATAAGTCCTTCCAGAAGAGCTTTGATTAAATGCGCTCTTGTTATTCTCATCGATTTGCGATGTTACCCAGTCAGAACGCGTAGTGCGATACCCCAATTCCCCCTTGCTACCATCAGCGTTGACGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca                                                                                                                                                                                   <  2:292242/239‑1 (MQ=255)
                                ttCTTTTTAAAATCGGGGGTATAATAAGTATTTCCTAAAGGCCCTTTATTAAATGCTCCTTTTTTTTTCCCATCGATTTGCGATTTTAACCACCCCGGACGCGGATGTGAATAACCCAATTTCCCCTTGCTACCATCAGCGTTGAACCGTTTGATCTTTCCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAg                                                                                                                                                                       <  2:603489/237‑1 (MQ=255)
                                  cTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTaa                                                                                                                                                                        >  1:966079/1‑234 (MQ=255)
                                             cgtgggAGTAATAAGTCCTCCCTGAGGTGCCTTAATTAAATGCGTTTTTGTTATTCCCTTCGTTTGGGGTTGTTACCCCGGCCAACCTGTAATTGGAAAACCCCATTTTCCCCTTCCTATCATCAGCGTGGAAGCTTTTGATCTTGGCAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa                                                                                                                                                        <  2:374100/236‑1 (MQ=255)
                                                 ggtgTATTAAGTACTTCCTGAAGAGCTTTGATTAATTGCGCTTTTGTTATTCTCATCGATATGCGATGTAACCCAGTCCGACCGCTTGGTGCGATAACCCATTTCCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAc                                                                                                                                                     <  2:1179654/235‑1 (MQ=255)
                                                             tccTTCGTGAAGAGCCTGGATAAAATCCCTTCTTTTTATTCTCATCGATTGCGGTTGTTACCCCGTCCGACCGCGTGGTGCGAAACCCAAATTCCCCCTTCCTATCATAAGCGTTGAAGCGTTTGATTTTTGCAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATc                                                                                                                                         <  2:966079/237‑1 (MQ=255)
                                                                ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATccc                                                                                                                                       >  1:1138731/1‑238 (MQ=255)
                                                                ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATccc                                                                                                                                       >  1:1206655/1‑238 (MQ=255)
                                                                ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGACATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAAAAc                                                                                                                                                     >  1:1283308/1‑224 (MQ=255)
                                                                ttCCTAAAGAGCTTTGATTAAACGCCCTCTTGTTTTTCCCATCGTATGGGATTTTTAACCAGCCCGAACGGGTAGTGCGATAACCCAATTACCCCTTCCTATCATCAGCGTAGAAGCTTTTATTTTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCg                                                                                                                                      <  2:332368/239‑1 (MQ=255)
                                                                          tccTTGATAAAAGGCTCTCTTGTTATTCTCATCGATTGGCGATGTTAACCAGCCCGAACGCGTAGTGCTATCACCCAATTTCCCCTTGCTACCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAg                                                                                                                                                                       <  2:897502/195‑1 (MQ=255)
                                                                           tcTTGATTAAATGCCTTTTTTTTTTTTCCACCGATTTGCGATTTAAACCAGCCCGAACGCGTATTGCGAAACCCCAATTCCCCCTTCCTATCATCAGCGTTGAGGCTTTTGATCTTTGCAACATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca                                                                                                                                                                                   <  2:274462/182‑1 (MQ=255)
                                                                           ccTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGt                                                                                                                           >  1:59927/1‑239 (MQ=255)
                                                                           ccTTGATAAAATCCCTTCTTGTTATTCTCAGCGATAGGCGGTGTAAACCGGCCCGAACGCGTGGTGCAATACCCCAATTTCACCTTGTTATCATCGGCGTTGAAGCGTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGa                                                                                                                                                           <  2:719550/207‑1 (MQ=255)
                                                                             ttgtatAAATGCGTTCTTTGTATGCTCATCGATTTGCGATTTTAACCATTCCGAACGGGTGTGGCAATAACCCAATTTCCCCTTGCTAACATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTg                                                                                                                         <  2:59927/234‑1 (MQ=255)
                                                                                aTTAAATGCGTTCTTTTTATTCTCATCGATAGGCGATGTTACCCAGTCCGAACGCTGAGGGCGATAACCCAATTTCCCCATCCTACCACAAGCGTTAAAGGTTTTTATCTTGCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             <  2:1420500/200‑1 (MQ=255)
                                                                                   aattGCGCTTTTTTATTCTCCACCGATTGGCGTTGTTACCCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTCCTATCATCAGCGTTGAAGCGTTTGATCTTGCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactac                                                                                                                   <  2:234206/236‑1 (MQ=255)
                                                                                   aatggcgcTCTTGTTCTGCTCATGGATTTGCGATGTTACCAAGTCCGAGCGCGTAGTGCGATAACCCATTTCCCCCTGCCTACAACCAGGGTTAAACCTTTTGATTTCTACATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactac                                                                                                                   <  2:1126309/235‑1 (MQ=255)
                                                                                    aaTGCGCTTTGATTATTTCAATCAATTTGCGTTGTCAACCAGTCCAAACGCGTAGTGCAAAACCCAAATTCCCCCTTCCTATCTGCAGCGTTGAAGCTTTTGATCTTTCCAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactact                                                                                                                  <  2:1206655/239‑1 (MQ=255)
                                                                                    aaTGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             >  1:897502/1‑196 (MQ=255)
                                                                                     aTGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTa                                                                                                                                                                         >  1:274462/1‑183 (MQ=255)
                                                                                     aTGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCt                                                                                                                                                 >  1:719550/1‑207 (MQ=255)
                                                                                     aTGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCact                                                                                                                     >  1:457347/1‑235 (MQ=255)
                                                                                          ctctTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAccc                                                                                                                                                   >  1:1420500/1‑200 (MQ=255)
                                                                                          ctctTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCg                                                                                                                                      >  2:326264/1‑213 (MQ=255)
                                                                                          ctctTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCg                                                                                                                                      <  1:326264/213‑1 (MQ=255)
                                                                                                        cATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTc                                                                                                  >  2:456631/1‑236 (MQ=255)
                                                                                                        aatcaatgggCGATGTTAACCAGTCCGAACGCGTGGTGCGATACCCCATTTACCCCTGGCTACCATCGGCTGGAAGCTGTTTGACTTTTCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCt                                                                                               <  2:756580/230‑1 (MQ=255)
                                                                                                            cgattggagATGTCAACCAGCCCGAACGGGTGGTGCAATAACCCACTTTCCCCTTCCTACCACCGCCTTTAACGCTTTTTATCTTTGCACCATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTgg                                                                                                                            <  2:793142/199‑1 (MQ=255)
                                                                                                                    gATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCAtt                                                                                                                                                                                                                                 >  1:734496/1‑98 (MQ=255)
                                                                                                                       tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCactact                                                                                                                  >  1:793142/1‑207 (MQ=255)
                                                                                                                         ttACCCAGTCGGAGCGCATAGTGCGAAACCCCAATTTCCCCTCGCTACCATCAGCGTTGAAGCGTTTGATTTTAGCATCAATTACCCCGCCATTGAAGCGACCAAATACAAc                                                                                                                                                                                                               <  1:813097/112‑1 (MQ=255)
                                                                                                                          tAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAg                                                                               >  1:125797/1‑239 (MQ=255)
                                                                                                                           accattctcaaaCTCGTAGGGCGAACCCCCAATTCCAACTTGTTTACAACGGCGTCGACGTTTTAGATCTTTACATAAATAACCCCGCCATTGAAGCGACCAAATTCAACCGGAACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGt                                                                              <  2:640481/229‑1 (MQ=255)
                                                                                                                                 tCCGAACGGGTAGGGCAATACCCCAATTTCCCCTTTCTATCATAAGCGTGGAAGCGTTTGATCTTTACATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCa                                                                                                                       <  2:498145/192‑1 (MQ=255)
                                                                                                                                 tCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACt                                                                        >  1:803442/1‑239 (MQ=255)
                                                                                                                                   cGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa                                                                                                                                                                                                     >  2:813097/1‑112 (MQ=255)
                                                                                                                                           tAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc                                                                                                                                                                                                               >  2:1075543/1‑94 (MQ=255)
                                                                                                                                           tAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc                                                                                                                                                                                                               <  1:1075543/94‑1 (MQ=255)
                                                                                                                                           tAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCAt                                                                                                             >  1:498145/1‑192 (MQ=255)
                                                                                                                                           tAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATAAAGATAATACCCCTGTGCCATCCAGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGtca                                                              >  2:397114/1‑239 (MQ=255)
                                                                                                                                                 cGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATAGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAACACCGtccatcaa                                                         >  2:1050295/1‑239 (MQ=255)
                                                                                                                                                  gATAACCCATTTCCCCCTTCCTACCATCGGCGTTGAAGCGTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa                                                                                                                                                        <  2:1208604/143‑1 (MQ=255)
                                                                                                                                                   aTAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAAt                                                       >  1:363127/1‑239 (MQ=255)
                                                                                                                                                    tAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATa                                                      <  1:1118173/239‑1 (MQ=255)
                                                                                                                                                        cccAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTGGCACTACTGGCATCGGACTCATTCGCTgg                                                                                             >  1:1228939/1‑197 (MQ=255)
                                                                                                                                                         ccAAGTACACCTTGCTACCATCAGCGTTGAAGCGTTTGATTTTGTCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATGCCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTGATTAGTTGCACTAATACCGTCAATCAAATAGGCa                                                  <  2:521779/239‑1 (MQ=255)
                                                                                                                                                         caaaTTTCACCTGCCTATCATCAGCGGTGAAGCGTTTGATCTTGACATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCa                                                  <  2:125797/237‑1 (MQ=255)
                                                                                                                                                          cAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGttatta                                                                                >  1:10287/1‑208 (MQ=255)
                                                                                                                                                             tttCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCtgtg                                                                                                                                              >  1:1208604/1‑143 (MQ=255)
                                                                                                                                                                       cTATCATCAGCGTTGAAGCTTTTGATCTTTACATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACAt                                                                                                                                          <  2:639174/137‑1 (MQ=255)
                                                                                                                                                                            accaGCGTGAAAGCTTTTGATCTGTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCAc                               <  2:1138731/237‑1 (MQ=255)
                                                                                                                                                                                 cGTTGAGGCTTTTATTCTTAGCATAATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCa                                                  <  2:154544/215‑1 (MQ=255)
                                                                                                                                                                                 cGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGc                                                                                                                                                                      >  2:620485/1‑99 (MQ=255)
                                                                                                                                                                                 cGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGc                                                                                                                                                                      <  1:620485/99‑1 (MQ=255)
                                                                                                                                                                                 cGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAata                                                                                                                                >  1:639174/1‑137 (MQ=255)
                                                                                                                                                                                 cGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGa                          >  2:330818/1‑239 (MQ=255)
                                                                                                                                                                                    tGAAGCGTTTAATTTTTGACTCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATa                       <  2:1356212/239‑1 (MQ=255)
                                                                                                                                                                                      aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGca                                                                                                                                                                                                          <  2:967850/58‑1 (MQ=255)
                                                                                                                                                                                      aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGca                                                                                                                                                                                                          >  1:967850/1‑58 (MQ=255)
                                                                                                                                                                                      aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             >  2:787598/1‑103 (MQ=255)
                                                                                                                                                                                      aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                                                                                                                             <  1:787598/103‑1 (MQ=255)
                                                                                                                                                                                      aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATaga                     >  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)
                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
CGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGCTGGGTTCAGGTTATTAGTTGCACTAATACCGTCAATCAAATAGGCATTCTGGTAGGGCGACGCACCGTGAATAGAGATTTTCTCAGGGCGAATA  >  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: