Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7190TA28.9% 158.7 / 28.1 76intergenic (+52/‑45)ytcB/nrfAprotein YtcB/cytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base T (48/6);  new base A (0/22);  total (48/28)
Fisher's exact test for biased strand distribution p-value = 5.18e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTG  >  NC_000913/4287472‑4287962
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
aaGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                                                                                                                                                                                                                                    >  1:431575/1‑250 (MQ=255)
  gTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaaaaa                                                                                                                                                                                                                                                 <  1:837289/251‑1 (MQ=255)
    aGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAc                                                                                                                                                                                                                                                <  2:817929/250‑1 (MQ=255)
             cAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGc                                                                                                                                                                                                                                      >  2:150174/1‑250 (MQ=255)
              aaaTAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGca                                                                                                                                                                                                                                     >  1:167482/1‑250 (MQ=255)
                aTAAACGGCAGGGTAATTTTTGAAGGACAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGGGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaac                                                                                                                                                                                                                                   >  1:933166/1‑250 (MQ=255)
                           ggTAATTGTTGAAGTCCAGCAAAAAAAGTGTATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTGGAAAAAGAGAAAGATACTGACTAACTCTAAAGTGGTATTTTACAGGCACTTACATTGTGTTAAATACAACATTCACAGTGTGGTTTTTTGTTCAACTTGGGGGCGAGAAATGTAATGACAGTTGAGGGGCGGTTCCTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAtg                                                                                                                                                                                                                       <  1:788685/251‑2 (MQ=255)
                            gTAATTTTTGAAGGTAAGCAACAAAATTTGATTAATTCTTTGAGGAACATGCAGTTATGTATGTGTTGAAAAAAGAGAAAGATACTGACTACCTCAAAAGTGTTTTTTTACATCCATTTAAAATTGATAAAAGACAAAATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTGGGTGGGTTTACTGTCGGGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                                                                                                                                                      <  2:1184274/251‑1 (MQ=255)
                              aaTTTTTGAAGGTCAGAAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTAGGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCATTTACAATTGATTAAAGACAACATTCACAGTGTGTTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                                                                                                                                     <  2:69662/250‑1 (MQ=255)
                                      aagttCAGAAACAAAAGTTAATTAATTTTTGGGGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATCTGGACTAACTCAAAATTGGTATTTTCCTGGCACTTAAAATTAATAAAAGAAACACTTCACAGTGTGGTTATTTTTGACACATAGGGCCGAGAAATGTCATGAGATTGTAGGGGCGGTTTCTATCTTGAAAACCAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGcc                                                                                                                                                                                                            <  1:865410/247‑1 (MQ=255)
                                           aagaaaaaaaaGTTGATTAATTTTTTGTGGAACATGAATTATTGCATGTGGTTGAAAAATAGAAGGATACTTCCAACCTCTAAAGTGGTTTTTTACAGGCCCTTACAATTGATAAAAACAACAATCCAGAGTGTGTTTTTTTGTTACACATAGGGGCGAGCATTGTCATGACGGTGTGGTCGCGGTTACTGTCGTAAAAAACAAAAAAAGCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                                                                                                                       <  2:568320/244‑1 (MQ=255)
                                               aacaaAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATg                                                                                                                                                                                                                 >  2:511105/1‑237 (MQ=255)
                                               aaaaaaaGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTTTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCCCTTACAATTGATTAAAGACAACATTCACAGTGTGTTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTGGGGGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATg                                                                                                                                                                                                                 <  1:511105/234‑1 (MQ=255)
                                                  aaaaGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAgg                                                                                                                                                                                                >  1:555832/1‑251 (MQ=255)
                                                    aatttGATTATTTCTGTGAGGAACATGATGTTAGGAATGTGGTAAAAAAAGAGGAAAATACTGACAAACTCTAAAGTGGTTTTTTAATTGCACTTAAAATGGATTAAAGAAAACATTCACGGGGTGGTTATTTGTACCACATGGGGGCGAGCAATGTCATGACAGTGTAGGTGGGTTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGAt                                                                                                                                                                                              <  2:787110/248‑1 (MQ=255)
                                                        tgcttAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGACAACGAGGATGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaac                                                                                                                                                                                           >  2:354147/4‑249 (MQ=255)
                                                         gATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaaa                                                                                                                                                                                          >  1:170412/1‑250 (MQ=255)
                                                              aTTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac                                                                                                                                                                                    >  1:323416/1‑251 (MQ=255)
                                                              aTTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac                                                                                                                                                                                    >  2:114701/1‑251 (MQ=255)
                                                                                    atgcatgcTGTTAAAAAAGAGAAAAATACTGATTAACTCAAAAGTGGTTTTTTACATGCATTTAAAATTAATAAAAGACAACATTCACAGTGGGTTTTTTTGTTACACATGGGGGCGAGAATGGTCATGACGGTGTAGGGGGGGTTACTTTCTTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTa                                                                                                                                                               <  1:1140069/250‑1 (MQ=255)
                                                                                       gcatgcTGTTGAAAAAGTGGAAGATACTGACCACCTCTAACGCGGTATTTTACATACACTTACAATTGATTAACAACACCATTCAAAGTGTGGCTATTTGTTACACATCGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTAGTGAAAAACAATAAAAAACGCCATTGCAACAATGGCGCAATTCGGCTGAAGCCACTCTGACAAGGATAACAATAAACGCACGCCGTATCTCCAGCTTATTg                                                                                                                                                            >  2:308513/1‑251 (MQ=255)
                                                                                        catgcTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                                                                                                                            >  1:589382/1‑250 (MQ=255)
                                                                                                  aaaaaaagaggagaatacaactaaaCTCAAAAGTGGTTTTTTACATTCGTTTAAAATTTATTAAAAACACCTTCAACAGTGGGGTTTTTTTTTAAAAATGGGGGGGAGGAATGTCATGACGGGGTGGGTGCGTTTATTGTCAGAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt                                                                                                                                                   <  2:871541/228‑1 (MQ=255)
                                                                                                    aaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                                                                                                                               >  2:58446/1‑251 (MQ=255)
                                                                                                     aaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctt                                                                                                                                              >  2:142766/1‑251 (MQ=255)
                                                                                                            aaGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc                                                                                                                                                       >  2:631933/1‑236 (MQ=255)
                                                                                                             aaaatcctgacaacctcaaaaGGGGTATTTTACATGCATTAAAAATTGATAAAAGAAAACATTCACAGGGTGTTTTTTTTTTACACATGGGGGCGAGAATTGTCATGACAGTTTAGGTGGGTTTCCTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc                                                                                                                                                       <  1:631933/218‑1 (MQ=255)
                                                                                                                tACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGCCAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                                                                                                                                    >  2:881086/1‑251 (MQ=255)
                                                                                                                   tGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGAGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTATTTCTTTTTCACTTCTGtt                                                                                                                                 >  2:1154970/1‑251 (MQ=255)
                                                                                                                    gACTAACTCAAAAGTGGATTTTTACAGTCATTTACAATTGATAAAAGAAAACTTTACAAGTTTGGTTTTTTTTTCCACTGGGGGGGGACAAATGTCATGACAGTGAGGGGGGGTTTCCTGTGGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt                                                                                                                                 <  1:789518/250‑1 (MQ=255)
                                                                                                                         aaCTCTAAAGTGTTATTTTACATGCATTTACAATTTATTAAAGACAACATCCAAAGTGTGTTTATTTGTTACACATAGGGGCGAGCATTGTAATGAAAGTGTGGTTGCGTTTCCTGCCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACg                                                                                                                             <  1:931567/250‑1 (MQ=255)
                                                                                                                          aCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCt                                                                                                                           >  2:756483/1‑251 (MQ=255)
                                                                                                                             ataagttGGTATTTTACATGCCTTAAAAATTGATAAAAGAAAAAATTCAAAGGGGGGTTTTTTGTTAAACAGAGGGGAAAAAAAGGTAAGGACGGTGTGGGTGGGGTTACTGTTGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGaa                                                                                                                        <  2:434213/245‑1 (MQ=255)
                                                                                                                               aaaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGAAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAAc                                                                                                                       >  2:692908/1‑250 (MQ=255)
                                                                                                                                aaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAAAATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCACTTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACa                                                                                                                      >  2:981770/1‑250 (MQ=255)
                                                                                                                                 aGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa                                                                                                                     >  2:390263/1‑250 (MQ=255)
                                                                                                                                 aGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa                                                                                                                     >  1:293194/1‑250 (MQ=255)
                                                                                                                                         tttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTAt                                                                                                                                                                         >  2:337711/1‑190 (MQ=255)
                                                                                                                                         tttACATGCACTTACAATTGATTAAAGACAACATTCACAGGGTGGTTATTTGTTACACATGGGGGCGAGCATTGTAATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTAt                                                                                                                                                                         <  1:337711/190‑1 (MQ=255)
                                                                                                                                           tACATGCACTTACAATTGATTAAAGAAAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAAAAATAAAAACCGCCATTGCAACAATGGCGAAATTCGGATGAAGCCCCTATGAAAAGGATAACAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttta                                                                                                                                          >  2:635082/1‑218 (MQ=255)
                                                                                                                                             cATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCg                                                                                                        >  2:762024/1‑251 (MQ=255)
                                                                                                                                               tGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGc                                                                                                       >  1:993158/1‑250 (MQ=255)
                                                                                                                                                       cAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGt                                                                                              >  2:654647/1‑251 (MQ=255)
                                                                                                                                                            tGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAt                                                                                                                                                                                                                        >  2:609668/1‑125 (MQ=255)
                                                                                                                                                            tGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTCCTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAAt                                                                                                                                                                                                                        <  1:609668/125‑1 (MQ=255)
                                                                                                                                                                ttAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAAc                                                                                                                       >  2:981058/1‑219 (MQ=255)
                                                                                                                                                                 tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTgg                                                                                      >  2:154476/1‑251 (MQ=255)
                                                                                                                                                                 tAAAGACAACATTCACAGTGTGGGTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGAATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACtgggg                                                                                      >  2:374166/1‑248 (MQ=255)
                                                                                                                                                                      acaacaTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCg                                                                                 >  1:448348/1‑251 (MQ=255)
                                                                                                                                                                       caacaTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCg                                                                                 >  1:710072/1‑250 (MQ=255)
                                                                                                                                                                       aaacaTTCACAGTGTGTTTATTTGTTACACATGGGGGCGAGAAATGCCATGACAGTGTGGGTGCGTTTACTGTCGTAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCga                                                                                <  1:483858/250‑1 (MQ=255)
                                                                                                                                                                        aacaTTCACAGTGTGGTTATTTGTTACCCATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGACGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTAACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCgc                                                                                >  2:265190/1‑249 (MQ=255)
                                                                                                                                                                          caTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCga                                                                                >  2:209168/1‑248 (MQ=255)
                                                                                                                                                                             ttCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAAt                                                                           >  2:759743/1‑251 (MQ=255)
                                                                                                                                                                              tCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCCCAAAACCTGTAACTGTGGACGCGAAGAATg                                                                          >  2:800575/1‑251 (MQ=255)
                                                                                                                                                                                      tgtgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCt                                                                    >  2:688371/1‑251 (MQ=255)
                                                                                                                                                                                          gTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGc                                                                >  1:573149/1‑251 (MQ=255)
                                                                                                                                                                                                 gTTACACATAGGGGCGAGCAATGTCATGAGAGTGTAGGTGGGTTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCgcag                                                         <  1:654647/251‑1 (MQ=255)
                                                                                                                                                                                                  tcacacATTGGGGCGAGCAATGTCATGACAGTGTAGGTGCGATTACTGTCGTGAACAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACCAGGATAAAAATAAACGCACGCCGAATCTTCAGCTTATTGATGCCTTTTTTCTTTTTCATTTCTGTTCACGCTGGACAAACGGCTGCTCCCGCAAAACCTGTAACGGGGGAAGCGAAGAATGAAACCTTTGCCCCgcagc                                                        >  2:651792/3‑251 (MQ=255)
                                                                                                                                                                                                           ggggCGAGAAATGTCAAGACAGTTGAGGTGCGGTTTTTTTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATc                                               <  1:800575/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCAATGTCATGACAGTGTAGGTGCGTTTACTGTCGAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATc                                               <  1:64918/244‑1 (MQ=255)
                                                                                                                                                                                                                  gCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATc                                               >  2:64918/1‑244 (MQ=255)
                                                                                                                                                                                                                          taatgacggggtgggggggtttACTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCt                                   <  1:762024/230‑1 (MQ=255)
                                                                                                                                                                                                                          gagtccAGTGTAGGGGCTGTTGCTGGCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAACCGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTg                                  <  1:1026485/245‑1 (MQ=255)
                                                                                                                                                                                                                           attaccattttgggtcgttttACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTg                                  <  1:398722/231‑1 (MQ=255)
                                                                                                                                                                                                                                     aGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTc                      >  1:652644/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gggggggtttACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTc                      <  1:154476/243‑1 (MQ=255)
                                                                                                                                                                                                                                       gTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCgg                    >  2:510000/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gcgtttACTGTCGAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGaa                  <  1:58446/247‑1 (MQ=255)
                                                                                                                                                                                                                                              ttACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGt              >  1:1109058/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTc             >  1:407013/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  cgttCGTAAAAAAAAAAAAAACCCGCTATGGCAACATTGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATGAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCag           <  1:374166/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCag           >  2:1034607/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     tggtaaaaaaCAAAAAAAGCCGCCAGTGCAACAATGGCGCAATGCGGATGAAGCCCATATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCagag         <  2:833021/245‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAAAACAATAAAAACCGCCATTGAAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCgtgt    >  2:816272/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           aaaaCAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTg  >  1:820432/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
AAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTGTTG  >  NC_000913/4287472‑4287962

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