Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,8310TA28.8% 76.9 / 44.7 59A41A (GCT→GCAychSputative uncharacterized protein YchS
Reads supporting (aligned to +/- strand):  ref base T (15/22);  new base A (10/5);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAA  >  NC_000913/1286541‑1287129
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                           
acaTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                                                                                                                                                                                                                    >  1:128865/1‑300 (MQ=255)
acaTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTAGCTTATCCGAACTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                                                                                                                                                                                                                    >  1:304985/1‑300 (MQ=255)
      aCCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCTGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTg                                                                                                                                                                                                                                                                                              <  1:98385/300‑1 (MQ=255)
        ccTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                                                                                                                                                                                                                                            <  1:439713/300‑1 (MQ=255)
                      gTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTg                                                                                                                                                                                                                                                                              >  1:223098/1‑300 (MQ=255)
                                  gCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGttgggaagttcaggga                                                                                                                                                                                                                                                                  <  1:628130/300‑15 (MQ=255)
                                   acaggTGTGGGGTGAGTTGGCGTATGCGCTTCGCTGATCCGGCCGACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGttgggaagttcagggac                                                                                                                                                                                                                                                                 <  2:593129/297‑16 (MQ=18)
                                       gTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTNCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAGTTCCGGGGCtttt                                                                                                                                                                                                                                                            >  2:580593/1‑300 (MQ=255)
                                       gTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGttgggaagttccgggacctgt                                                                                                                                                                                                                                                             >  1:603017/1‑281 (MQ=255)
                                           ggggTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTAATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGaaa                                                                                                                                                                                                                                                        >  2:355035/1‑300 (MQ=255)
                                           ggggTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTTGAAAGTTCAGGGACTTTTGaaa                                                                                                                                                                                                                                                        >  2:702485/1‑300 (MQ=255)
                                                gAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGag                                                                                                                                                                                                                                                   >  1:447152/1‑299 (MQ=255)
                                                  gTTGTCGGATGCGCTTCGCTTATCCGACCTAAAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTAGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGGAGTTCAGGGGCTTTTGAAAGTGAtgg                                                                                                                                                                                                                                                 >  2:410429/1‑300 (MQ=18)
                                                             cgcttcgcttATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGCAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATAGTGGTGGGGGaa                                                                                                                                                                                                                                      >  1:298924/1‑300 (MQ=255)
                                                             cgcttcgcttATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                                                                                                                                                                                                                      >  1:685740/1‑300 (MQ=18)
                                                             cgcttcgcttATCCGACCTACAGGGGAGGATATTGTAGGCCCGGGAAACGTGGTGCCGTCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGTAAGCGATGTTGGTGGTGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCTTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTTATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGaa                                                                                                                                                                                                                                      <  2:606209/300‑1 (MQ=25)
                                                                      tattcGACCTACAGGGGAGGATATTGTATGCCCGGTAAACGCGGTGCCATCGGGCAAAATTATCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                                                                                                                                                                             <  2:262530/296‑1 (MQ=255)
                                                                      tATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAcgc                                                                                                                                                                                                                             <  1:750987/300‑3 (MQ=255)
                                                                      tATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                                                                                                                                                                             <  1:463038/300‑1 (MQ=17)
                                                                      tATCCGACCTACAGGGGAAGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTTAGAATTCCGGGGCTTTTGAAAGTGATAGTGGTGGTGGAAGGATTACTc                                                                                                                                                                                                                             <  2:269055/300‑1 (MQ=255)
                                                                        tCCGACCTACAGGTGAGGAGATTGTATTCCCGGTAAACGTGGGGCCATCGGGCAAAATTAGCTCGTAGTAACAGGTTTTGATTACTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGTGGAAGGAttcgaacc                                                                                                                                                                                                                           <  2:128433/300‑7 (MQ=14)
                                                                           tacctacatttgaggaTATTGTAGGCCCGGTAAACGTAGTGCCATCGGGCAAAATTAGCTCGGATTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCNCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGTCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAANTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                        <  2:562295/289‑1 (MQ=255)
                                                                           gaccgacgggggAGGATATTGTAGGCCCGGGAAACGTGGTACCATCGGGCAAAATTAGTTCGGGGTAACGGGTTTTGATGATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                        <  2:347462/292‑1 (MQ=21)
                                                                           gACCTACAGTGGAGGATATTGTAGTCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGTGGCTTTTGAAATTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                        <  2:403600/300‑1 (MQ=255)
                                                                           gACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCNCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGGATAGTTGATTGCGGATTCGTTTGAGAATTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTccgctc                                                                                                                                                                                                                        >  2:183081/1‑295 (MQ=255)
                                                                           gACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                                                        <  1:471585/300‑1 (MQ=255)
                                                                            aCCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcgaaccttcg                                                                                                                                                                                                                       >  1:679993/1‑290 (MQ=14)
                                                                            aCCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGNGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcganccttcg                                                                                                                                                                                                                       >  2:638895/1‑290 (MQ=17)
                                                                                acatgttaggaTATTGTATGCCCGGTAAACGTGGTGCCAGCGTGCAAAACCAGCCCGGACTAACAGGTTTTGACTCTTCACGGAATCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGGCGTTGCCTGCGGCAACGCCCTCTCGCTGGCGCTCAAGTCGAACTTTGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATGCGATTGCTGATTCGTGGGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                                                                                                                                                                                                                   <  1:20139/293‑1 (MQ=17)
                                                                                aCAGTGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                                                                                                                                                                                                                   <  1:416235/300‑1 (MQ=255)
                                                                                aCAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                                                                                                                                                                                                                   <  1:500456/300‑1 (MQ=255)
                                                                                   ggggAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACTCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCCGCGCTACGCGCt                                                                                                                                                                                                                >  1:679377/1‑300 (MQ=255)
                                                                                        ggaTATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCTCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTTGGGAATTCAGGGACTTTTGAAAGTGATGGTGGTgggggtgtacctcatccacatctcccactccacc                                                                                                                                                                                                           >  1:706520/1‑271 (MQ=18)
                                                                                             ttGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcaccccctctacgtccctcactcccca                                                                                                                                                                                                      >  2:429467/1‑273 (MQ=14)
                                                                                                     ccGGTAAACGCGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGTTTTTGATTATTTAAGGAGGCGGTGTTGGTGTGGGAAGAATTACTCAGCGCTGCGCGCTGTGCCCGTCGGGTCGTTGCCTGCGTCAANGCTCGCGAGCTGGCGCTCAAGTCGAACCTTGGTCGAAGGTTCTCATCCTTCCCCGCGTGGGCAGAGTATTTGATTGCGGAGTCGTTTGAGAATTCAGGGGCTTTTGAAAGTGATGTTGGTGGGGTAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtgg                                                                                                                                                                                              <  2:183830/300‑3 (MQ=18)
                                                                                                      cGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGc                                                                                                                                                                                             >  1:20957/1‑300 (MQ=255)
                                                                                                                                            tAACAGGTTTTGATTATTTAAGGAATCGATGGTGGTGGGGGAAGGACTGCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTGGATTGCGGATTCGTTTGAGAATTCCGGGGCTTGTGAAATTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAcc                                                                                                                                                       <  2:351988/300‑1 (MQ=255)
                                                                                                                                                                aGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATcc                                                                                                                                   <  1:275649/300‑1 (MQ=255)
                                                                                                                                                                          gttggtgtGGTAAGGATGATTCAGCACTGCGCGCTTCGCCCTTTGGTTCGTTGCCTGCGGAAACGCTCTCTCGCTGGCGCTCAATTCGAACCTGGGGCGAAGCTTCTCATCCTTCCCCGCTTGTGCAGAATATGTGAGTGCTGATTCGTTGGAGAATTCCGGGGCTTTTGAAATTGATGGTGGTGGTGGAAGGATTACTCAGCGCTGCACACTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCCGACGCTCGAGCCGAACCTTATTCGAAGCTTCTCATCCTTCCCCGCac                                                                                                                         <  2:262038/298‑2 (MQ=17)
                                                                                                                                                                                   gggaGGATTACTCATCGCTGCGTGCTTCGCTCTGCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTGGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                                                <  2:740561/297‑1 (MQ=255)
                                                                                                                                                                                   ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCTGGGCTTTTGAAAGTGATGGTGGTGGTGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaac                                                                                                                <  1:697562/300‑2 (MQ=255)
                                                                                                                                                                                       ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                                            >  1:451427/1‑300 (MQ=255)
                                                                                                                                                                                       ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCCGCGCTGCGCGCTTAGCCCTTAGGGTCGTTGCCTGCGGCAACGTTCTCTCGCTGACGCTCGAGCAGAACCTTAGTCGAGGCTTCTAATCCTTcacagtatgggccgaacagtt                                                                                                            >  1:185039/1‑280 (MQ=14)
                                                                                                                                                                                            aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTAGAACCTTGGTCGAAGCTTCTAATCCTTCACCGAATGGGCAGAATATTTCATTg                                                                                                       >  2:588613/1‑300 (MQ=14)
                                                                                                                                                                                                 gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAACTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGGGAGTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGTTTAGCCCTTAGGGTCGTTGCCTGTGGCAACGCTCTCTCGCTGACGCTCGAACCGAACCTTAGTCGCAGCTTCTCATCCTTCCAAGAATGGGAAGAATAtttagctgcagat                                                                                                  >  2:385925/1‑290 (MQ=14)
                                                                                                                                                                                                   gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCCCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                                                >  1:138475/1‑300 (MQ=255)
                                                                                                                                                                                                          gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                                         >  1:728979/1‑300 (MQ=255)
                                                                                                                                                                                                                     cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCAAGTCTAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGCCAGAATATTTGATTGCGGATTCGTTGGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCNTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATNTGATTGCGGATTCGCTTGAGAGTTCAGGGAc                                                                              <  2:552603/300‑1 (MQ=255)
                                                                                                                                                                                                                              gacTGCGGCAACGCGCTCTCGCTGGGGATCAATTCGAGCTTTGGTCGAAGCTTCGCATCCTTACCCGCTGGGTCAGAATATTGGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAATTGATGGTGTTNGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                                     <  2:187418/298‑1 (MQ=17)
                                                                                                                                                                                                                                             ctctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                      >  1:158258/1‑300 (MQ=255)
                                                                                                                                                                                                                                                       cgcTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                                            >  1:691469/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        gcTCAATTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGTGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                           <  1:281327/300‑1 (MQ=255)
                                                                                                                                                                                                                                                        gcTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAATTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                           <  1:609958/300‑1 (MQ=255)
                                                                                                                                                                                                                                                        gcTCAAGTCAAACCCTGGTCTAAGCTTCTCATCCTTCCCCGCTTGTGCAGAATATTTGATTGCTGATTCGGGTGAGAATGCCGGGGCTGTTGAAAGTGATGGTGTTGGGGGAAGGACGACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                           <  2:137964/300‑1 (MQ=17)
                                                                                                                                                                                                                                                             ggtcggcatttgttcggattttctacgccTTCCCCGCGTGGGTAGAATATCTGACGGCTGATTCGTGGTAGAGTTCCGGTGTTTGCGAAATTGATTGTAGTGTGTNAAGGACTACTCAGCGCTGCGCGCTTCTCCCTTCGTGTCGTTGCCTGCGGCAGCGCTCTCGCGTTGACGCTCGGGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAGTATGGGATTGCGGATTCGCGTTAGAGTTCGGGTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                      <  2:687408/273‑1 (MQ=9)
                                                                                                                                                                                                                                                                           cGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGccc                        <  1:467852/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  tctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCATGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTCGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGt                 >  1:209225/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                          ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGATTCTCACCCTTCCCCGCCTGGGCAGAATATTTAATTGCGGACTCGCTCGAGATCCCAGGGACTTTCGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGTGCTTCGCCCTTAGGGTCGTTGCCt         >  2:765350/1‑300 (MQ=17)
                                                                                                                                                                                                                                                                                                 cTTGGGCAGAATATTTGATTGCGGATTCGTGGGAGAATTCCGTGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGTTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGCTGGAGAGTTCAGGGACTTGTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa  <  1:472571/300‑1 (MQ=17)
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                           
ACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAA  >  NC_000913/1286541‑1287129

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