Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0610GA25.6% 88.7 / 11.3 43intergenic (‑124/+79)ves/spyHutD family protein Ves/ATP‑independent periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base G (29/3);  new base A (0/11);  total (29/14)
Fisher's exact test for biased strand distribution p-value = 6.33e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTT  >  NC_000913/1824813‑1825307
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
gatgCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaa                                                                                                                                                                                                                                                        <  1:137637/251‑1 (MQ=255)
                    gTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                                                                                                                                    >  2:115060/1‑251 (MQ=255)
                             aCGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGAAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGACCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCCACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGCGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGAATTACTg                                                                                                                                                                                                                           >  1:539127/1‑251 (MQ=255)
                               gTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                                                                                                                          >  1:125944/1‑250 (MQ=255)
                                   tGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGAGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACATCTGTTTGAGAAGAGAAAAGACAACAGCCGATCCTGTCCACCGCCTTACTGCCAGGt                                                                                                                                                                                                                     >  2:50020/1‑251 (MQ=255)
                                       ggTGGGCACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTGCGCCACAGATTCACCGACATTTTACGCATATCAAAGAATTCCATACTCGCTCCCGGCCCGTTCTTTTAATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACAGCATTACTGCACGGTAGTg                                                                                                                                                                                                                 >  2:487294/1‑251 (MQ=255)
                                       ggTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg                                                                                                                                                                                                                 >  1:121361/1‑251 (MQ=255)
                                       ggTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg                                                                                                                                                                                                                 >  1:61874/1‑251 (MQ=255)
                                       ggTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg                                                                                                                                                                                                                 >  1:295001/1‑251 (MQ=255)
                                           ggAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTCCGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAAAAAAGAAACCCCCCGACCTTGTCCACCGCTTTACTGCaggaaaggggac                                                                                                                                                                                                              >  2:54007/1‑239 (MQ=255)
                                            gAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATATTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGTGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                                                                                                                             >  2:2934/1‑250 (MQ=255)
                                                   cAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                                                                                                                                     >  1:330447/1‑251 (MQ=255)
                                                                 ttCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                                                                                                                               >  1:399803/1‑243 (MQ=255)
                                                                         cagctttcccccACAGATTCCCCGATATTTCACGCATAACAAATTTTTCCATCCTCGCTCCCGGCTTGTTCTTTTCATTGTTATCGTTCAGGTACCTGGGGTTGGTTATCGGCAACTGGGGTGATTTACGTTAGTGTTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                                                                                                                               <  1:40586/241‑1 (MQ=255)
                                                                                 cGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCACGGCCCGTTCTTTTCATTGATCTCGTTCAGGTATCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGCGGTGTTCCGGCTTTGTTTACAACTGTTTGAGGAGAGAAGAGAAAACCGCCGATCCTGTACACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGACTTATGGTTTTTgtgtgaacgagtctt                                                                                                                                                                       >  2:298159/1‑237 (MQ=255)
                                                                                             aCCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGGTTAAAAATTATTctgcaggttccg                                                                                                                                                          >  2:631331/1‑245 (MQ=255)
                                                                                               cGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGCAggg                                                                                                                                                        >  1:358503/1‑250 (MQ=255)
                                                                                                          gCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                                                                                                                             >  1:94108/1‑251 (MQ=255)
                                                                                                               acAAATTATTCCATCTCCGCTCCCGCCCTTTTTTTTTTATTGTTATCGTTAAGGTACCCGTGGATGGTTATCGGCAACTCGGTTTATTTACTTTGGTGGTGATCAGCCTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCg                                                                                                                                        <  1:279516/250‑1 (MQ=255)
                                                                                                                     tATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGCTTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTTTTCCGCAGTTGCAGGAATTTTAACTTTTGCCGCTGGa                                                                                                                                   >  1:546689/1‑250 (MQ=255)
                                                                                                                     tATTCCATACTCGATCCCAGCACGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGCTTTACGCTAGTGGTGTTCAGGCTTTTTTTACAACTGTTTGAGAAGCGAAAAGAAAACCGCCGATCCTGTCAACCGCATTACTGCACGGTAGTGGACAAGACCGGCGGCATTCAGGTTTTTTACTTCAAGGTTATTGCACAGTTGCAGGCATTTTACCTTTTTCCGCTggg                                                                                                                                   >  2:484177/1‑249 (MQ=255)
                                                                                                                                  cctcccggcctttttttttttttGTTATCGTCCAGGTACCCGGGGATGTTTTTCGGCATCTGGTGTGATTTACTTGAGGGGTAATGAGGTTTTTTTTACATTGTTTTGAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                                                                                                                     <  1:356575/230‑1 (MQ=255)
                                                                                                                                          cccttttttttCATTGTTTGCGTTCGGGTAGCGGTGGTGGGTTTCCGGCACCTCGGGTGATTTATGTTGGGGGGGATCGGCTTTTTTTTACAATTTTTGGAAAAGAGAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCa                                                                                                             <  1:120361/244‑1 (MQ=255)
                                                                                                                                                     cATTGTTATCGTTGGGGTAACGGTGAATGTTTATCGGCAACTCGGGTTTTTTACGTATGTGGTAATCAGTTTTTTTTTACACCTTTTGAAGAAGAAAAAAAAAAACCGCCGATCCTGGCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCAtt                                                                                                  <  1:163473/251‑1 (MQ=255)
                                                                                                                                                              tgttggggTAACGGTGTATGGTTACCGGCAATTCGGGTGATTTGCGTTGGTGGTGAGCAGGTTTTTTTTAAAATTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCt                                                                                          <  1:177393/244‑1 (MQ=255)
                                                                                                                                                                ttCAGGTACCGGTGGATGTTTCTCGGCACCGGGGGTGATTTACTTTAGTGGTAATCAGTTTTTTTTAACATTTGTGGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCttt                                                                                        <  1:520442/250‑1 (MQ=255)
                                                                                                                                                                    ggTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTGCGCTGGACGTTTTGTCCGACGCTTCTCAAAATTTGCATTAAATTGCtttttttt                                                                                   >  2:398718/1‑248 (MQ=255)
                                                                                                                                                                                 ggttaccggaatctggggTGTTTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTGTTGAAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGAATTTTACCTTTTGCCGATGGAAGTTATGTGCCAGGCAGCGCCGAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAg                                                                       <  1:582885/235‑1 (MQ=255)
                                                                                                                                                                                           aaCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGGTTTTTTGCTGAACGGATATTTAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGGCGTTCTTTCAGACGCTTTTCCAAATTATCCTTAAATTTCTTTTTCTTTTTCGGCGTCCGGAtgttgaaa                                                             >  2:151909/1‑247 (MQ=255)
                                                                                                                                                                                                  gTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGATGGACGTTCTGTCAGACGCTTCTCAAAATTATCATTAAATTGCTTTTTCTTTTCCGGCGTCAGGCTGTTGTTAATTTTGtt                                                     >  2:134614/1‑251 (MQ=255)
                                                                                                                                                                                                    gttttcggttagtgggcaccagtcTTTTTTTACAACTTTTTGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTc                                                    <  1:343173/227‑1 (MQ=255)
                                                                                                                                                                                                        tctgtgaggggggatcgtctttttttaACATCTGTTAGAGAAGAAAAAAAAAAACCGCCGATCCTGGCCACCGCATTACTGCAAGGGAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTggg                                                <  1:586333/230‑1 (MQ=255)
                                                                                                                                                                                                        tACGTTAGTGGTGATCAGGCTTTTTTTACAAATGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGAAGGCATTTTACCTTTTTCCGCTGGACGTTCTTTCAGACGCTTCTCCAAATTAGCATTAAATTGCTTTTTCTTTTCCGGCGTCCGGATGTTGTAAATTTTTTTCTGGGt                                               >  1:171394/1‑251 (MQ=255)
                                                                                                                                                                                                           gTTAGTGCTGAACAGGATTTTTTTCCAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAAACGACCGGCGGATTTAAGTTTTTTGGGTGAAAGAATAATCAGTAATTGCAGGCATTTTACCTTTTGCCGCCGGAAGTTCGTCCAAAAGCCTCTCAAAAATTGCATTAAATTGCTTTTTCTTTTCACGCGTCAGGCAGATGTAATTTTTGtttttgttttt                                            >  2:536041/1‑242 (MQ=255)
                                                                                                                                                                                                               gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCCGAAAAATTATTCAGCAGTTGCAGGCCTTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAAATTTGTTCTGGGTTTTCAtt                                        >  1:287396/1‑250 (MQ=255)
                                                                                                                                                                                                                         gggtttttttttcatttGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTACATGTGCGCCAGc                               <  1:247611/235‑1 (MQ=255)
                                                                                                                                                                                                                              tttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTTCCGGCATTTTATCTTTTTCCTCTGGGCGTTTTGTCAGACGCTTCTCAAAATTAGAATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTTTCAATTTTGTTTTGGGTTTCTATGTGCGCCAGCCTGttt                         >  2:55981/1‑250 (MQ=255)
                                                                                                                                                                                                                              tttttAAAACTGTTTGAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGtt                          <  1:596441/250‑1 (MQ=255)
                                                                                                                                                                                                                                  tACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCtt                     >  1:285530/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaaaCTGTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCttt                    <  1:631331/249‑1 (MQ=255)
                                                                                                                                                                                                                                    aCAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCCTTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCttt                    >  1:192082/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aaCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTgc                  >  1:429176/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gtttagaaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCt               <  2:285530/243‑1 (MQ=255)
                                                                                                                                                                                                                                                aaaaaaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttct          <  1:50020/243‑1 (MQ=255)
                                                                                                                                                                                                                                                   aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGGTGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTcc       >  1:100663/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCAtttt  >  1:277927/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
GATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTT  >  NC_000913/1824813‑1825307

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