Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,394,8180TC34.3% 136.8 / 13.8 70D625G (GAT→GGT) csrDRNase E specificity factor CsrD
Reads supporting (aligned to +/- strand):  ref base T (44/2);  new base C (0/24);  total (44/26)
Fisher's exact test for biased strand distribution p-value = 9.26e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTACAAATCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAC  >  NZ_CP009273/3394572‑3395065
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
gTTACAAATCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGAGTAAGGACATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcca                                                                                                                                                                                                                                                       >  1:668420/1‑249 (MQ=255)
     aaaTCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGc                                                                                                                                                                                                                                                 >  1:19673/1‑250 (MQ=255)
        tctcTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCg                                                                                                                                                                                                                                              >  1:492680/1‑250 (MQ=255)
             cccTGAGTCCAGAATGGCTTAATGCACCAGCCGCATTCCCCGCCTGTTGGCAACCGCGTAAGTACATTCATCTGCATATGTGCGCAAACAATCTGGCGAGCAGACAAGCACACTCCCATGAAACGCAGCGCGCATTATTCAACTTAAAAACGGATAAAAGGGGAGGTTAAACGGATTATTTTTGTAAATTTTTTTTCCCTTTTTTTTCAGGTGGTTGTGAGGACGCGAAAAATCCCCCTTGCCCGCCTGt                                                                                                                                                                                                                                         <  2:533084/249‑1 (MQ=255)
              ccTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCt                                                                                                                                                                                                                                           <  1:158567/247‑1 (MQ=255)
                  aGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAac                                                                                                                                                                                                                                      <  1:129405/248‑1 (MQ=255)
                               ttAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAg                                                                                                                                                                                                                       >  1:556823/1‑250 (MQ=255)
                                     cATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTg                                                                                                                                                                                                                 >  1:601694/1‑250 (MQ=255)
                                      aTCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGc                                                                                                                                                                                                                >  2:237175/1‑250 (MQ=255)
                                           aagctttcaccGCTTGTTTGCAGCCGCGAAAGTACATTCATCTGCATATGTGCGGCAACATTCGGGCGAGCAAACAGGCCCACTCCCATGAGACGAAGCGCGAATTTTTCTACTTAAAAACGGTTTAAACGGAAGTTAAAACCAATTATTTTGGTAAATATTTTTTCGCGTTAGTATAAAGGGGCGGGGAGGACGCAAAAAACCCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTc                                                                                                                                                                                                           <  2:740002/240‑1 (MQ=255)
                                             gcaTTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGc                                                                                                                                                                                                         >  2:79399/1‑250 (MQ=255)
                                               aTTCCTCGCCTGTTTGCAACGGCGTAAGGTCATCCATCTGCATATGTGCTGAAACAATCTGGCGAGCAGACAAGCACACCCCCATGAGCCGCAGCGCTTATTATTCTACGTGAAAACGGATAAAACGGCAGTTTAAACCGAGTATCTTTGTGAATATTTTTCCACGTTAGTATCAGGGGGCTGGGAGGACGCAAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCtt                                                                                                                                                                                                       <  2:121818/250‑1 (MQ=255)
                                                ttCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTc                                                                                                                                                                                                      >  2:180093/1‑250 (MQ=255)
                                                               cAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCgg                                                                                                                                                                                       >  1:163923/1‑250 (MQ=255)
                                                                   cgcgAAAGGTCATTCATCTGCATAGGTGCTCCACAAACCTGGCGAGCAAACAAGCACACCCCCATGACACCCAGCGCGCTTTTTTCTACGTAAAAACGGATTAAACGGCAGTTTAACCCGAGTATCTTTGTAAATATTTTTTCACTTGAGTATCAAGGGGCTGGGAGGCCGGGAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGc                                                                                                                                                                                   <  2:12896/250‑1 (MQ=255)
                                                                    gcgTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCg                                                                                                                                                                                  >  2:582367/1‑250 (MQ=255)
                                                                    gcgTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCATTGCCCGCCTGTAACACCGCGCTGGATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTggag                                                                                                                                                                                  >  2:582402/1‑248 (MQ=255)
                                                                                tcatCTGCATATGTGCTGCAACAATCGGGCGAGCAGACAAGCACCCCCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGAATTAAACGGCAGTTTAACCCGAGTATCTTTGTGAATTTTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGt                                                                                                                                                                       <  1:613329/249‑1 (MQ=255)
                                                                                 catCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGc                                                                                                                                                                     >  1:152943/1‑250 (MQ=255)
                                                                                 catCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGATTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGc                                                                                                                                                                     >  2:611153/1‑250 (MQ=255)
                                                                                       cATATGTGCTGAAACAATCGGGCGACCAGAAAAGCACACACCCAGGAGACGCAGCGCGCTTTTTTCTACGTGAAAACGGATTAAACGGCAGGTTAACCCGAGTATCTTTGTGAATATTTTTTCACGTTATTATAAAGTGGCTGTGAGGCCGCGAAAAACCCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTc                                                                                                                                                                <  2:707013/249‑1 (MQ=255)
                                                                                       cATATGGGCTACAACAACCGGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGAGAAAACGGATTAACCGGCGGTTTAAACCGAGTATCTTGGTAAATATTTTTTCAGGTTAGTATAAAGTGGCTGTGAGGACGCAAAAAACCCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTgg                                                                                                                                                                  <  2:246280/247‑1 (MQ=255)
                                                                                           tgtgCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTCAACCTGGGTGCTGGTCCCggc                                                                                                                                                           >  2:456915/1‑249 (MQ=255)
                                                                                           tgtgCTGCAACAATCGGGCAACCAGACAAGCACCCCCCCATGAGCCGCAGCCCGCTTTATTCTCCGTAAAACCGGTTAAACCGGCAGGTTAACCCGAGTATCTTTGTAAATATTTTTTCACTTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCgg                                                                                                                                                            <  2:656453/249‑1 (MQ=255)
                                                                                              gctgcAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg                                                                                                                                                          >  1:427332/1‑248 (MQ=255)
                                                                                               ctgcAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAg                                                                                                                                                       >  2:682739/1‑250 (MQ=255)
                                                                                                           ggCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAgg                                                                                                                                           >  2:236298/1‑250 (MQ=255)
                                                                                                                      aaGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcag                                                                                                                                >  1:217914/1‑250 (MQ=255)
                                                                                                                      aaGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcag                                                                                                                                >  2:92524/1‑250 (MQ=255)
                                                                                                                         cacaCTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagcag                                                                                                                             >  1:334031/1‑250 (MQ=255)
                                                                                                                           cccTCCCATGAGACGCAGCGCGCATTATTCTACGTAAAAACGGATTAAACGGCAGGTTAAACCGATTATCTTTGTAAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGCCGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCt                                                                                                                           <  1:301330/248‑1 (MQ=255)
                                                                                                                           cacTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGACGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCt                                                                                                                           >  2:139213/1‑250 (MQ=255)
                                                                                                                                ccATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGt                                                                                                                        <  1:194895/248‑1 (MQ=255)
                                                                                                                                  aTGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTtctc                                                                                                                    >  1:80873/1‑250 (MQ=255)
                                                                                                                                       aCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCCTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTc                                                                                                               >  2:204889/1‑250 (MQ=255)
                                                                                                                                           aGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCtt                                                                                                           >  1:63142/1‑250 (MQ=255)
                                                                                                                                             cgcgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTtct                                                                                                         >  1:296849/1‑250 (MQ=255)
                                                                                                                                                gcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCaa                                                                                                      >  2:671720/1‑250 (MQ=255)
                                                                                                                                                 cATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAAt                                                                                                     >  1:121896/1‑250 (MQ=255)
                                                                                                                                                        cTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGCGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCt                                                                                               >  1:611357/1‑249 (MQ=255)
                                                                                                                                                          aCGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGAc                                                                                            >  1:691660/1‑250 (MQ=255)
                                                                                                                                                           cGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATAAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACAACGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGAc                                                                                            >  2:502023/1‑249 (MQ=255)
                                                                                                                                                           cGTAAAAACGGATTAAACGGCAGGTAAACCCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGTTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGAca                                                                                           <  1:312964/250‑2 (MQ=255)
                                                                                                                                                              gAAAACGTATTAAACGGAAGGTAAATCGGAGTTTCTTTGTGAATTTTTTTTCAGTAGAGTATAAAGGGCCGTTGAGCCCGCAAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGc                                                                                        <  2:405734/250‑1 (MQ=255)
                                                                                                                                                               aaaaCGTATAAAACGAAGTGTAAACCCGATTTTTTTTGTAATTTTTTTTCCCGTTTATTACAAGTGTGTTGTGAGGCCGGAAAAAACCCCCTTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGcc                                                                                       <  2:608275/250‑1 (MQ=255)
                                                                                                                                                                 aaggaaTTAACGGGCAGTTTAACCGGGGTATTTTTGTAAATTTTTTTCCCCGTTAGTAGAAGGTGGCTGGGAGGACGCAAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGcccc                                                                                     <  2:560447/245‑1 (MQ=255)
                                                                                                                                                                 aaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGcccc                                                                                     >  1:257707/1‑250 (MQ=255)
                                                                                                                                                                       ttAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATg                                                                                >  1:62566/1‑249 (MQ=255)
                                                                                                                                                                       ttAAACGGAAGTTTAAACCTAGTATCTTTGTAAATATTTTTTCCCTTTATTTTCAAGTGGCTGTGAGGACGCGAAAAATCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATgg                                                                               <  1:455426/250‑1 (MQ=255)
                                                                                                                                                                         aaaCGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGCGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGGTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAg                                                                             >  2:247027/1‑250 (MQ=255)
                                                                                                                                                                             ggCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCtt                                                                          >  1:25542/1‑249 (MQ=255)
                                                                                                                                                                              gCAGGTAAAACGGAGTATTTTTGTAAATTTTTTTTCACGTTATTATAAAGGGGCTGGGAGGCCGCAAAAAAACCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGa                                                                        <  2:675187/250‑1 (MQ=255)
                                                                                                                                                                                 ggTTAAACCGAGTATCTTTGTAAATATTTTTTCGCGTTATTATCAAGTGGCTTTGGGGCCGGAAAAAAACCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTa                                                                     <  1:48448/250‑1 (MQ=255)
                                                                                                                                                                                     aaaCGGATTATTTTTGTAAATTTTTTTTCGCTTTATTAACAAGTGGCTGGGAGGACGAGAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTc                                                                 <  2:511401/250‑1 (MQ=255)
                                                                                                                                                                                              aTCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAg                                                        >  1:465561/1‑250 (MQ=255)
                                                                                                                                                                                               tCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGt                                                       >  2:655624/1‑250 (MQ=255)
                                                                                                                                                                                                cTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGACTCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGCTTCTGACCAGCCCCGGATGGTGCTTGAGTAACCCAACATTAAGtt                                                      >  1:63507/1‑250 (MQ=255)
                                                                                                                                                                                                cTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGtt                                                      >  1:315340/1‑250 (MQ=255)
                                                                                                                                                                                                   tgtaAATATTTTTTCCCGTAAGTACAAGGGGGCGGTGAGACCGCAAAAAACCCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCt                                                    <  2:42831/249‑1 (MQ=255)
                                                                                                                                                                                                        atatTTTTTCACGTTAGTTTAAAGTGGCGGTGAGGACGCAAAAAACCCCCCTGCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATc                                              <  2:45636/250‑1 (MQ=255)
                                                                                                                                                                                                        atatTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATc                                              >  2:260665/1‑250 (MQ=255)
                                                                                                                                                                                                         tatTTTTTCGCGTTATTATAAGGTGGCTGTGAGGCCGCAAAAAACCCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATcc                                             <  2:644960/250‑1 (MQ=255)
                                                                                                                                                                                                           ttttttCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCCGCCCCGGATGGAGCTTGAGTAACTAAACATTACGTTCTTTGTTCCaa                                           >  2:732714/1‑250 (MQ=255)
                                                                                                                                                                                                              tttCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTg                                        >  2:409793/1‑250 (MQ=255)
                                                                                                                                                                                                                    gTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACtt                                  <  1:359591/250‑1 (MQ=255)
                                                                                                                                                                                                                           taaaGGGGGTGTGAGGACGCAAAAAACCCCCCTCGCCCGCCTGTAACACCGCGCTGAAGCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCATCg                           <  2:174135/248‑1 (MQ=255)
                                                                                                                                                                                                                             aaGTGGCTGGGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTc                         <  1:112456/250‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAAc                     >  2:593021/1‑248 (MQ=255)
                                                                                                                                                                                                                                        gAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCt               >  1:734693/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     caTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTCAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAc  >  2:577619/2‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
GTTACAAATCTCTGCCTGAGTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAC  >  NZ_CP009273/3394572‑3395065

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