Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,3970AC35.4% 211.0 / 25.7 136L198F (TTA→TTCyiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base A (58/26);  new base C (1/45);  total (65/71)
Fisher's exact test for biased strand distribution p-value = 3.74e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 138 total aligned reads displayed.

GGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCT  >  NC_000913/3745169‑3745629
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
ggTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTGGCATGATGATTGGGGCAAAGCTGATGGTTATCTGGGGGGGGGAGGCCAGTCGCCTTAATTTACCCCGCCAGc                                                                                                                                                                                                                                >  1:3547584/1‑239 (MQ=255)
    ccGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCCGGCATGATGATGGGCGCAACGCTGATGTTTACCTTGGGGGGGGAGGGCCGTCGCCTTAATTTATCCCGCCAGCaaaa                                                                                                                                                                                                                            >  1:244331/1‑239 (MQ=255)
    ccGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGCGCAAAGCTGATGCTTACCTGGTGGGGGGAGGCCCGTCGCCTTAATTTACCCCGCCCGCaaaa                                                                                                                                                                                                                            >  1:827839/1‑239 (MQ=255)
      gatgatGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTGGCCTGATGATGGGGGCAACGCTGATGCTTTCCTGGTGGTGGGCGGCCCGGCGGCTTAATTTTCCCCCCCACCAAAAAg                                                                                                                                                                                                                          >  1:2342016/1‑239 (MQ=255)
      gatgatGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCCTGATGGTTGGCGCAACGCTGATGGTTAACTTGTGGGGGCAGGCCAGGCGCCCTTATTTACCCCGCCCGCACAAAg                                                                                                                                                                                                                          >  2:582527/1‑239 (MQ=255)
           ggcgcAGCGCCACTTATCCGGTCACCCGGGGGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCCCCAATTATCCCCCCTTCCATTCCGTTTATTATCTTCGGCATTTCCAGCGGATTATACATCAGTAAGCTGTTTATGGCGGGCATTGCCCGTGGCATGAGGATGGGCCCAGCGCTGAGGCTTACCTGGGGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAAc                                                                                                                                                                                                                      <  2:2578995/237‑1 (MQ=255)
                  cgcCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATATCCCCCGGCGTGGAGATGGGCGCAACACTGATGCCTACCTGGTGGTGGCGGGCCAGGCGCCTTCAATTTTCCCGCCCCCAAAAAGCAACAATGCAga                                                                                                                                                                                                              >  2:3478804/1‑238 (MQ=255)
                         tacccGGTCACCCGGGGGGCGGGTCTGTTTGCCTCTGGCGGCATTATCGCGCCATTTATCCCCCCTTCCTTTCCTTTTATTTTCTTCGGGTTTTCCAGAGAATTATCCATCAGTAGGCTTTTTATGGCGGACATTACCCCGGGCATGATGATGGGCGCAACCCTGATGCTTGCCTGGTGATGGCAGGCAGGTCGCTTTAATTTACCCCGCCAGCAAAAAGCAACGATGCa                                                                                                                                                                                                                <  2:2114350/227‑1 (MQ=255)
                          aTCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTgg                                                                                                                                                                                                      <  1:3140186/239‑1 (MQ=255)
                                   aCCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCTTTTCCCCCCGGCTTCGGGTGGGGGCACAGCTCATTCTTACCCAGGGGGGGGAGGCCCGGCGCCTTCATTTATCCCCCCCGCACAAAAAAACGACGAGGGAAATCTGGCCCTCcttt                                                                                                                                                                                             >  2:2478836/1‑239 (MQ=255)
                                      gggCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtc                                                                                                                                                                                          <  1:456581/239‑1 (MQ=255)
                                                 cTGATTGCCTCTGGCGGCATTATCGCGCAAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCTTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTCCCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCgg                                                                                                                                                                                     <  1:2979763/233‑1 (MQ=255)
                                                     ttcccTCTGGGCGCATTTTCGCGCCAATTATCCCCCCTTCCATTCCTTTTATTATCTCCGGCTTTTCCAGCGGTTTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCACCGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                                           <  2:2911509/236‑1 (MQ=255)
                                                     tgtcttctGGCGGCATTATCGGCAAAATAATTCCCCTTTCCATCCCTTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGATTTTTATGGCGGGCATTGCCCCTGGCATAATGATGGGCGCACCGCTGTTGTTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgcg                                                                                                                                                                           <  2:244331/234‑1 (MQ=255)
                                                                  cATTATCGCGCCAATTATCCCCCCTTCCATTCCGTTTATTATCTTCGGCTTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCa                                                                                                                                                               <  1:2328277/238‑1 (MQ=255)
                                                                   aTTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCCCGCATGGTGGTGTGCGCAACACCGCTGCCTACCCCGGGGCGGCAGGTCAGTCGCCTTTCTTTTCCCCCCCAGCCAAAAACAACCATGCAGGAAAATCAGCGCTCTTCTTTCTCCCGTATCTGGGCGGTGtttttttccg                                                                                                                                                              >  2:3295195/1‑234 (MQ=255)
                                                                          cgcCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGAt                                                                                                                                                                                                                   <  2:1579155/178‑1 (MQ=255)
                                                                              aatttctCCCCCTCCGCTTCCTTTTATTATCTTGGGCTTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCGTTGCCCGTGGCATGATGAGGGGCGCACCGCTGATGCTTACCGGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTg                                                                                                                                                  <  2:3547584/233‑1 (MQ=255)
                                                                                ttattCCCCCTTCCATTCCTTTTATTATCTTCGGCGTTTCCAGGGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGc                                                                                                                                                <  1:3437727/239‑1 (MQ=255)
                                                                                      ccccTTCCATTCCGTTTATTATCTTCGGCTTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTCCCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGc                                                                                                                                                                                                                 <  1:2145314/168‑1 (MQ=255)
                                                                                         cTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGGGCAAAGCTGATGGTTTCCTGGGGGGGGCCGGGCAGGCGCCCTTATTTACCCCCCCAGCAAAAAACAACGGTGCAGGAAATCTGGCACTCTTTTGTCTCCGGGATCTGGGCGCCGGTTCCTTCAGTCATCATTATTTGCggggg                                                                                                                                           >  2:3437727/1‑233 (MQ=255)
                                                                                                tCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGGCTTTCCCCCGGCATGATGATGGGCGCCACGCTGATGCTTACCTGGTGGTGGGGGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAACAACGATGCAGGAAAACtt                                                                                                                                                                                                       >  2:2145314/1‑167 (MQ=255)
                                                                                                tCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTACTTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTTTCTCCGGTATATTGGCGCTGTTTCTTCCCGTCATCATTATTGGCGGGgtttttgtctgc                                                                                                                                 >  2:2107539/1‑229 (MQ=255)
                                                                                                    tttATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGGTTGTGGCGGGCATTGCCCCGGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGGGGCAGGCGAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCAc                                                                                                                                                                                                   <  1:636399/168‑1 (MQ=255)
                                                                                                              tCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCCTGATGATGGGGGCAACGCTGATGCTTACCTGGGGGGGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCCAAAAGCAACGATGCAGGGAAACTGGCACTCTTTTTTCTCCCGGATCTGGgc                                                                                                                                                                            >  2:552961/1‑181 (MQ=255)
                                                                                                              tCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCCCGATGATCGGCGCCACGCCGATGCTTGCCTGTGGGGGGCAGGCGACTCGCCCTACTATCACCCCCCCGCACAAAGCAAAGATGCGGGGAAACTCGTCCTCTTTTGtgt                                                                                                                                                                                         >  2:636399/1‑166 (MQ=255)
                                                                                                                               ttATCCTGCAGTAAGCTGTTGAGGGGGGACATTGCCCCTGGCATGATGAGGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGCCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTtgc                                                                                                 <  2:827839/239‑1 (MQ=255)
                                                                                                                                tATCCATCAGTAAGTTGTTTATGGCGGGCATTGCCCCGGGCATGATGATGGGCGCAACGCTGATGCTTCCCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTtgct                                                                                                <  1:571255/239‑1 (MQ=255)
                                                                                                                                  tCCATCAGTAAGCTTTTTATGGGGGGCATTGCCCCTGGCATGATGATGGGCGCACCGCTGATGCTTACCGGGTGGGGGCGGGCCAGTCGCTTTAATTTACCCCGCCAGCaaaa                                                                                                                                                                                                                            <  2:3172185/113‑1 (MQ=255)
                                                                                                                                       cAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAgg                                                                                                                                                                                                              <  1:1367871/122‑1 (MQ=255)
                                                                                                                                                 tttATGGCGGGCATTGCCCCTGGCATGATGATGGGGGCAACGCTGATGCTTACCTGGTGGGGGCAGGCCAGTCGCCTTAATTTAACCCGCCAGCAAAAAGCAACAATGCAGGAAAACTGGCa                                                                                                                                                                                                    >  2:1367871/1‑122 (MQ=255)
                                                                                                                                                 tttATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTCCCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCg                                                                                                                                        <  1:1027078/182‑1 (MQ=255)
                                                                                                                                                 tttATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCa                                                                                                                                                                                                    <  1:177858/122‑1 (MQ=255)
                                                                                                                                                 tttATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGtt                                                                                                   <  1:366383/219‑1 (MQ=255)
                                                                                                                                                 tttATGGCGGGCATTGCCCCTGGCATGATGAGGGGCGCACCGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACg                                                                                                                                                                                                                     <  1:2883732/105‑1 (MQ=255)
                                                                                                                                                    aTGGCGGGCATTGCCCCGGGCATGATGATGGGCGCAACGAGGATGTTTCCCTGGGGGTGGCAGGCAAGTCGCTTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGt                                                                            <  1:1067163/239‑1 (MQ=255)
                                                                                                                                                        cGGGCTTTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGc                                                                         <  1:582527/238‑1 (MQ=255)
                                                                                                                                                           tgaattccccTGGTCATGATGGTGGGCCAAACGTAGGTGTTTACCTGGTGGTGGAGGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCca                                                                       <  2:1468142/232‑1 (MQ=255)
                                                                                                                                                           gCATTTCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTACCCGGGGGGGGCAGGGCAGTCGCCTTAATTTACCCCCCCAGCCAAAAACAAAGAAGCAGGGAAACTTGCACTCTTTTGTCTCCCGGATCTTGgc                                                                                                                                                                            >  2:1296633/1‑136 (MQ=255)
                                                                                                                                                           gCATTGCCCCTGGCATGATGATGGGGGCAAAGCTGATGGTTACCTGGGGGGGGGGGGGCAGGCGCCTTTATTTACCCCCCCCGCAAAAAGCAACGATGCAGGAAAACTGGCACTCTTTTTtc                                                                                                                                                                                          >  2:177858/1‑122 (MQ=255)
                                                                                                                                                            catgcccccTGGCATGATGATGGGCGCAACGCTGATGCTTCCCCGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCcaca                                                                     <  1:1296633/233‑1 (MQ=255)
                                                                                                                                                            cATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCcaca                                                                     <  1:16422/238‑1 (MQ=255)
                                                                                                                                                                 ccccTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:2290306/1‑130 (MQ=255)
                                                                                                                                                                    cTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:1149653/1‑127 (MQ=255)
                                                                                                                                                                      ggCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtc                                                                                                                                                                                          >  1:1433221/1‑111 (MQ=255)
                                                                                                                                                                           tgatgatgGGCGCAACGCGGATCCTTCCCTGGGGGGGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCaaaa                                                                                                                                                                                                                            <  2:2413375/73‑1 (MQ=255)
                                                                                                                                                                            gatgatgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:3392890/1‑120 (MQ=255)
                                                                                                                                                                              tgatgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGtt                                                                                                                                           <  1:2585494/151‑1 (MQ=255)
                                                                                                                                                                              tgaggGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGc                                                                         <  1:1744335/217‑1 (MQ=255)
                                                                                                                                                                                 tgGGCGCACCGCGGATGCTTCCCGGGTGGTGGCGGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGTTGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGAc                                                 <  1:3115445/238‑1 (MQ=255)
                                                                                                                                                                                 tgGGCGCAACGCTGATGCTTACCTGGTGGGGGCAGGCCATTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGAc                                                 <  1:1759549/238‑1 (MQ=255)
                                                                                                                                                                                  gggCGCAACGCTGATGCTTACCTGGTGGGGGCAGGCCAGTCGCCTTACTTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:1881039/1‑114 (MQ=255)
                                                                                                                                                                                     cgcAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATg                                                                                                                                                                                                                  >  1:2413375/1‑73 (MQ=255)
                                                                                                                                                                                       cAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGaa                                                                                                                                                                                                            >  2:2919984/1‑77 (MQ=255)
                                                                                                                                                                                       cAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGaa                                                                                                                                                                                                            <  1:2919984/77‑1 (MQ=255)
                                                                                                                                                                                        aaCGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCtttt                                                                                                                                                                                             >  1:1552393/1‑91 (MQ=255)
                                                                                                                                                                                        aaCGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCtttt                                                                                                                                                                                             <  2:1552393/91‑1 (MQ=255)
                                                                                                                                                                                        aaCGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:1744335/1‑108 (MQ=255)
                                                                                                                                                                                        aaCGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGcc                                          <  1:3543017/238‑1 (MQ=255)
                                                                                                                                                                                          cGCTTATCCTTCCCGGGTGGGGGCAGCCCATTCGTCTTATTTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCcaca                                       <  1:2981743/239‑1 (MQ=255)
                                                                                                                                                                                            gtGATGCTGACCGGGTGGTGGCAGGCAAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACt                                     <  1:1926979/238‑1 (MQ=255)
                                                                                                                                                                                            cTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCatta                                                                                                                                                     <  2:1149653/127‑1 (MQ=255)
                                                                                                                                                                                             tGATGCTTCCCTGGTGTGGGCAGGCCAGTCGCTTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACt                                     <  1:62419/238‑1 (MQ=255)
                                                                                                                                                                                             tGATGCTTCCCTGGTGGGGGCGGCCCAGCCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACt                                     <  1:1766783/238‑1 (MQ=255)
                                                                                                                                                                                             tGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACt                                     <  1:552961/238‑1 (MQ=255)
                                                                                                                                                                                                 gtttccttgttggtggCGGGCAGGCCGCTTTATTTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGc                                <  1:1570164/229‑1 (MQ=255)
                                                                                                                                                                                                 gCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGc                                <  1:2478836/239‑1 (MQ=255)
                                                                                                                                                                                                 acttcctgggtggtggCCGGCCATCCCCTTAAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGc                                <  2:590718/231‑1 (MQ=255)
                                                                                                                                                                                                   ttACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCca                                                                       <  1:2225735/198‑1 (MQ=255)
                                                                                                                                                                                                    tACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:3547900/1‑96 (MQ=255)
                                                                                                                                                                                                     aCCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:94898/1‑95 (MQ=255)
                                                                                                                                                                                                     aCCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:2962337/1‑95 (MQ=255)
                                                                                                                                                                                                         ggtggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:1678937/1‑91 (MQ=255)
                                                                                                                                                                                                         ggtggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:1279104/1‑91 (MQ=255)
                                                                                                                                                                                                         ggtggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:26670/1‑91 (MQ=255)
                                                                                                                                                                                                          gtggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:210062/1‑90 (MQ=255)
                                                                                                                                                                                                           tggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:106382/1‑89 (MQ=255)
                                                                                                                                                                                                           tggtggCAGGCCAGTCCCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCgg                                                                                                                                  <  1:3435966/131‑1 (MQ=255)
                                                                                                                                                                                                            ggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:3094216/1‑88 (MQ=255)
                                                                                                                                                                                                             gtggCAGGCCAGTCGCTTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCca                                                                       <  1:188316/188‑1 (MQ=255)
                                                                                                                                                                                                              tggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:405949/1‑86 (MQ=255)
                                                                                                                                                                                                               ggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:119672/1‑85 (MQ=255)
                                                                                                                                                                                                                gCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:508614/1‑84 (MQ=255)
                                                                                                                                                                                                                   ggCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGaa                                                                                                                                                                                                            >  2:777343/1‑49 (MQ=255)
                                                                                                                                                                                                                   ggCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGaa                                                                                                                                                                                                            <  1:777343/49‑1 (MQ=255)
                                                                                                                                                                                                                    gCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGtt                                                                                                                            <  1:1394796/128‑1 (MQ=255)
                                                                                                                                                                                                                    gCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:2459291/1‑80 (MQ=255)
                                                                                                                                                                                                                    gCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:2867394/1‑80 (MQ=255)
                                                                                                                                                                                                                    gCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:570834/1‑80 (MQ=255)
                                                                                                                                                                                                                    gCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            <  1:1200453/80‑1 (MQ=255)
                                                                                                                                                                                                                    gCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:1200453/1‑80 (MQ=255)
                                                                                                                                                                                                                     ccAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:3435966/1‑79 (MQ=255)
                                                                                                                                                                                                                     ccAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:2217894/1‑79 (MQ=255)
                                                                                                                                                                                                                       aGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:188316/1‑77 (MQ=255)
                                                                                                                                                                                                                          ccccTTAAGTCACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAac       <  1:2146646/237‑1 (MQ=255)
                                                                                                                                                                                                                          cGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAac       <  1:1519707/239‑1 (MQ=255)
                                                                                                                                                                                                                              ttAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  2:1394796/1‑70 (MQ=255)
                                                                                                                                                                                                                               tAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGTAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCt  <  1:2681211/239‑1 (MQ=255)
                                                                                                                                                                                                                                 aaTTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGgc                                                                                                                                                                            >  1:3251553/1‑68 (MQ=255)
                                                                                                                                                                                                                                  aTTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGc                                                                                                       <  1:3256443/136‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
GGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCCCCCTTCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCT  >  NC_000913/3745169‑3745629

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