Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,851,6540TG28.4% 175.4 / 10.2 75G218G (GGT→GGGhypBhydrogenase isoenzymes nickel incorporation protein HypB
Reads supporting (aligned to +/- strand):  ref base T (11/42);  new base G (21/0);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 8.52e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGAAGGCGACCAGCAAACCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGAC  >  NC_000913/2851409‑2851902
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
aTTGAAGGCCACCAGCACACCGGGAACGATGCCGCTCGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGACCACCGGTAAAGGCTGCCATCTTGACGCACAGATGAGTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCGCGTATGCCCGGCCAGCTGCGATCTCGGGGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAg                                                                                                                                                                                                                                                         <  2:304187/250‑1 (MQ=255)
    aaGGCGACCAGCAAACCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACaaa                                                                                                                                                                                                                                                    <  1:155131/251‑1 (MQ=255)
    aaGGCAACCAGCAAACCGTTAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACaaa                                                                                                                                                                                                                                                    <  1:355830/251‑1 (MQ=255)
       gCGACCAGCAAACCGTGAACGATGCCGCACGCATTCTCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCa                                                                                                                                                                                                                                                 <  2:331945/251‑1 (MQ=255)
                aaaCCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAatat                                                                                                                                                                                                                                        <  1:338031/251‑1 (MQ=255)
                  aCCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATcc                                                                                                                                                                                                                                      <  1:246344/251‑1 (MQ=255)
                     gTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTTCCTTACCGACGGGGAAAACAAACCACTGAAATATCCGc                                                                                                                                                                                                                                    >  1:76097/1‑250 (MQ=255)
                      tGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTTCCGAAGGGGAAGACAAACCACTGAAATATCCGCat                                                                                                                                                                                                                                  >  1:27737/1‑251 (MQ=255)
                       caaCGATGCCTCACGCATTCGCGCTACCGGCCCACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCat                                                                                                                                                                                                                                  <  1:243298/249‑1 (MQ=255)
                         aCGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATg                                                                                                                                                                                                                               <  1:292519/251‑1 (MQ=255)
                         aCGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCCTTTCCGCAGGGGGAGACAAAACACTGAAATATcaccaaaat                                                                                                                                                                                                                               >  2:437699/1‑243 (MQ=255)
                          cGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGCGCTTTCCGTTACCGAAGGTGAAGACAAACCACGGAAAAATCCGCCTATGt                                                                                                                                                                                                                              >  1:560109/1‑251 (MQ=255)
                           gATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGtt                                                                                                                                                                                                                             <  2:76097/251‑1 (MQ=255)
                            attccGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGttt                                                                                                                                                                                                                            <  2:6007/248‑1 (MQ=255)
                            aTGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGATTTTCGTTACCGAAGGGGGAGACaaaaaaacgaaaaatcaccatatgttt                                                                                                                                                                                                                            >  2:240223/1‑227 (MQ=255)
                                     gcaTTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTTCCGCAGGGGAAGACAAACCACTCAAatataccccaatgttttccgccccc                                                                                                                                                                                                                   >  2:95877/1‑230 (MQ=255)
                                          cgcgCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCt                                                                                                                                                                                                              <  1:240223/251‑1 (MQ=255)
                                             gcTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGAt                                                                                                                                                                                                           <  1:377794/251‑1 (MQ=255)
                                                    gCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGTCGCACAGATGATTGCTGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTc                                                                                                                                                                                                    <  1:328435/251‑1 (MQ=255)
                                                      acacCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGGTTTCCCGTACCCAAGGGGGAGAAAAAACACCGAAAAATCCGCATATGTTTTCCGCCCCCCCGCTGAtgcttcgccc                                                                                                                                                                                                  >  2:107135/1‑245 (MQ=255)
                                                            gCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGt                                                                                                                                                                                            <  1:173164/251‑1 (MQ=255)
                                                             cGATTCAGGTGAACACCGGTAAAGGCTTCCATCTTGACGCACAGATGATTGCCGACGCCGCTCCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGtt                                                                                                                                                                                           <  1:151911/251‑1 (MQ=255)
                                                                  cAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGGGAAAACAAACCACCGAAAAATCCGCATTTGTTTGCCGCCCCCCCGCTGATGCTGCTCCACAAAGTTTAccc                                                                                                                                                                                      >  1:339462/1‑250 (MQ=255)
                                                                  cAGGGGAACAGCGATAAAGGCTGCCCTCTTGCCTCACAGATGATTGCCGACGCCGCGCCGTGTCTGCCACTTGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTGGGTGAAAAACTAAAAGGGGCGGTGCTACCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAGAGTTGACCt                                                                                                                                                                                     <  2:426466/251‑1 (MQ=255)
                                                                        aaCACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGcc                                                                                                                                                                               <  2:270104/251‑1 (MQ=255)
                                                                          cacCGGTAAAGGCTGCACTGGACGCACAGAGGAGTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTGCCGTTACCGAAGGTGAAGACAAACCACGGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTAt                                                                                                                                                                           <  2:44006/251‑1 (MQ=255)
                                                                             cGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTTCCTTACCGAAGGGGAAAACAAACCACTGAAATATACGCCTAAGTTTTCCGCCGCCCCGCCGATTCTGCTACACCAAGTTGACCTGTTGCCggttc                                                                                                                                                                          >  1:376775/1‑247 (MQ=255)
                                                                              ggTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTTCCGAAGGGGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCCCGCTGATGCTGCTCCACAAAATTGAACTGTTGGCGTAtct                                                                                                                                                                         >  1:330581/1‑251 (MQ=255)
                                                                               gTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCCGTACCCGAGGTGAAGACAAACCACCGAAATATCCGCATCTGTTTGCCGCCGCCCCCCTTATGCTGCTCAACAAAGTTTACCTGTTGCCggttctc                                                                                                                                                                        >  1:462089/1‑245 (MQ=255)
                                                                                 aaaGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCa                                                                                                                                                                      >  1:51037/1‑250 (MQ=255)
                                                                                        gCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTTCGTTACCCAAGGGGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCCCCGATGCTGCTCCACAAAAATGACCTGTTGCCGTATCTCAAATTTGAc                                                                                                                                                             >  1:394654/1‑251 (MQ=255)
                                                                                                    cacaGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATAGGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTGGCCGTATCTCAACTTTGACGTTGAAAAGTGc                                                                                                                                                 <  1:343069/251‑1 (MQ=255)
                                                                                                      caGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCAt                                                                                                                                               <  2:394654/251‑1 (MQ=255)
                                                                                                      caGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCAt                                                                                                                                               <  1:242815/251‑1 (MQ=255)
                                                                                                       aGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATc                                                                                                                                              <  1:241653/251‑1 (MQ=255)
                                                                                                        gatgatTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCCTTCCGTTACCCAAGGGGAAGACAAACCACTGGAAAATCCCCCTATGTTTTTCGCCGCCTCCCTCATGCTCCTCAACAAAATTGACCTGTTGCCGGTTATCCACTTTGACGTTGAAAAGGGCATCg                                                                                                                                             >  1:97506/1‑251 (MQ=255)
                                                                                                           gatTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCCTATGTTTGCCGCCCCCTCGCCGATGCTGCTCAACAAAGTTGACCCTGTTCCGGATCTCAACTTTTACGTTGGAAAGAGCATCGccc                                                                                                                                          >  1:456748/1‑250 (MQ=255)
                                                                                                                  gACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTTCCGAAGGGGAAGACAAAACACTGAAATATCCGCATATGTTTGCCGCCCCCTCGCTTATGCTGCTCAACAAAGTTGACCCGTTTCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGcccc                                                                                                                                   >  1:43611/1‑250 (MQ=255)
                                                                                                                      ccgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGaa                                                                                                                               <  1:95877/251‑1 (MQ=255)
                                                                                                                       cgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGGCTTCCCTTACCCCAAGGGAAGAAAAACCACTGTAAAAACCGCATATGTTTGCCGCCCCCTCGCTGATGCTGCTCAACAAAATTTATCTTTTGCCGGAACTTCAATTTTACGTTGAAAAATGCATCGCCTGCGCCcgcgga                                                                                                                               >  2:475909/1‑248 (MQ=255)
                                                                                                                                  tGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTTCCGCAGGGGAAGACAAAACACTGAAATATCCGCATATGTTTGCCGCCGCCCCGCTGATGCTGCTCAACAAAATTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGAGCATCGCCTGCGCCCGCGAAGTCAATTCAGa                                                                                                                    >  1:431757/1‑250 (MQ=255)
                                                                                                                                   gCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAg                                                                                                                     <  2:79557/248‑1 (MQ=255)
                                                                                                                                   gCCACTGGACGATAACGGTATTCGGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGa                                                                                                                    <  1:382850/249‑1 (MQ=255)
                                                                                                                                                aaCGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATcc                                                                                                      <  2:560109/250‑1 (MQ=255)
                                                                                                                                                aaCGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCt                                                                                                     <  2:135604/251‑1 (MQ=255)
                                                                                                                                                    gTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGaa                                                                                                                                                       <  2:7814/197‑1 (MQ=255)
                                                                                                                                                    gTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGGGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAATTGACCTGTTGCCGTATCTCAACTTTTACGTTGaa                                                                                                                                                       >  1:7814/1‑197 (MQ=255)
                                                                                                                                                      aTTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCGCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTAttt                                                                                                <  1:23470/250‑1 (MQ=255)
                                                                                                                                                       ttCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGGGAAGACAAAACACTGAAATATCCGCATATGTTTGCCGCCGCCCCGCTGATGCTGCTCAACAAAATTGACCTGTTGCCGTCTCTCAACTTTGACGTTTAAAAGTGCATCGCCTGCGCCCGCGAAGTCAAACCAGAAATTTAAATCATCCTTATTTcc                                                                                              >  1:430170/1‑251 (MQ=255)
                                                                                                                                                             ttATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGccac                                                                                         <  2:339462/250‑1 (MQ=255)
                                                                                                                                                                 cGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTTCTTTCCTTTACCCAAGGGGAAGACAAACCACTGAAAAATCCGCATATGTTTTCCGCCCCCTCGCTGATGCTGCTCAACAAAGGTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTTCATCCCCTGCGCCCGCCAAGTCCATTCCGGAATTGAAATCATCCTTATTTTCGCCCCCAgcg                                                                                    >  1:29814/1‑251 (MQ=255)
                                                                                                                                                                        gTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTc                                                                                                                                                                                                                 <  2:62671/122‑1 (MQ=255)
                                                                                                                                                                        gTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGAGCTTTCCGTTACCGAAGGTGAAGAAAAACCACTGAAATATCCGCATATGTTTTCCGCCGCCTc                                                                                                                                                                                                                 >  1:62671/1‑122 (MQ=255)
                                                                                                                                                                          tGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAggg                                                                            <  1:112257/250‑1 (MQ=255)
                                                                                                                                                                           ggCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCCTTTCCCTTCCGACGGGGGAGGAAAACCAAATCAATAATCGCATCTGTTTGCCGCCGCCTCGCTGGTGCTGCTCAACCAAGATGAACCGCTGCCCTATCCCAAACTTTACGGTTAAAAGGGCATCACCTGCGGCCGCGAAGGCAATCCAGAAATTGAAATCATCCTTTTTTCCGCCCCCAccgggccaggggt                                                                          >  2:409314/1‑242 (MQ=255)
                                                                                                                                                                                  tCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGGGAAAACAAAACACTGAAATATTCGCCTATGTTTTGCGCCGCCTCCCCGATGCTGCTCAACAAAATTGACCTGTTGCCGTATCTCAACTTTTACGGTGAAAAGTGCCTCGCCTGCGCCCGCGAAGACAATCCACAAATTGAAATCATCCTTATTTTCGCCCCCCACGGCGGAGGGGTGGaaca                                                                    >  2:126052/1‑247 (MQ=255)
                                                                                                                                                                                   cGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGt                                                                  <  1:5877/251‑1 (MQ=255)
                                                                                                                                                                                     tATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGGGGAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAATGCATCGACTGCGCCCGCGAAGCCAATCCAGAAAATGAAATCATCCTTATTTCCGCCACCAACGGCGAAGGGATGGACCAGTgg                                                                >  1:160388/1‑251 (MQ=255)
                                                                                                                                                                                         cccGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAAAATCCGCATATGTTTGCCGCCGCCCCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGGAAAGGGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCAACCTTATTTCCGCCACCAGCGGCGAAGGGAAGGaacaagggctgg                                                            >  1:24563/1‑240 (MQ=255)
                                                                                                                                                                                            ggCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAAct                                                         <  2:456748/251‑1 (MQ=255)
                                                                                                                                                                                                     cGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCCCAGGGGAAGAAAAACCACTGAAAAAACCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACCAAGTTGACCTGTTGCCGTATCTCAACTTTTGCGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATTCTTATTTCCGCCCCCCGCGGCGAAGGGATGGACCAGTGGCTTAActggctgggg                                                 >  2:368484/1‑248 (MQ=255)
                                                                                                                                                                                                           cGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGacac                                             <  1:107135/248‑1 (MQ=255)
                                                                                                                                                                                                           cGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAg                                           <  1:2035/250‑1 (MQ=255)
                                                                                                                                                                                                                aaaaaCACAAAGTGGCGGTGCTTTCCGTTACCGAAGGGGAAAAAAAACCACCGGAATATCCGCCTATTTTTGCCGCCGCCCCGCCGATGCTGCTCCACAAAGTTGACCCGTTGCCGTATCTCAACTTTGACGTTGGAAAATGCATCGCCCGCGCCCCCGCAGTCAATCCAGAAAATGAAAACATCCTTATTTCCGCCCCCCCCGGCGAAGGGATGGACCAAAGGCTGAActggctggggaacccacgaggg                                     >  2:234651/1‑237 (MQ=255)
                                                                                                                                                                                                                     cacaAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATAGCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCGCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATa                                <  1:140422/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaaGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGc                             <  2:148439/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaaGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGc                             <  2:339578/251‑1 (MQ=255)
                                                                                                                                                                                                                          aggtgCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGTACCAGTGGCTGAACTGGCTGGAGACGCAGCGATTTGCATAGGCGt                           <  2:462089/247‑1 (MQ=255)
                                                                                                                                                                                                                             ggCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTccc                       <  2:27737/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTccc                       <  1:344103/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGc                    >  1:225222/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTAGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGGGCATAGGAGTTCCCGGc                    >  2:346782/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGcc                   >  2:558098/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ccGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATcc             >  1:106380/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     ggTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGAc  >  2:528117/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
ATTGAAGGCGACCAGCAAACCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGAC  >  NC_000913/2851409‑2851902

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