Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0610GA31.1% 117.4 / 20.4 61intergenic (‑124/+79)ves/spyHutD family protein Ves/ATP‑independent periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base G (38/4);  new base A (0/19);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 2.98e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTC  >  NC_000913/1824818‑1825298
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
aaTGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                                                                                                                          >  1:143354/1‑251 (MQ=255)
   gCTGGCACGCCAGTAAAAAGCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATCCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACGGTTTGAGAAGAGAAAAGAAAAccgcc                                                                                                                                                                                                                                        <  2:115335/250‑1 (MQ=255)
     tGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGAt                                                                                                                                                                                                                                     >  1:689331/1‑251 (MQ=255)
                    aaTCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACCACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                                                                                                                      >  2:115980/1‑251 (MQ=255)
                    aaTCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc                                                                                                                                                                                                                         >  1:476483/1‑248 (MQ=255)
                       cACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTTATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCAGATCCTGTCCACCGCATTACt                                                                                                                                                                                                                   >  2:65938/1‑251 (MQ=255)
                        ccGTTTTGCGGGTGAGACCGTCAAAATTCCGCGTGTTCCCCCGGCACATTTCCGCCACAGATCCCCCGACATTTTACCATAATAAAAGTATTCCATACGCGCCCCCGGCCCTTTTTTTTTATTTTTACCGTTCAGGTACCTGTGGATGGTTATCGGCAACTGGGGTAATTTACGTTAGTGGTGATCGGGTTTTTTTTCCAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTg                                                                                                                                                                                                                  <  2:40509/250‑1 (MQ=255)
                          gTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAGCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                                                                                                                 >  2:269154/1‑250 (MQ=255)
                          gTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGGTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCAACCGCATTACTGc                                                                                                                                                                                                                 >  2:884190/1‑250 (MQ=255)
                                   gTGGAAACGTGAAAATTTCGCGCGTTCCCCGGGCAGCATTCCGCCACAGATCCACCGACATTTTACGCATACAAAAGTTTTCCATCCCCGCCCCCGGCCCTTTCTTTTCATTGTTATCGTTCAGGTAACCGGGGATGGTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGACCAGGTTTTTTTTACACCTGTTTGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                                                                                                                       <  1:36686/251‑1 (MQ=255)
                                       gAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                                                                                                                    >  2:646041/1‑250 (MQ=255)
                                        aaCGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                                                                                                                    >  1:611751/1‑249 (MQ=255)
                                         aCGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                                                                                                                    >  1:329352/1‑248 (MQ=255)
                                            tGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                                                                                                              >  1:329444/1‑251 (MQ=255)
                                                   ttCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATCTCAAAGTATTCCATCCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGCTGGTTATCGGCAACTCGGGTGAGTTACGTTCGTGGTGATCAGGCTTTTTTTACAAATGTTTGAGTAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGtgt                                                                                                                                                                                       >  2:145741/1‑249 (MQ=255)
                                                   ttCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTAACCGACATTTTACGGATATCAAAGTATTCCATACTCCCTCACGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGCTTTACGTTAGTGGTGATCAGGCTTTTTTTACAACCGTTTGAGAAGAGAAAAGAAAACCGCCGCTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTcc                                                                                                                                                                                       >  1:311726/1‑250 (MQ=255)
                                                     cgcgcgTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACACGACCGGCGGTCttt                                                                                                                                                                                     >  2:302064/1‑250 (MQ=255)
                                                     cgcgcgTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcggctttt                                                                                                                                                                                     >  2:765167/1‑246 (MQ=255)
                                                                     agcaTTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                                                                                                                     >  1:725308/1‑250 (MQ=255)
                                                                           tCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGCTGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGATTTAAGATTTTTTGGTGGAAGatttt                                                                                                                                                              >  2:618618/1‑249 (MQ=255)
                                                                             cgccacatgtccACCAACATTTTGCGTATATCAAAGTTTTCCATACTCGCTCCCGGCCCGTTTTTTCCATTGTTATGTTTCAGCTAACGTGGATTGTTTACGGCCACCTGGGTGGATTTAGTTAAGTGGTAATCGGGTTTTTTTTACTACTGTTGAAGAAAAAAAAAAAAATCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                                                                                                                             <  1:751176/240‑1 (MQ=255)
                                                                               ccACAAATTCACCGACATTTTACGCATATAAAATTTTTCCAACCCCCCCCCCGGCCCTTTTTTTTATTGTTTATCGTTCAGGTAACGGGGGATGGTTATCGGCAACGGGGGTGTTTTACGTTAGGGGTGATCAGGTTTTTTTTACAGTTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTca                                                                                                                                                           <  1:269154/251‑1 (MQ=255)
                                                                                              cATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTCGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtt                                                                                                                                            >  1:308637/1‑251 (MQ=255)
                                                                                               tttttACGCATATAAAAGTTTCCCATCCTCGCCCCCGGCCCGTTTTTTTATTTGTTATCGTTCAGTTACCCGTGGATGGTTATCGGCACCTGGGGTTTTTTACGTTGGTGGTGATAGGGTTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtt                                                                                                                                            <  1:646041/249‑1 (MQ=255)
                                                                                                     cGCATATCAAAGTTTCCCATCCCCGCTCCCGCCCCGTTCTTTTCTTTTTTATCGTTCAGTTACCCGTGTATGTTTATCGGCAACGCGGGTGTTTTCCGTTAGTGGTGATCAGGTTTTTTTTACACCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                                                                                                                     <  1:293011/251‑1 (MQ=255)
                                                                                                               aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACCGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg                                                                                                                            >  1:251610/1‑251 (MQ=255)
                                                                                                                       catccgcgcccccGCCCCTTTCTTTTCTTTGTAATCGTTCAGGTGACGGTGGTTGGTTATCGACAACGCGTTGATTTTATGTTTGTGGTGATCAGTCTTTTTTGACAACTGTTTGAAAAGAAAAAGGAAAACCGCCGATCCTGTCCACCGCACTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCATTTGCCGCTGGACGTTCt                                                                                                                    <  2:188318/241‑1 (MQ=255)
                                                                                                                         cactcccgcccGGCCTGTTTTTTTCATTTTTATCTTGCGGGTAACGGTGGAGGTTTAGCGGAACAGGGGGTGTTTTACGTTAGTGGTATACAGGCTTTTTTTAAAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTTAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                                                                                                                  <  1:4176/243‑1 (MQ=255)
                                                                                                                                 cccggcctttttttttCATGTTTATCGTTGGGGTACCGGTGGTTGGTTTTCGGCAACCGGGGTTATTACCTTTGGTGGTGAGCAGGTTTTTTTAACAATGTTTGAAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCt                                                                                                          <  2:143354/239‑1 (MQ=255)
                                                                                                                                   cGGCCCGTTCTTTTTATTTTTTTCTTTCAGGTAACCGTGGATGGTTATCGGAACCTGGGTTGTTTTACGTTAGTGGTGATCGGTTTTTTTTTACAACTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc                                                                                                        <  1:529471/251‑1 (MQ=255)
                                                                                                                                         ttcggttctaTTGTTATCGTTCAGGTACCGGGGGTTGTTTATCGGCACCTGGGGTGATTTCTTTTAGTGTTGATCGGGTTTTTTTAACAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa                                                                                                   <  1:98970/241‑1 (MQ=255)
                                                                                                                                                  aTTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTTTCAGACGCTtctccaa                                                                                                   >  2:98970/1‑238 (MQ=255)
                                                                                                                                                         ttcgttcaggaagccttggttttttATCGGCAACTGGGGTGATTTAGGTTGGGGGTGATCGTGTTTTTTTTACAACTGTTAGAGAAAAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc                                                                                             <  2:102455/218‑1 (MQ=255)
                                                                                                                                                         aTCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc                                                                                             >  1:102455/1‑241 (MQ=255)
                                                                                                                                                           gtttCAGGTACCGGGGGATGTTTATCGGAACCTGGGGTTATTTACTTTTGTGGTGATCAGTTTTTTTTTACACCTTTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCtt                                                                                 <  1:175818/249‑1 (MQ=255)
                                                                                                                                                           ggttcagttaccgggggATGGTTATCGACACCTCGGGTGTTTTACGTAAGGGGTGACCGGGCTTTTTTTAAAACTGTGGGAGAGAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCtt                                                                                 <  1:507844/236‑1 (MQ=255)
                                                                                                                                                            gTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGCGGTGATAAGGCTTTTTTTACAACTGTTTGAGAAGAGAAACGAAAACCGCCGATCCTGTCCACCGCATTACTGCCAGGTAGTGGACAAGACCGGCGGTCCTCAGGTTTTTGGCTGAAAGATTATTTAGCAGTTGCCAGCCTTTTACCTTTTGCCGCTGGTCGTTCTCTCAGACCCTTCTCAACATTAGCCTTTAATTGCttt                                                                                >  2:281205/1‑251 (MQ=255)
                                                                                                                                                              tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGGGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCCTTTTACCTTTTGCCGCTGGACGTTCTGTCACACGCTTTTCAAAATTAGCATTAAATTGCttttt                                                                              >  1:817585/1‑251 (MQ=255)
                                                                                                                                                                     aCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGCAGTGGACAAGACCGGCGGTCTTAAGCTTTTTGGCTGAAAGATTATTCAGAAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCCGAAGCTTCTCAAAATTTGCATTAAATTGCTTTTTTTGTTcc                                                                       >  1:404411/1‑251 (MQ=255)
                                                                                                                                                                        gTGGATGTTTATGGGCATCTGGGGTATTTTCGTTTAGTGGTGACCGGGCTTTTTTTAAAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGc                                                                    <  1:842251/251‑1 (MQ=255)
                                                                                                                                                                        gTGGATGGTTATCGGCAGCGCGGTTTATTTCCGTTAGTGGTGATCAGTTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGc                                                                    <  1:215168/251‑1 (MQ=255)
                                                                                                                                                                         tGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTTAGCAGTTGCAGGCCTTTTACCTTTTGTCGCTGGACGTTCTGTCAGACGCTTTTCAAAATTAGCTTTAAATTTCTTTTTCTGtttagggg                                                                   >  2:442422/1‑245 (MQ=255)
                                                                                                                                                                           gATGTTTATCGGCACTTCGGGTGATTTACGTTAGTGGTGATAAGGTTTTTTTTAAAACTGTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTc                                                                 <  1:320024/251‑1 (MQ=255)
                                                                                                                                                                              ggTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAgg                                                              >  1:52986/1‑251 (MQ=255)
                                                                                                                                                                               gTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTTTGTTCCGGCGTCAgg                                                              >  1:847783/1‑250 (MQ=255)
                                                                                                                                                                                     ggCAACTCGGGTGATTTACGTTAGTGGGGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTTCCTTTTGCCGCTTGACGTTCTGTCAGACGCCTCTCAAAATTAACATTAAATTGCCTTTTCTGTTCCGGCGTCCGGGtgttgt                                                       >  2:405286/1‑251 (MQ=255)
                                                                                                                                                                                                  gATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCCTTTGCCGCTGGACGTTCTGTCAGACGCTTTTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCt                                           >  1:340743/1‑251 (MQ=255)
                                                                                                                                                                                                    aTTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTg                                          >  1:536052/1‑251 (MQ=255)
                                                                                                                                                                                                      ttACGTTAGAGGTGATCAGGCTTTTTTTACAACTGTTTGAGAGGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCGGAAAGATTTTTCGGCAGTTACAGGCATTTTTCCTTTTGCCGCGGGACGTTTTGTCAGACGCTTCTCAAAATTAGCTTTAAATTGTTTTTTTTGTTCCGGCGTCAGGGTGTTGTAAATTTTGTTCTggg                                        >  1:798942/1‑251 (MQ=255)
                                                                                                                                                                                                           gTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCGGGTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTc                                    >  1:740803/1‑251 (MQ=255)
                                                                                                                                                                                                           gTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGACTTAAATTTTTTGGATGAAAAAATATTTAGCAGTTGCAGGCATGTTTCCTTTTGCCGGTGGGAGGTCTAGCAGACGGTTCTCAAAATTAGCATTACATTGCTTTTTCTGTTCCGGCGTTCTGATGTTTTTAATTTTTTTCTGGGtttt                                    >  2:488026/1‑250 (MQ=255)
                                                                                                                                                                                                                gagtgttgacaagttttttttaACAACCTTTGGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTcc                                   <  2:251610/235‑1 (MQ=255)
                                                                                                                                                                                                                    ggtgattggcttttttttACACTTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCAtgtg                              <  1:775729/240‑1 (MQ=255)
                                                                                                                                                                                                                          cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTTAATTGCTTTTTTTGTTCCGGCGGCAGGATGTTGTAAATTTTGTTCTGGGTTTTCATGTGCGcccg                       >  2:570584/1‑249 (MQ=255)
                                                                                                                                                                                                                          cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGcc                         >  1:511992/1‑249 (MQ=255)
                                                                                                                                                                                                                             gCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTTAGCCATTGCAGGCATTTTAACTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAAACTTAAATTGCTTTTTCTGTTCCGGCGTCAGGGTGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCAt                    >  2:645886/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tttttACAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGtt                 <  1:884190/250‑1 (MQ=255)
                                                                                                                                                                                                                                 ttttaACAACTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGtt                 <  1:818512/250‑1 (MQ=255)
                                                                                                                                                                                                                                    ttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTTGCCTAAAGATTATTCAGCAGTTTCAGGCCTTTTAACTTTTGCCGCTGGACGTTTTTTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTTAATTTTGTTCTGGGTTTCCCTTGTCGCCAGCATGTTCGCt             >  2:450941/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTTAATTGCTTTTTCTGTTCCGGCGTTAGGATGTTTTAAATTTTGTTCTGGGTTTCCATGTTCGCCAGCATGTTTGCtt            >  2:771478/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gtggaagaaaaaaaaaaaaacCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCGGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTTGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTg     <  2:740803/239‑1 (MQ=255)
                                                                                                                                                                                                                                                ggaaaagagaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttc  <  1:65938/247‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
AATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTC  >  NC_000913/1824818‑1825298

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