Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,560,495	0	C	T	33.3%	129.1 / 16.8	72	T301T (ACG→ACA)	dosC	diguanylate cyclase DosC
Reads supporting (aligned to +/- strand):  ref base C (9/39);  new base T (24/0);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 4.85e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAG  >  NZ_CP009273/1560249‑1560743                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          aGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGt                                                                                                                                                                                                                                                         >  1:874545/1‑249 (MQ=255)              aTAAATTCATCGCCCCCGTAGCGGAAAACATAATCTCTTCTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATa                                                                                                                                                                                                                                            <  1:938030/249‑1 (MQ=255)                      aTCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTc                                                                                                                                                                                                                                    <  2:751714/249‑1 (MQ=255)                       tCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGCTTAATAATATTTCCAGTTAATCCATACCCACCTCGt                                                                                                                                                                                                                                  >  1:931505/1‑250 (MQ=255)                         gCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTg                                                                                                                                                                                                                                 <  1:318507/249‑1 (MQ=255)                         gCCCCCGTAGCGGAAAACATAAGCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTgg                                                                                                                                                                                                                                <  1:96662/250‑1 (MQ=255)                          cccccGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGAGCTCTTTGAATTTATCAACGTCAACAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGgc                                                                                                                                                                                                                               <  2:122047/250‑1 (MQ=255)                               gTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTTTCAGTAAATCCATACCGACTTCGTGGggcgg                                                                                                                                                                                                                           >  2:776255/1‑248 (MQ=255)                                         aCATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCaa                                                                                                                                                                                                                <  2:738329/250‑1 (MQ=255)                                            tAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCCTACCGACTTCGTGGCGCGCTACTTCTTCAAAc                                                                                                                                                                                                              >  2:176043/1‑249 (MQ=255)                                                         ggACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAATAATTTTGTCAATAAATCCATAACCACCTCGGTGCGCGATACTTTTTTAAACaaattaacccgg                                                                                                                                                                                                  >  1:350756/1‑238 (MQ=255)                                                           aCGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAgg                                                                                                                                                                                              <  1:416422/250‑1 (MQ=255)                                                           aCGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATGTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATGCTTCTTCAAACAATTCACGCAGCAAgg                                                                                                                                                                                              <  1:437938/250‑1 (MQ=255)                                                             gTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTTGTTAGGAAATCCATACCGACTTCGGGGGGGGGTACTTTTTTTAAAAAATTACGCAGCAAggg                                                                                                                                                                                             >  1:368320/1‑248 (MQ=255)                                                              ttGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGtaa                                                                                                                                                                                           <  2:56965/250‑1 (MQ=255)                                                               tGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTCGTCCGTACAACCATAACGAATTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGgaa                                                                                                                                                                                           >  1:493572/1‑249 (MQ=255)                                                                      aaaGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGGAAAATCGTCAGTAAATTCATAAAGACATAGTGGGGCGCGAATTCTTCAAACAAATTACGCAACAAGGGCATCATTTgt                                                                                                                                                                                   >  2:348409/1‑249 (MQ=255)                                                                       aaGGCCGGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCg                                                                                                                                                                                  <  2:892384/250‑1 (MQ=255)                                                                             tGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATCCATAACGACATCGTGGGGCGGTACTTCTTTAAACAATTCACCCAGCAAAGTAATAATTTTCGATACgg                                                                                                                                                                            >  1:200823/1‑250 (MQ=255)                                                                                gagaCTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTAt                                                                                                                                                                         <  2:368141/250‑1 (MQ=255)                                                                                gagaCTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTAt                                                                                                                                                                         <  1:920114/250‑1 (MQ=255)                                                                                     tttACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATACATACCGACTTTGTGGCGCGGTACTTATTCAAAAAATTTACGCAGCAAGGGAATAAATTGCGATACAGGATTTCtt                                                                                                                                                                    >  1:392964/1‑250 (MQ=255)                                                                                       tACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTTAGTACATTCATACCCACTTCTTGGGGCGCGACTTTTTTCAACAAATTACCCCACAAGGTAATAATTTGCGATACGGGATTTCtttt                                                                                                                                                                  >  2:272929/1‑248 (MQ=255)                                                                                           cAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATCCATAACGACTTCGTGGCGCGGGACTTTTTCAAACAATTCACGCAGCACAGTAATAATTTGCGATACGGGATTTCTTATCTg                                                                                                                                                               >  1:579414/1‑249 (MQ=255)                                                                                              aaTTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGtaat                                                                                                                                                           <  2:988382/250‑1 (MQ=255)                                                                                                tttCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGtaat‑aa                                                                                                                                                        <  1:255277/250‑1 (MQ=255)                                                                                                tttCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGGTTAGGAAATTTTTCAGGACATCCATACCGACTTTGTGGCGAGGTACTTCTTCAAAAAATTTTCGCAACAAAGTAATAATTTGTGGTACGGTATTTTTTTTCTGtaat‑a                                                                                                                                                         >  2:959676/1‑249 (MQ=255)                                                                                                           gtgtTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTTGTCAGTACATCCATACCGACTTCGTGGGGCGGTACTTTTTCTAACAATTTACGCAACAAAGGAATAATTTTCGATAAGGGATTTTTTATTTTTTAT‑AAAAAattttccc                                                                                                                                             >  2:234583/1‑245 (MQ=255)                                                                                                                 tGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGa                                                                                                                                        <  1:557339/249‑1 (MQ=255)                                                                                                                             tCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCa                                                                                                                           <  1:886837/250‑1 (MQ=255)                                                                                                                                ttGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATCCCTAACGACTTCGTGGGGCGGTAATTCTTTAAACAATTTACGCAGCAAGGTAATAATTTGTGATAAGGTATTTTTTATTTGTTAT‑AAAAAATTTAACCCCAGAATTCTGTTATTCAAAt                                                                                                                        >  2:334150/1‑250 (MQ=255)                                                                                                                                 tGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGAAATTTCGTCAGTAAATTCATAACGACTTCCTGGGGCGATAATTCTTTTAACAAATCACCCAGAAAGGGAATAATTTGCGATACGGTATTTCTTTTCTGTTAT‑AAAAATTTTAACCCCAGACTATTTTTTTTCAAAtt                                                                                                                       >  2:278529/1‑250 (MQ=255)                                                                                                                                   aTCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCTGTACATCCCAACCGACTTCGTGGGGCGGTACTTTTTCAAACAATTCACGCACCCAGGTAATAATTTTCGATTCGGTATTTCTTTTCTGTAAT‑AAAAATTTTACCCGCCGACTTCTGTTTTTTTAAtttt                                                                                                                     >  2:118129/1‑249 (MQ=255)                                                                                                                                    tctcTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTTGTTAGTACATCCATAACGACTTCCTGGGGCGGTACTTCTTTCAACAAATCAAGCCACAAAGGAAAAATTTTCGGTAAGGTATTTTTTTTTTTTTAT‑AAAAAATTTACCCGCCAAATTTTGTTTTTTTAAttttt                                                                                                                    >  2:945620/1‑248 (MQ=255)                                                                                                                                        tttGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTgg                                                                                                                  <  1:292958/248‑1 (MQ=255)                                                                                                                                                 aTCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCTGTAAATCCATAACGCCTTCGTGGCGCGAGACTTCTTTAAAAAATTTAACCAGCAAGGGAATAATTTTCGATAAGGTATTTCTTATTTGTAAT‑AAAAATTTTAAACGCCAACTTCTGTTATTTAAATTTTTTGTGTTTCgtttg                                                                                                       >  2:764461/1‑246 (MQ=255)                                                                                                                                                      cGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTTTCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTTTGGTGTTACGCATGGTttg                                                                                                  >  1:418455/1‑250 (MQ=255)                                                                                                                                                          aataatCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTttgattg                                                                                              <  2:268776/250‑1 (MQ=255)                                                                                                                                                          aataatCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTttgattg                                                                                              <  1:331826/250‑1 (MQ=255)                                                                                                                                                            taatCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATCCATACCGACTTCGTGGCGCGGTACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTTCGATAAGGTATTTTTTATTTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTTAAATTTCTGGTGTTACGCATGGTttgattga                                                                                             >  1:779005/1‑249 (MQ=255)                                                                                                                                                            taatCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGaa                                                                                            <  2:5231/250‑1 (MQ=255)                                                                                                                                                                 aGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTCTTTGTACATCCATAACGGCTTCGTGGGGCGAAACTTTTTCAAACAATTCACCCAGCAAAGGAATAATTTTCGATAACGGATTTTTTTTCTTTAATAAAAAATTTCCCCCCCA‑ACTTCTGTTTTTTTAATTTCTTGTGTTTCACACGGTTTGTTTTTAAATAc                                                                                       >  2:440440/1‑250 (MQ=255)                                                                                                                                                                  gCACTGACAGCGTTGTACCGGTCCGGTTGGCATGGGCAAGTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATAcc                                                                                      <  2:725844/250‑1 (MQ=255)                                                                                                                                                                   cACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCg                                                                                     <  2:236500/250‑1 (MQ=255)                                                                                                                                                                    aCTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGGAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTTAAATTTTTTGTGTTACGCATGGGTTGATTTGAAATACCGt                                                                                    >  1:64043/1‑250 (MQ=255)                                                                                                                                                                    aCTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGGAATTTTGTCAGTACATCCCTAACGACTTCGTGGGGCGGGACTTCTTCCAAAAAATCCCGCCGCAAAGGAATAATTTTCGATACGGGATTTCTTTTCTTTAAT‑AAAAAATTTACCCCCCAAATTCTGTTATTTAAAATTTTTGTGTTTATCATGGTTTGTTTTGAAATTCCGt                                                                                    >  2:683507/1‑250 (MQ=255)                                                                                                                                                                                 tACCGGTACGGTTGGCATGGGCAATATCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGaa                                                                        <  1:807544/249‑1 (MQ=255)                                                                                                                                                                                  aCCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATCCATACCGACTTCGTGGGGCGGTACTTCTTCAAACAATTCAAGCAGCAAGGGAATAATTTGCGATAAGGTATTTTTTATCTGTAAT‑AAAAAATTCACCCGCAGACTTCTGTTATTCAAAATTCTGGGGTTACGCATGGGTTGATTGAAAATACCGTCGAAATCCTGAAt                                                                       >  1:211945/1‑249 (MQ=255)                                                                                                                                                                                     ggtccggtTGACATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGGACATCCATAACGACTTCGTGGGGCGATACTTTTTTAAACAATTTACGCAGCAAGGGAATAATTTTCGATACGGGATTTTTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTTAAATTTTTTGGGTTTTGCATGGGTTGTTTTAAAATACCGTTGAAATCCTGAAttaaa                                                                   >  1:79785/1‑246 (MQ=255)                                                                                                                                                                                            ttGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTTGTCAGGACATCCATACCGACCTCGTGGGGCGGTAATTCTTTTAAAAAATTACGCAGCAAGGTAATAATTTTCGATACGGTATTTTTTTTCTTTAAT‑AAAAATTTTACCCCCAGACTTCTGTTTTTTTAATTTTTTGTGTTATGCATGGGTTTTTTTTAAATACCGTCGAAAACCTGAATCAAACgg                                                                >  2:277577/1‑246 (MQ=255)                                                                                                                                                                                               gCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGata                                                          <  2:910291/249‑1 (MQ=255)                                                                                                                                                                                                  tGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGc                                                      <  2:132304/250‑1 (MQ=255)                                                                                                                                                                                                    ggCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGcc                                                     <  2:874545/249‑1 (MQ=255)                                                                                                                                                                                                     gCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTGACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTa                                                   <  1:682942/250‑1 (MQ=255)                                                                                                                                                                                                             gcgTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGc                                            <  1:985803/249‑1 (MQ=255)                                                                                                                                                                                                              cgTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTTGTCAGTACATCCATACCGACTTCGTGGGGCGGTACTTCTTTAAAAAATTTACGCAGCAAGGGAATAATTTTCGATACGGTATTTTTTTTTTGTAAT‑AAAAAATTCACCCCCAGACTTCTTTTATTTAAATTTTTTGTGTTACGCATGGTTTGATTGAAAAAACCGTCGAAATCCTGAATTAGACGGGGGGTATGGCCTACTTCTGCa                                           >  2:394513/1‑249 (MQ=255)                                                                                                                                                                                                               gTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGGACATCCATAACGACCTCGTGGCGCGATACTTCTTCAAAAAATTTCCGCAGCAAGGTAATAATTTTCGATACGGGATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGAATTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCaa                                          >  1:736349/1‑249 (MQ=255)                                                                                                                                                                                                               gTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAAt                                         <  2:4201/250‑1 (MQ=255)                                                                                                                                                                                                               gTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAAt                                         <  1:944874/250‑1 (MQ=255)                                                                                                                                                                                                                      gATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTa                                  <  2:64195/250‑1 (MQ=255)                                                                                                                                                                                                                       aTAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTaa                                 <  2:368320/250‑1 (MQ=255)                                                                                                                                                                                                                             ggTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAaataa                            <  2:891736/249‑1 (MQ=255)                                                                                                                                                                                                                                 ggAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCg                       <  1:452014/250‑1 (MQ=255)                                                                                                                                                                                                                                       ggcggTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCttt                 <  2:579784/250‑1 (MQ=255)                                                                                                                                                                                                                                       ggcggTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCtt                  >  2:196872/1‑249 (MQ=255)                                                                                                                                                                                                                                        gcggTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTa                <  2:112532/250‑1 (MQ=255)                                                                                                                                                                                                                                         cggTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCCCGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTAt               >  2:130568/1‑250 (MQ=255)                                                                                                                                                                                                                                           gTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATgg             >  2:18106/1‑250 (MQ=255)                                                                                                                                                                                                                                            ttAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGt            <  2:418455/250‑1 (MQ=255)                                                                                                                                                                                                                                                   aTTTCGTCAATACATCCATACCGACTTCGTGGCGCGATCCTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGAGACGGTATTGCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGTTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACca     <  2:931505/250‑1 (MQ=255)                                                                                                                                                                                                                                                    tttCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACcac    <  2:754997/250‑1 (MQ=255)                                                                                                                                                                                                                                                      tCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAGTGTCGACCTTTATGGTTAAACCACAg  >  1:586587/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          AGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAAT‑AAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAG  >  NZ_CP009273/1560249‑1560743

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