Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,401,4520GA28.3% 91.0 / 13.4 53F627F (TTC→TTTcsrDregulator of CsrB and CsrC decay
Reads supporting (aligned to +/- strand):  ref base G (28/10);  new base A (0/15);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 5.23e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTCCTTTGTTACAAATCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTC  >  NC_000913/3401204‑3401683
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
gTTTCCTTTGTTACAAATCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGa                                                                                                                                                                                                                                        >  2:147532/1‑250 (MQ=255)
        tGTTACAAATCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATccc                                                                                                                                                                                                                               >  1:225451/1‑251 (MQ=255)
           tACAAATCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcccatt                                                                                                                                                                                                                            >  1:202327/1‑248 (MQ=255)
            aCAAATCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTg                                                                                                                                                                                                                           <  2:61306/251‑1 (MQ=255)
                    ctGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTg                                                                                                                                                                                                                   >  2:179478/1‑251 (MQ=255)
                        cTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAac                                                                                                                                                                                                               >  1:125035/1‑251 (MQ=255)
                          gAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCACCCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGTTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGGGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAac                                                                                                                                                                                                               <  1:235852/249‑1 (MQ=255)
                            gTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcg                                                                                                                                                                                                           >  1:507522/1‑251 (MQ=255)
                            gTCCAGAAAGGCTTAATGCATCAGCAGCATTCCCCGCCTGTTTGCAACCGCTTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGACG‑‑AGCAGACAAGCACACCCCCATGAGACGCAGCGCGCTTTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATTTTTTTTCACGTTAGTATCAAGTGGCTGGGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACcg                                                                                                                                                                                                           <  1:522929/251‑1 (MQ=255)
                                             gcagcagcagCATTCCTCGCCTGTTTCCACCCGCGAAAGGTGATTCATCGGAATAGGTGCGGCAAAAATCGGGCG‑‑AGAAGACAAGCCCCCCCCCATGAGACGCAGCGCGCATTTTTCTACGTGAAAACGGATGAAGCGGCAGGTAAACCCGATTATTTTGGTAAATATTTTTTCCCGTTAGTATCAAGTGGCTGTGGGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGTGTCTg                                                                                                                                                                                          <  2:4189/247‑1 (MQ=255)
                                             gCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTg                                                                                                                                                                                          >  2:325845/1‑251 (MQ=255)
                                                         ttCCTCCCCTGTTTGCAACCGCTAAGGGTCATTCTTCTGCATATGTGCTGCAACAA‑‑TGGCGGTAGCAGACAAGCAACCACCCAGGATACGCAGCGCGCATTTTTCTACGAGAAAACGTATAAAGCGGCGGTTAAACCGGGGTATTTTTGTAAATATTTTTCCCCGTTAGTATCAGGTGGCTTTGGGGGCGCAAAAAAATCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTc                                                                                                                                                                               <  2:1432/250‑1 (MQ=255)
                                                                      tGCACCCGCTTAAGGTCATTCATCTCCATATGTGCTGCAACAAGTTGGCG‑‑AGCAGAAAAGCACACTCCAAGGAGCCGCAGCGCGCATTATCCTGCGTGAAAGCGGATTAAACGGCGGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATGAATGGGTCGGGAGGACGCAAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGcc                                                                                                                                                                  <  1:206185/250‑1 (MQ=255)
                                                                        cAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGt                                                                                                                                                               >  1:416297/1‑251 (MQ=255)
                                                                         aaCCGAGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTg                                                                                                                                                              >  2:435278/1‑251 (MQ=255)
                                                                                 aaGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAAc                                                                                                                                                      >  1:270603/1‑251 (MQ=255)
                                                                                 aaGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAAc                                                                                                                                                      >  2:312367/1‑251 (MQ=255)
                                                                                     ntcattaatCTGCATATGTGCTGAAACAATCTGGTG‑‑AGCAGACAAGCACAATCACATGAGACGAAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTAACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATACACTTGCCCGCCTGTAACACAGCGCTGAATCAGGGTCTGACAATCGCTTCGCGAACGCAAGCCGGTGGCGTAAACCTg                                                                                                                                                   >  2:247989/2‑251 (MQ=255)
                                                                                          tcatCTGAATATGTGCGGCAACAATCGGGCG‑‑AGCAGACAAGCACCCCCCCATGAGACGCAGCGCGCATTATTCTACGTGAAACCGGATTAACCGGCAGTTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAGGTGGCTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTg                                                                                                                                               <  1:332002/250‑1 (MQ=255)
                                                                                                  atatGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGCAAACCTGGGTGCTGGTCCCg                                                                                                                                      >  2:168524/1‑251 (MQ=255)
                                                                                                      gtgcgccaacaaTCTGGCG‑‑AGCAGACAAGCACACCCCCATGAGACGCAGCGCGCATTATTCTACGTAAAACCGGATAAACCGGCAGTTTAAACCGAGTATTTTTGTGAATATTTTTTCCCTTTATTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCTTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg                                                                                                                                   <  1:555527/244‑1 (MQ=255)
                                                                                                           gcAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAggg                                                                                                                                                                                              >  2:97180/1‑186 (MQ=255)
                                                                                                           gcAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAggg                                                                                                                                                                                              <  1:97180/186‑1 (MQ=255)
                                                                                                                     ggCG‑‑AGCAGACAAGCACACTCCCAGGAGACGCAGCGCGCATTATTCTACGTAAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTGGTAAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGCCGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAgg                                                                                                                    <  1:422714/250‑1 (MQ=255)
                                                                                                                           ataaGACAAGCACACCCCTATGAGACGCAGCGCACATTATTCTACGTGAAAACGGATAAAGCGGAAGGTAAAACCGAGTATTTTTGTGAATATTTTTTCCCGTTATTATCAAGTGGCTGTGGGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                                                                                                               <  2:160326/248‑1 (MQ=255)
                                                                                                                           aGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                                                                                                               >  1:93517/1‑251 (MQ=255)
                                                                                                                           aGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                                                                                                               >  1:443112/1‑251 (MQ=255)
                                                                                                                                 aaaGCACACTCCCATGAACCGCAGCGCTCATTTTTTTAGGTGAAACCGGATTAAACGGCAGGTAAAACCGATTATTTTTGTGAATATTTTTCCACGTTAGTAGAAAGGGGCTGTGAGGACGCAAAAAACCCCCCTCGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcag                                                                                                         <  2:49048/250‑1 (MQ=255)
                                                                                                                                  aaGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagc                                                                                                        >  1:357078/1‑251 (MQ=255)
                                                                                                                                  aaGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagc                                                                                                        >  2:285741/1‑251 (MQ=255)
                                                                                                                                  aaGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGACTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagc                                                                                                        >  1:99692/1‑251 (MQ=255)
                                                                                                                                       cacTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTACACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTg                                                                                                   >  1:35800/1‑251 (MQ=255)
                                                                                                                                        acTCCCATGAGACGCAGCGCGCATTATTCTACGAGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATACCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGACCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCCGCTgg                                                                                                  >  2:534363/1‑251 (MQ=255)
                                                                                                                                              aTGAGACGCGGCGCGCATTATTCTCCGTAAAAACGGATTAAACGGCAGGTTAAACCGATTATCTTTGTAAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGCCGCAAAAAAACCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTcc                                                                                            <  1:510770/251‑1 (MQ=255)
                                                                                                                                                                     tACGTGAAACGGGATTAAACGGCAGGTTAACCCGGTTATTTTTTGGAATATTTTTTCACGTTATTAGCAAGTGGGGGTGGGGCCGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGa                                                                      <  2:186177/250‑1 (MQ=255)
                                                                                                                                                                            aaagttataaaaCGGCAGTTTAAACCGAGTTTCTTTGTGAATTTTTTTTCCTGTTATTATCAAGTGGTTGGGAGGGCGCAAAAAAACCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGcccc                                                              <  1:465385/242‑1 (MQ=255)
                                                                                                                                                                                ggATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGAt                                                          >  1:548617/1‑251 (MQ=255)
                                                                                                                                                                                ggATAAAACGGCAGGTTAAACCGATTATTTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGAt                                                          <  1:247989/251‑1 (MQ=255)
                                                                                                                                                                                   ttAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATgg                                                        >  2:40313/1‑250 (MQ=255)
                                                                                                                                                                                      accGGCAGGTTAAACCGAGTATTTTTGTGAATATTTTTTCAGGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCt                                                    <  1:264452/249‑1 (MQ=255)
                                                                                                                                                                                        agggcagttaaaacgggttttttttGTAAATTTTTTTTCCCGTTAGTATCAAGTGGCGGGGGGACCGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCtt                                                   <  2:270603/229‑1 (MQ=255)
                                                                                                                                                                                            aGGTTAACCCGAGTTTCTTTGTGATTTTTTTTTCACTTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTtct                                                                                              <  1:378684/203‑1 (MQ=255)
                                                                                                                                                                                            aGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTtct                                                                                              >  2:378684/1‑203 (MQ=255)
                                                                                                                                                                                                 naaCCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTc                                          >  2:246126/2‑250 (MQ=255)
                                                                                                                                                                                                 aaacccttttatttttGTAAATATTTTTTCTTGTTATTATCAAGTGGCTGTGAGGACGCGAAAAACTCCCCTGCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGTTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACt                                           <  1:496693/237‑1 (MQ=255)
                                                                                                                                                                                                 aaaCCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACATGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTc                                          >  1:371408/1‑250 (MQ=255)
                                                                                                                                                                                                                gtggAATATTTTTTCATGTGAGTACCAAGTGGGTGTCAGGCCGAGAAAAATTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCtt                            <  2:35800/247‑1 (MQ=255)
                                                                                                                                                                                                                   aaagttttttccccGTTATTATCAAGTGGCTGTGAGGACGCAAAAAAGCCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGAt                        <  1:59687/238‑1 (MQ=255)
                                                                                                                                                                                                                       atTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCa                     >  2:534091/1‑250 (MQ=255)
                                                                                                                                                                                                                         tttttCACGTTATTATCAAGTGGCTGTGAGGACGCAAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAAc                   <  1:147532/250‑1 (MQ=255)
                                                                                                                                                                                                                             ccactttattATCAAGTGGCTGTGAGGCCGCAAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGt               <  1:325845/241‑1 (MQ=255)
                                                                                                                                                                                                                                    aGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTAcc        >  1:105225/1‑250 (MQ=255)
                                                                                                                                                                                                                                          aatgggggtggggggACGCAAAAAAACCCCCTCGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTc  <  2:125035/238‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
GTTTCCTTTGTTACAAATCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCG‑‑AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTC  >  NC_000913/3401204‑3401683

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: