Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,570,451	0	C	T	36.0%	145.1 / 26.1	75	E277K (GAA→AAA)	hsdS	type I restriction‑modification system specificity subunit
Reads supporting (aligned to +/- strand):  ref base C (2/46);  new base T (27/0);  total (29/46)
Fisher's exact test for biased strand distribution p-value = 2.21e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGAC  >  NZ_CP009273/4570210‑4570690                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   tGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCt                                                                                                                                                                                                                                             <  1:497924/249‑1 (MQ=255)  gACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTg                                                                                                                                                                                                                                            <  1:46783/249‑1 (MQ=255)    ccAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCAACTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTCTCttttttaa                                                                                                                                                                                                                                          >  2:47160/1‑244 (MQ=255)             gttTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTACTTTGTGGCGGTTTAGATTACTTTTTTAACaattttta                                                                                                                                                                                                                                 >  1:663670/1‑242 (MQ=255)              ttttCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGaa                                                                                                                                                                                                                                <  2:615272/249‑1 (MQ=255)                ttCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATAATTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAAc                                                                                                                                                                                                                              <  1:748652/249‑1 (MQ=255)                 tCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAAcc                                                                                                                                                                                                                             <  1:595886/249‑1 (MQ=255)                           catcatTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGtt                                                                                                                                                                                                                   <  1:135619/249‑1 (MQ=255)                            atcatTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGGTAATCGATCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGttt                                                                                                                                                                                                                  <  2:565400/249‑1 (MQ=255)                             tcatTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTCCTTTTTTAAACATCTATAAAACCAAAATTGttttt                                                                                                                                                                                                                >  2:327689/1‑249 (MQ=255)                                tttcATGTCGTGCTTAGGGTGATGAAAAAAATATGTCTATATATTCTGGTAAAGCATCTTTGTTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTTATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAATATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTGTCTGAACATTCTATAAACCGAATATCGTTTTGATc                                                                                                                                                                                                             <  2:435584/247‑1 (MQ=255)                                       tCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATg                                                                                                                                                                                                       <  2:252477/249‑1 (MQ=255)                                       tCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATg                                                                                                                                                                                                       <  2:151290/249‑1 (MQ=255)                                        cGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGGAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATgg                                                                                                                                                                                                      <  2:569051/249‑1 (MQ=255)                                         gTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTCTTAACATTCTAAAAACCAAATATCGTTTTGATTTACATGGcc                                                                                                                                                                                                    >  2:326998/1‑250 (MQ=255)                                         gTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGcc                                                                                                                                                                                                    <  1:666823/250‑1 (MQ=255)                                            cTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGc                                                                                                                                                                                                 <  1:450665/250‑1 (MQ=255)                                             tGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTCTTAACATTCTATAAAACGAATATCGTTTTTATTTACATTGCCAGCa                                                                                                                                                                                                >  1:487755/1‑250 (MQ=255)                                               aGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCAc                                                                                                                                                                                               <  1:51835/249‑1 (MQ=255)                                                        aaaaaaaTATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAAATCTAGAAAACGAATATTGTTTTGATTTACAAGGCCAGCAcagacaaaaa                                                                                                                                                                                      >  1:386698/1‑240 (MQ=255)                                                         aaaaaaTATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTTATTTTTTTACATTTTATAAAACCAATATATTTTTTATTTTCATTGCCAGCACGGACAGaaaa                                                                                                                                                                                     >  1:477623/1‑247 (MQ=255)                                                          aaaaaTATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTTTTTACATTCTATAAAACGAATATCGTTTTTTTCTTTATTGCCAGCACGTACAGAACTa                                                                                                                                                                                    >  2:216247/1‑249 (MQ=255)                                                             aaTATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAATTTACTTTCTGAACACTTTAGAAACCCAAAATCGTTTTTATTTACATGGCCCACACGGAAAGaaaaaaagt                                                                                                                                                                                >  1:716026/1‑243 (MQ=255)                                                              atatTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTTTTAACATTCTATAAACCCAATATTGTTTTTATTTATATTGCCAGCACGGACAAAACTAATACg                                                                                                                                                                               >  1:42263/1‑250 (MQ=255)                                                               tatTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTATTTACTTTTTGAAAATTTTATAAAAAGAAAATCGTTTTTATTTACATTGCCAACCCCTTAAGAACTaaattg                                                                                                                                                                               >  1:677822/1‑245 (MQ=255)                                                               tatTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgt                                                                                                                                                                              <  1:533879/250‑1 (MQ=255)                                                                      atatatTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTg                                                                                                                                                                      <  2:416140/250‑1 (MQ=255)                                                                          atTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTTCTTTTTTAACATTCTATAAACCCCCTTTTGTTTTTATTTTCATTGCCACCCCCTACCGAAACAATGACTAG‑TAttttatg                                                                                                                                                                  >  2:73433/1‑245 (MQ=255)                                                                            tCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTCCTTTTTTTACATTCTAAAAAACCAAAATAGTTTTTTTTTATATTGCCAACACGTACAGAAATAAATCGTAG‑TATTGGGTGAc                                                                                                                                                                >  1:756058/1‑250 (MQ=255)                                                                            tCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTCCTTTTTTAAAATTTTAAAAAACGAATATCGTTTTTATTTTCATGGCCAGCACCTACAAAACTAATTGGGAG‑TATTGGATggc                                                                                                                                                                >  1:466381/1‑248 (MQ=255)                                                                             cTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTTTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGAc                                                                                                                                                                <  2:605827/249‑1 (MQ=255)                                                                                gTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCaa                                                                                                                                                             <  2:705420/249‑1 (MQ=255)                                                                                 tAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAaca                                                                                                                                                           <  1:464257/250‑1 (MQ=255)                                                                                       aTCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACtt                                                                                                                                                    <  2:580725/250‑1 (MQ=255)                                                                                            tGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACttt                                                                                                                                                   <  1:325498/246‑1 (MQ=255)                                                                                                  aTCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACT‑CCACTTTCATTTGGCtt                                                                                                                                         <  2:617119/250‑1 (MQ=255)                                                                                                   tCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCttt                                                                                                                                        <  1:100608/250‑1 (MQ=255)                                                                                                            aaTAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGa                                                                                                                               <  2:261109/250‑1 (MQ=255)                                                                                                             aTAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTATTTACTTTTTGAACATTCTAGAAAACAAAAATTGTTTTTTTTTATATTGCCAGCACGGACAAAAATAAAACGTAG‑GATTTGATTAGCAAAA‑CCCCCTTTATTTTGATTTTTTGaaaaaa                                                                                                                              >  2:219074/1‑247 (MQ=255)                                                                                                                          gATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTc                                                                                                                <  1:307701/250‑1 (MQ=255)                                                                                                                            tataGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTGAAAATTCTAAAAACCAAATATCGTTTTTTTTTTCATTGCCCCCCCCAACAGAAATAATACGTAG‑TATTTTATAACAAAAA‑CCCCCTTCTTTTTGCTTTTTTTAAAAAACATTTAGTAA‑TTTAg                                                                                                              >  2:212994/1‑250 (MQ=255)                                                                                                                              taGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGtt                                                                                                            <  2:304467/250‑1 (MQ=255)                                                                                                                                gCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TAtt                                                                                                                                                                       <  1:55741/191‑1 (MQ=255)                                                                                                                                      tttGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTATTTTTTTAAAATTCTACAAAACCAATATTTTTTTTTTTTTCATTGCCCGCAAGATCAACAATAATTCGTAA‑GATTATATTGTCAACA‑CAACCTTCCTTTTGATTTTTTGAAAAACAATTATGTGA‑TTAATTTaataaata                                                                                                    >  2:186840/1‑244 (MQ=255)                                                                                                                                        tGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTTTGAAACTTTTAGAAAACCAAAATTTTTTTTTTCTTAATTGCCAGCACGGACACAACTAATGAGGAGTTATTGG‑TGACCAAAA‑CCACTTTCATTTTTCTTTTTTTAAAAAACATTACGAAATTTTAGTTAAGA‑AGAtt                                                                                                  >  1:769115/1‑250 (MQ=255)                                                                                                                                               gtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTATGAAAAATTTAGAAAACGAAAATAGTTTTTATCTAAATAGACAGACCGCACAGAAAAAATTCGTAG‑GATTGGCTGAGCAAAA‑CCACTTTCATTTGTCTTTGATGAAAAAACATTACGTAA‑TTCAGTTAATATAGAATCtaaaat                                                                                           >  2:648342/1‑248 (MQ=255)                                                                                                                                               gtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAAt                                                                                            <  1:720542/249‑1 (MQ=255)                                                                                                                                                          aaTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTaa                                                                                <  2:386730/250‑1 (MQ=255)                                                                                                                                                          aaTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTaa                                                                                <  1:342257/250‑1 (MQ=255)                                                                                                                                                            tAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTAATTTTTTTAAATTCTTACAAAAGAAAATAGTTTTGTTTTTAATTGGCCGAAAGGACAAAAATAAAAACGAA‑TATTTGGTGGACAACA‑CCACCTTCATTTTTCTTTTTTTGAAGAACAATATGTAA‑TTTACTTAATATATATTTAAAATTTatattatttaa                                                                               >  2:674401/1‑239 (MQ=255)                                                                                                                                                               cccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACa                                                                           <  1:713132/250‑1 (MQ=255)                                                                                                                                                                   caAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTCCTTTTTTAACATTCTAGAAAACGCATATTGTTTTTATCTACATGGCCAACACGTACAAAACAAATTCGTAG‑TTTTGGGTGACCAAAA‑CCACATTTATTTTGCGTTTTTTAAAAACCATTTAGGAA‑TTCAATTTAGATAAATTCAAAATTTTAATTCTTTAATaaaaaaa                                                                       >  1:27894/1‑244 (MQ=255)                                                                                                                                                                   caAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAAt                                                                       <  1:264545/250‑1 (MQ=255)                                                                                                                                                                           aacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTCTGAACATTCTAGAAACCGAATATCGTTTTTATCTACATTGCCAGCACGTACAAAAATAATGCGGAG‑TATTTGGTGACCAACA‑CCACTTTCATTTTGGTTTTGTGAAAGACCATTACGTAA‑ATTAGTTAAGATAGATTTAAAAATTAAATTTTTTAATAAAAAATGtttgggc                                                               >  1:769168/1‑245 (MQ=255)                                                                                                                                                                               aaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGC‑TTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCaaa                                                          <  2:327691/250‑1 (MQ=255)                                                                                                                                                                                   cTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTATTTTACTTTTTTAACATTTTTGAAAACGAATATCTTTTTTATCTACATTGCCCGCACCTACAGAACTAAAACGTAG‑TATTTGATGACCCACA‑CCACTTTCATTTGGCTTTTATTAAAGACCATTACGTAA‑TTTAGTTTAGATAGATTTAAAATTTAACTTCTTTAATACAGAATGTTGCGGCTCAAAAt                                                        >  1:729219/1‑249 (MQ=255)                                                                                                                                                                                      aaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTTCTTTTTTTAAATTTTTGAAAACGAAAATCGTTTTTTTTTTCATTGCCCGCCCGTACAGAACTAATGCGTAG‑TATTTGGTGACCAACA‑CCCCTTTTATTTTGTTTTTTTTAAAAACCAATACATAA‑TTTAGTTAAGATTAGTTTAAAAATTAAATTCTTTAATACAAAATGTTTTGGCTCAAAATTaaa                                                    >  1:290702/1‑248 (MQ=255)                                                                                                                                                                                      aaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTAc                                                     <  1:92965/249‑1 (MQ=255)                                                                                                                                                                                      aaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTAc                                                     <  1:23502/249‑1 (MQ=255)                                                                                                                                                                                       aCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGc                                                   <  1:64890/250‑1 (MQ=255)                                                                                                                                                                                        cTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGcc                                                  <  1:140608/250‑1 (MQ=255)                                                                                                                                                                                             gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTCTTAACATTCTAGAAAACGAATATCGTTTTTATCTACATTGCCAGCACGTACAGAACAAAAGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGGTTTGATGAAAAACCAATACGTAA‑TTTAGTTAAGATAGATTTTAAAATTAAATTCTTTAATACAAAAAGTTTCGGCTCAAAATTACACCAtttt                                             >  1:125395/1‑250 (MQ=255)                                                                                                                                                                                             gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCAtttt                                             <  2:386682/250‑1 (MQ=255)                                                                                                                                                                                              ttATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGCTTTGTTTTCTTTCTCGACATTCTAAAAAACGAATATCTTTTTTTTCTATATAGCCAGCACGAACAGAACCAATTGGTAG‑TATTTGTTGATCAACA‑CCCCTTTCTTTTTGCTTTTTTGATAAAACATTACGCAA‑TTCATTTCAGATAGAATCTAAATTTTACTTTTTAAATAAAAAATATTGCGGCTCAAAATTACCCCCgtttt                                            >  1:309842/1‑250 (MQ=255)                                                                                                                                                                                                 taGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCt                                          <  1:187555/249‑1 (MQ=255)                                                                                                                                                                                                 taGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTg                                         <  1:629081/250‑1 (MQ=255)                                                                                                                                                                                                   gCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTCTTTACATTTTAAAAAACGCATATCGTTTTTATTTTCATTGCCAGCACGGAAAAAAAAAAATCGTAG‑TAATGGATGACCCACACCCACTTT‑ATTTGGCTTTTATGTAAGAACATTACGTAA‑TTTAGGTTAGATAAATTTAAAAATTAACTTCTTTAATACAGAATGTTTTGGCGCAAAAATACCCCAttttttttt                                        >  1:380628/1‑245 (MQ=255)                                                                                                                                                                                                       gTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAAtt                                   <  1:644997/250‑1 (MQ=255)                                                                                                                                                                                                             aaaaGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTTTTTACATTCTTTAAAACAAATATTGTTTTTTTTTTCATTGCCAGCCCGTACACAAATAAATAGTAG‑TATTTGATTAACAACA‑CAACTTTCCTTTTGGTTTTATTTAAGACCCTTACGTAA‑TTTCGTTTATATAGATTTAAAATTTTAATTTTTTAATACAAAATTTTTTGGGTCAAAATTTAACCCTTTTTTTTTCTATTTTTCAt                             >  2:577530/1‑250 (MQ=255)                                                                                                                                                                                                                 gATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAAc                         <  1:149389/250‑1 (MQ=255)                                                                                                                                                                                                                       cATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACt                        <  2:693186/245‑1 (MQ=255)                                                                                                                                                                                                                        aTCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGccc                  <  1:114604/250‑1 (MQ=255)                                                                                                                                                                                                                              tAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAAt            <  1:506121/250‑1 (MQ=255)                                                                                                                                                                                                                                 tttGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATAcc         <  1:327689/250‑1 (MQ=255)                                                                                                                                                                                                                                         ggTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACATACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGAc  >  1:91679/1‑249 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   TGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAG‑TATTGGATGACCAACA‑CCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAA‑TTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGAC  >  NZ_CP009273/4570210‑4570690

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