Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,136,4670AC31.5% 80.9 / 18.5 54E120D (GAA→GACyijEcystine exporter
Reads supporting (aligned to +/- strand):  ref base A (14/23);  new base C (17/0);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.06e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGAGTTATAGCTGGATTTTCATGAAGCAAGTCACCAGTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTGGGC  >  NC_000913/4136170‑4136763
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                        
cTGGAGTTATAGCTGGATTTTCATGAAGCAAGTCACCAGTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACg                                                                                                                                                                                                                                                                                                        <  1:840028/300‑1 (MQ=255)
       tataGCTGGATTTTCATTAAGCAAGTCACCATTTACATCGTTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGa                                                                                                                                                                                                                                                                                                 <  1:368377/300‑1 (MQ=255)
       tataGCTGGATTTTCATGAAGCAAGTCACCAGTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGa                                                                                                                                                                                                                                                                                                 <  1:369030/300‑1 (MQ=255)
                             aaGTCACCAGTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTc                                                                                                                                                                                                                                                                           <  1:831692/300‑1 (MQ=255)
                                       ttACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGc                                                                                                                                                                                                                                                                 <  1:708114/300‑1 (MQ=255)
                                                           ttACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGgtgt                                                                                                                                                                                                                                             >  2:118614/1‑300 (MQ=255)
                                                                gCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAg                                                                                                                                                                                                                                        >  2:404312/1‑300 (MQ=255)
                                                                gCCTTACGCTGCATTTTCGGCGCTCTCGTTTNATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGNTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAg                                                                                                                                                                                                                                        >  2:614918/1‑300 (MQ=255)
                                                                  cTTACTCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCTGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGcc                                                                                                                                                                                                                                      <  1:99081/300‑1 (MQ=255)
                                                                        cTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCa                                                                                                                                                                                                                                <  1:362884/300‑1 (MQ=255)
                                                                                                aTTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAg                                                                                                                                                                                                        >  2:301829/1‑300 (MQ=255)
                                                                                                        tCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTGTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGCGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGc                                                                                                                                                                                                <  1:349653/300‑1 (MQ=255)
                                                                                                            tttATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCAATGTTGGc                                                                                                                                                                                            >  2:248120/1‑300 (MQ=255)
                                                                                                                      ggTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTcc                                                                                                                                                                                  <  1:374225/300‑1 (MQ=255)
                                                                                                                       gTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCg                                                                                                                                                                                 >  1:493387/1‑300 (MQ=255)
                                                                                                                                   gcCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGCAGTGGGCAATGTTTAGCGATTCTGATTGCCGCTTTCGGATTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCAATGCTGGCAATCCTCTACGGCGTCAGTTggg                                                                                                                                                                     >  2:153682/1‑300 (MQ=255)
                                                                                                                                      cGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTgggggg                                                                                                                                                                  <  1:652967/300‑1 (MQ=255)
                                                                                                                                                      cacCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGACGTGGGCAATATTTAGCGATTCTGATTGCAGCTTTAGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGGTGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGc                                                                                                                                                  >  2:293479/1‑300 (MQ=255)
                                                                                                                                                        ccTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGCAGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCAATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTa                                                                                                                                                >  2:365799/1‑300 (MQ=255)
                                                                                                                                                           tAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGTTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCTGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAAc                                                                                                                                             <  1:280842/300‑1 (MQ=255)
                                                                                                                                                                   gCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTAt                                                                                                                                     <  1:339969/300‑1 (MQ=255)
                                                                                                                                                                      cTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGcc                                                                                                                                  <  1:61182/300‑1 (MQ=255)
                                                                                                                                                                             aaaCCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATc                                                                                                                           <  1:509427/300‑1 (MQ=255)
                                                                                                                                                                                               gtCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGCCGTGGGCAATATTTAGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGt                                                                                                         >  2:435534/1‑300 (MQ=255)
                                                                                                                                                                                               gtCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCAGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGt                                                                                                         >  1:338783/1‑300 (MQ=255)
                                                                                                                                                                                                      gcAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAAc                                                                                                  <  1:205957/300‑1 (MQ=255)
                                                                                                                                                                                                                  ggTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGCCGTGGGCAATATTTAGAGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGAt                                                                                      >  2:263911/1‑300 (MQ=255)
                                                                                                                                                                                                                    tCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGc                                                                                    <  1:249252/300‑1 (MQ=255)
                                                                                                                                                                                                                       gCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGt                                                                                 >  1:463187/1‑300 (MQ=255)
                                                                                                                                                                                                                          gAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGGCCGCCTTCTCCTTAGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGTTTATTTTTGGTGTTGCAGCAGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCCCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTCATTGTCGTTAACATCCTGGCAGATGCTGTAta                                                                              >  2:258758/1‑299 (MQ=255)
                                                                                                                                                                                                                          gAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGAAGTGGGCAATAGTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATg                                                                              >  2:404470/1‑300 (MQ=255)
                                                                                                                                                                                                                            ggTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATgcg                                                                            >  1:466971/1‑300 (MQ=255)
                                                                                                                                                                                                                                          ggCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAg                                                              <  1:293189/300‑1 (MQ=255)
                                                                                                                                                                                                                                              aTCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGANTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGtgtg                                                          >  2:625981/1‑300 (MQ=255)
                                                                                                                                                                                                                                               tCCTGAACTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCCCGGTGAACGCCTTCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGTGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGTTGATGAGTGTgg                                                         <  1:238395/300‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCg                                                      <  1:297462/300‑1 (MQ=255)
                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGGCCGCCTTCTCCTTGGGCAATGTTTCGCGATTCTGGTTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGTGCGATTGTTGCTAAACGTCTGTATGCACGTCATCCGCGTGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGc                                                     >  2:265934/1‑300 (MQ=255)
                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGGCCGCCTTCTCCTTAAGCAATATTTAGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCAGTGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGc                                                     >  2:436998/1‑300 (MQ=255)
                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTTCGACGTAGGCAATATTTCGAGATTCTGATTGCAGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGc                                                     >  2:769724/1‑300 (MQ=255)
                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGACGTGGGAAATATTTCGAGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGc                                                     >  2:511976/1‑300 (MQ=255)
                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGNCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGc                                                     >  1:611501/1‑300 (MQ=255)
                                                                                                                                                                                                                                                   gAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGc                                                     >  1:355036/1‑300 (MQ=255)
                                                                                                                                                                                                                                                    aGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCt                                                    >  1:414337/1‑300 (MQ=255)
                                                                                                                                                                                                                                                     gCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCtt                                                   >  1:335532/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        ataCCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTTCTCCTTAAGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTACGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGTGGCGCTGGTGATGAGTGTGGTCGCTTTAc                                                >  2:345058/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        ataCCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGACCGCCTGCGCCGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTAc                                                >  2:817697/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        ataCCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTAc                                                >  1:358875/1‑300 (MQ=255)
                                                                                                                                                                                                                                                          aCCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTg                                              >  1:797085/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                  gTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGGCCGCCTTCTCCTTGAGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCa                                      >  2:440220/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           ttgttTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCAcc             >  1:544307/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                             gttTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCgt           <  1:408618/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 ctcGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTc       <  1:348295/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      tGAACGCCTGCGACGTGTGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGTTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCCCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTGGgc  <  1:132359/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      tGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTGGgc  <  1:249454/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                        
CTGGAGTTATAGCTGGATTTTCATGAAGCAAGTCACCAGTTACATCGGTGCCTTCGACTTTACCGCCTTACGCTGCATTTTCGGCGCTCTCGTTTTATTCATCGTCCTTTTATTACGTGGTCGCGGAATGCGCCCGACACCGTTTAAATACACCTTAGCCATTGCCCTGTTACAAACCTGCGGGATGGTTGGTCTGGCGCAGTGGGCGTTGGTCAGCGGAGGTGCGGGGAAGGTGGCGATCCTGAGCTATACCATGCCGTTCTGGGTGGTGATTTTCGCCGCGTTGTTTCTCGGTGAACGCCTGCGACGTGGGCAATATTTCGCGATTCTGATTGCCGCTTTCGGCTTATTTTTGGTGTTGCAGCCGTGGCAACTCGATTTCTCTTCGATGAAAAGTGCCATGCTGGCAATCCTCTCCGGCGTCAGTTGGGGGGCGAGCGCGATTGTTGCTAAACGTCTGTATGCCCGTCATCCGCGCGTGGATTTATTGTCGTTAACATCCTGGCAGATGCTGTATGCGGCGCTGGTGATGAGTGTGGTCGCTTTACTGGTGCCGCAACGTGAAATTGACTGGCAGCCCACCGTGTTCTGGGC  >  NC_000913/4136170‑4136763

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: