Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,248,4240CA26.7% 128.1 / 18.5 60V37F (GTT→TTT) lysPlysine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (39/5);  new base A (0/16);  total (39/21)
Fisher's exact test for biased strand distribution p-value = 1.36e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAATCTCTTTTCTGGCG  >  NC_000913/2248186‑2248672
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          
gCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTg                                                                                                                                                                                                                                                 <  1:554855/250‑1 (MQ=255)
gCAACCAGGCCACCGCCGCTAGTCAGCGCCCAGTTGTAGTAGTATGTCGTGCACAGCGCGAAGCCAATGCCTTCTTCAACATAGTTCTGACCGTAAGCGGCAAACGACCCTGAAACCGGCATATATGCAGCCAGTTCACCGACACGGGTCATCAGGAATTAAACCATCAGGCTAAACTGCATAAACAAGAGCAATGCCCCGCCCGGTCCCGCCTGAGAAATCGTTTCGCCAGAGGccacaataatacccg                                                                                                                                                                                                                                                 <  2:716919/250‑10 (MQ=255)
           aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCCAAGCCTTCTTCAACATAGGTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGCGAAATCGTTGCGCCAAAGGCAACAAAAAGACCTGTACCGATGGAAc                                                                                                                                                                                                                                    >  2:552334/1‑251 (MQ=255)
           aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGAACCGATGGaa                                                                                                                                                                                                                                     >  2:833856/1‑250 (MQ=255)
                 gATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACAGCCa                                                                                                                                                                                                                              >  2:525895/1‑250 (MQ=255)
                                          tagtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACgcgc                                                                                                                                                                                                    >  1:107237/1‑251 (MQ=255)
                                               ccagcccagcGCGAAACCAAACCCTTCTTAAACATAGTCATGACCGTAAGGCGAAAACAAACCGGAAACCGTCATATAGGCACCCAGTGCGCCGACACGGGTAATCAAGAAGAAAACCCACGGGCCAACCAGCAAAAGGGAGAGAAATCCCCCGCCCGGACCTGCCTGAAAAAGCGTGGCCCGAGAGGAAAAAAAAAGCCCGGCCCCGACGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTa                                                                                                                                                                                               <  2:22371/251‑1 (MQ=255)
                                                  gcccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTTAGGTGACGCGCCTTTAAtt                                                                                                                                                                                            >  1:679642/1‑251 (MQ=255)
                                                      agcGCGAAGCCAAAGCCTTCTTCTACATAGTTATGACCGTAAGTGGCAACCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAAGCCAATCAGCATATACGAGAGCAATGCCCCGCCAGAACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAACAGACCTGTACCGATGCAACCGCCAATGGCAATCATCGTCAGGCGACGCGCCTTTAATTCACg                                                                                                                                                                                        >  2:103493/1‑251 (MQ=255)
                                                       gcgcGAACACAAAGCCTTCTAAAACATGGTTCGGACCGAAAGGGGAAAACGAACCGAAACCCGGCATATATGCAGCCAGTCCGCCGAGACTGGTCATCAGAAAGAAAACCATCGGGCAAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCGGCCGAAGAAATCGTTGCGCAAGGGGCAAAAAAAAGCCCTTTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACgg                                                                                                                                                                                       <  1:283513/251‑1 (MQ=255)
                                                              gCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGCCCTGCCTGAGAAAGCGTTGCGCCAGAGGCAACAAAAAGACCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGc                                                                                                                                                                                <  1:553318/251‑1 (MQ=255)
                                                                 aaaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCg                                                                                                                                                                             >  2:583269/1‑251 (MQ=255)
                                                                   aGCCTTTTTCAAAATAGTCAGGACCGAAAGGGCAAACCGAACCGTAACCCGGAAATTATGCACCCAGTTCACCGAGCGTGGTCATCAGGAATAAAACCATGAGCCAAACCACCATATGGTAGAGTAATGCCCCGCCCGGTCTGGCCTGAAAAATGGTTGCGCCAGAGGCAAAAAAAAGACCTGTTCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                                                                                                                           <  1:145680/251‑1 (MQ=255)
                                                                   aGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAatc                                                                                                                                                                                                                     >  1:651895/1‑209 (MQ=255)
                                                                   aGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                                                                                                                           >  2:710857/1‑251 (MQ=255)
                                                                   aGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                                                                                                                           >  1:574240/1‑251 (MQ=255)
                                                                    ccctttttaaaaatatttCTGACCGAAAGGGACAAAGAACCGGGCACCCGGCATATATGCAGCTGTTTCACGGAGCCGGGCCATAAGAAAGTAACCCACCGGGCCAACCACCATATACACGAGAAGTGCCCCGCCCGGCCTTGCCGGAAAATTCGTTGCGCAAGAGGAAAAAAAAAGCCTTGTCCCGAGGGATCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgcg                                                                                                                                                                          <  1:651921/235‑1 (MQ=255)
                                                                             aaCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAgg                                                                                                                                                                                                            >  1:44785/1‑208 (MQ=255)
                                                                                    tttggcccaaaaggggaaaaCGAACCGAAACCCGGCATATATGCCGCCAGTTCCCCGAGCCGGGTCAAAGGGAAGAAACCCACCAGGCAAATCAGCATATACGGGAGCAATCCCCCGCCCGGACCGGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAg                                                                                                                                                           <  1:740272/233‑1 (MQ=255)
                                                                                         acagTAAGAGGCAAACGAACCGGAAACCGGCATATATTCAGCCAGATCACCGAGACTGGTCATCCGGAAGTAAACCATCAGGCCACTCAGCATATACGAGAGCAATGCCCCTCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCTTTTACTTACCGCCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTctg                                                                                                                                                     >  2:673789/4‑249 (MQ=255)
                                                                                         aCCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCgg                                                                                                                                                     >  2:751950/1‑251 (MQ=255)
                                                                                               aGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAcc                                                                                                                                                >  1:821025/1‑250 (MQ=255)
                                                                                                             ccGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAACGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCACGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACACAAATGCtt                                                                                                                                  >  2:71590/1‑250 (MQ=255)
                                                                                                               ggAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc                                                                                                                                >  1:141650/1‑251 (MQ=255)
                                                                                                               ggAAACCGGCATATAGGCAGCCAGTCCACCAAGACGGGCAACAAGAAAGAAAACCACCAGCCCAATCAGCATATCGGAGAGCAAGCCCCCGCCCGGCCCTGCCTGAAAAATCGTGCCGCAAGAGGAAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc                                                                                                                                <  1:383288/251‑1 (MQ=255)
                                                                                                                gcaaCCGGCATATATGCACCCAGTTCCCCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGCGAGCAATGCCCCGCCCGGACCTCCCTGAGCAATCGTTGCGCCAGAGGCGACAAACAGACCTGTACCGCTGGAACCGCCAATGGCAGTCATCGTCAGGTGCCGCGCCTTTAATTACCCGCGTAAGCCCGGCGCTTCTGTGGTTTTAATTTCGGAATCCATACAAACATGATATcc                                                                                                                               >  2:399329/3‑251 (MQ=255)
                                                                                                                       ggCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                                                                                                                          >  1:469476/1‑250 (MQ=255)
                                                                                                                         cATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCa                                                                                                                              >  2:179041/1‑244 (MQ=255)
                                                                                                                          atatatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaa                                                                                                                       >  1:599894/1‑250 (MQ=255)
                                                                                                                                 cAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAACCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgc                                                                                                               >  1:262861/1‑251 (MQ=255)
                                                                                                                                   gccagtccaccgagctgtgtaagaaGGAAGAAAACCAGAAGGCCAATCAACAAATACGAGAGAATTCCCCCGCCGGGCCCTGCGTGAGAAATCGTTGCGCCAGAGACAAAAAAAAGACCTGCCCAGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcgc                                                                                                             <  2:454354/227‑1 (MQ=255)
                                                                                                                                   gCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcgc                                                                                                             >  1:43595/1‑251 (MQ=255)
                                                                                                                                     cAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGcgcgc                                                                                                           >  2:368668/1‑250 (MQ=255)
                                                                                                                                         tCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt                                                                                                       >  1:16905/1‑251 (MQ=255)
                                                                                                                                              gagaCTGGTCAGCAGGAAAAAAACCACCAGACAAATAAGAATATCGAAGAGCAACCCCCCGCCCGGCCCTGCCGGAAAAATCGTGGCGCGAGGAGCAAAAAAAAGCCTGGTCCCGGTGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAg                                                                                                  <  2:469476/251‑1 (MQ=255)
                                                                                                                                               gagactggcaaacaggaaaaaacccATCAGCCCAATGAGCTTATACGAGAGCATTGCCCCGCCCGGACCTCCCGTAAAAATCTTTGCGCAAGAGACAAAAAAAAGCCCTGTACCGGGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAg                                                                                                  <  1:44391/228‑1 (MQ=255)
                                                                                                                                                 aCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACg                                                                                               >  2:246007/1‑251 (MQ=255)
                                                                                                                                                        catcagaaagaaacccacaagccaaaTCAGCATATACGAGAGCAGTCCCCCCCCGGGCCCTGCATGAGAAATCGTTGCGCAAGAGACAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGc                                                                                         <  1:751950/227‑1 (MQ=255)
                                                                                                                                                         tcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCa                                                                                        >  1:456614/1‑250 (MQ=255)
                                                                                                                                                               aaGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAGTGGCACTCCtt                                                                                  >  1:168095/1‑251 (MQ=255)
                                                                                                                                                                 gaaaccaatcgggcaaaTCAGCATAGACGAGAGCAGTCCCCCGCCCGCCCCTCCGGGAGAATTGGTGGCGCGAGAGACAAAAAAAAGCCTTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTcc                                                                                <  1:751728/236‑1 (MQ=255)
                                                                                                                                                                    aaCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGc                                                                             >  2:268067/1‑251 (MQ=255)
                                                                                                                                                                           aGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGaa                                                                                                                                                                                                                                     <  1:316804/93‑1 (MQ=255)
                                                                                                                                                                           aGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGaa                                                                                                                                                                                                                                     >  2:316804/1‑93 (MQ=255)
                                                                                                                                                                                    aGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTtat                                                             >  1:365761/1‑251 (MQ=255)
                                                                                                                                                                                     aaataaaCGAGAGAAATGCCCCCCCCGGCCCTGCCTGAGAAATCGTGGCGCAAGAGGAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTtata                                                            <  1:570768/245‑1 (MQ=255)
                                                                                                                                                                                          atacgggggcagcgCCCCCCCCGGACCGGCCTGAAAAATCGGGGCGCAAGAAGCAAAAAAAACCCCTGCCCCGAGAAACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGac                                                        <  2:389696/238‑1 (MQ=255)
                                                                                                                                                                                            aagagacaaaTGCCCCCCCCGGCCCTGCCTGAAAAATCGTGGCGCCAGAGCCAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGacac                                                      <  1:805723/243‑1 (MQ=255)
                                                                                                                                                                                                    aaTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCACTGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATg                                               >  2:234596/1‑250 (MQ=255)
                                                                                                                                                                                                     aTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGAt                                             >  1:331693/1‑251 (MQ=255)
                                                                                                                                                                                                      gcccccGCCCGCACCTGCCTGAGAAACGGTGGCGCCGGAGGCAAAAAAAAGCCTTGTCCCGATGAATCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATc                                            <  1:50920/249‑1 (MQ=255)
                                                                                                                                                                                                       gCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCg                                           >  1:564247/1‑251 (MQ=255)
                                                                                                                                                                                                                ggACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGAGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGAt                                  >  1:814289/1‑251 (MQ=255)
                                                                                                                                                                                                                  aCCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGAtt                                 >  1:838818/1‑250 (MQ=255)
                                                                                                                                                                                                                    ctgcctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATtat                               >  2:557288/1‑250 (MQ=255)
                                                                                                                                                                                                                     cgccgaagaAACCGTGGCCCCAGAGGAAAAAAAAAGCCCTGTCCGGAGAGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATaa                             <  1:300161/245‑1 (MQ=255)
                                                                                                                                                                                                                       ccggaaaaaTCGTGGCGCAAGAGCCAACAAAAAGCCCTGTCCGGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGt                           <  1:660581/245‑1 (MQ=255)
                                                                                                                                                                                                                             aaaTCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                                                                                                                                          >  1:651499/1‑118 (MQ=255)
                                                                                                                                                                                                                              aaTCGTTGCGCCAGAGCCAACAAAAAGACCGGTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAagcag                    <  1:517153/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gTTGCGCCAGAGGAAACAAAACGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGCAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTACGGCAGAAATACACGGATATAAGACACCTTCATGATCGCCCAGGGCTTATCAGTAAAGCAGCCa                 >  2:397157/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAATCTCTTTTCTgg    >  2:190607/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 acaaAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAATCTCTTTTCTGGCg  >  1:537426/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          
GCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAATCTCTTTTCTGGCG  >  NC_000913/2248186‑2248672

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