Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,570,457	0	C	T	33.3%	80.1 / 14.3	45	E275K (GAA→AAA)	hsdS	type I restriction‑modification system specificity subunit
Reads supporting (aligned to +/- strand):  ref base C (6/24);  new base T (15/0);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.57e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAG  >  NZ_CP009273/4570210‑4570703                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             tGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCtt                                                                                                                                                                                                                                                            >  2:42245/1‑249 (MQ=255)        aaGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTGTTTTCTTTTTGAACAttt                                                                                                                                                                                                                                                     >  1:282242/1‑249 (MQ=255)                ttCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAAc                                                                                                                                                                                                                                              <  1:379725/249‑1 (MQ=255)                        gTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCg                                                                                                                                                                                                                                     <  2:186838/250‑1 (MQ=255)                             tcatTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTATTTCA‑TTTTTTCACATTTTTGAAACAGAataacaatttg                                                                                                                                                                                                                                >  1:96342/1‑242 (MQ=255)                                    aTTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTATTCACTTTTTGTAAAATTTATAAAACGAATAACGtttttttttaca                                                                                                                                                                                                                         >  1:1257/1‑242 (MQ=255)                                                        aaaaaaaTATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACt                                                                                                                                                                                                   <  2:172912/250‑1 (MQ=255)                                                           aaaaTATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGTGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAAt                                                                                                                                                                                                >  1:243263/1‑250 (MQ=255)                                                           aaaaTATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTGTTTTATTTTTTAACATTCTTTCAAACACAAATTGTTTTTTTTT‑ACATGT‑CCAGCACTCCCAGAACaaa                                                                                                                                                                                                 >  2:291462/1‑249 (MQ=255)                                                              atatTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCg                                                                                                                                                                                             <  1:290579/250‑1 (MQ=255)                                                                     gATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTTTTTACTTTTTTAAAATTTTAGAAAACCAAAATCGGTTTTATCT‑TCATGGCCCACCAC‑TACAGAACTAATGCGTAGTAtt                                                                                                                                                                                      >  1:575185/1‑250 (MQ=255)                                                                               ggTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑ATTAGTTTACTTTTTTAAAATTTTCTAAACCGCATATTGTTTTTTTCT‑ACATGG‑TCCGCACGTGCAGAACTCATGCGTAATATTGAATgcatc                                                                                                                                                                             >  2:182812/1‑245 (MQ=255)                                                                                            tGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTGTTTACTTTCTGAAAAATCTATAAACCGAAAATTGTTTTTATTT‑TAATTG‑CCAAGAAGAACAAAACTACAGCGTAGGATTTGATGA‑GCCAAAC‑CACATTTAtt                                                                                                                                                             >  1:237991/1‑250 (MQ=255)                                                                                                     gAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTtga                                                                                                                                                    <  2:300759/250‑1 (MQ=255)                                                                                                         tCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTTTTTACTTTTTTTACATTTTTGAAAACGCAAATCGTTTTTATTT‑ACATTG‑CCAGCACGTACAAAAATAAAGAGGAGTATTTGGTGA‑ACAACAC‑CAACTTTATTTTGCttttttga                                                                                                                                                 >  1:505622/1‑244 (MQ=255)                                                                                                              tAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑Acc                                                                                                                                          <  2:494100/250‑1 (MQ=255)                                                                                                              tAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑CTTTGTTCTCCTTCTTAAAACTCTAGAAAACCCAAATAGCTTTTTTCTAAAATGG‑CCA‑AACCCACAGAAACAAAATGTAGGATTGGGTGT‑CCAAAAA‑CACCTTCAATTTGCTTTTATGaaaaaa                                                                                                                                            >  1:447719/1‑246 (MQ=255)                                                                                                                     aTCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGaaa                                                                                                                                               <  1:535044/239‑1 (MQ=255)                                                                                                                           atatAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑CTTTGTTTACTTTTTTAAAATTCTAGAAAAAGAAAATCTTTTTTTTAT‑TAATGT‑CCCACACGAACAAAAACAATACGAATAATTTGTTTA‑ACAACACACACTTCAATTTGTC‑TTGATGAAAA‑ACCAATACAGAAttt                                                                                                                              >  2:521199/1‑248 (MQ=255)                                                                                                                               aGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTGTTTACTTTTTTAACATTTTTTATAACGAATAATGTTTTTTTTT‑TCATTT‑CCCGCCCCTTCAGAAAAAATACTGAGTATTGTATGG‑TCCACAC‑CCACTTCATTTTGCTTTTTTGAAAG‑AAAATTACTTAATTCCGtt                                                                                                                          >  2:210394/1‑249 (MQ=255)                                                                                                                                             ttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCaaa                                                                                                           <  2:54201/249‑1 (MQ=255)                                                                                                                                             ttgtAATTATTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCaaa                                                                                                           <  2:62689/249‑1 (MQ=255)                                                                                                                                              tgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTGTTCACTTTTTGAACATTTTAGAAAACGAATAACGTTTTTATCT‑TCATTG‑CCCGCACGTACAAAAATAATACGGAGGATTTGAAGA‑ACAACAC‑CACCTTCATTTTGCTTTGATGAAAGAACAATTAAG‑AAATCAATTAAA‑AAAGATTCaaaaa                                                                                                         >  1:204469/1‑249 (MQ=255)                                                                                                                                               gtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTGTTTACTTTTTTTAAAATTTAGAAAACGAAAATTGTTTTTATTT‑ACATGG‑CCCACACGGACAAAAAAAAAACGTGATATTTGGTTA‑ACAACAA‑CACATTTCTTTTGCTTTTTTTAAAG‑ACCAATAAGTAAATTAATTAAG‑ATAAATTTaaaaat                                                                                                        >  1:258337/1‑248 (MQ=255)                                                                                                                                                 aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTATTTCACTTTTTTAACATTTTAAAAAACGAATATCGTTTTTATCT‑TCATTG‑GCAGCACGTACAAAACTAATGCGTAGTATTTGATGA‑CCAACAC‑CACTTTCATTTTGCTTTTATGAAAA‑ACCATTAAGTAAATTAGTTAAA‑ATAGATTTAAAAatt                                                                                                       >  1:266192/1‑249 (MQ=255)                                                                                                                                                                 cacaAACACCAACAAATTCTAAACTTCCGTTATAGCGATTAAATAAAAGATCTCCATCTTGTAAGTTGTGGCGG‑TTTAGTTCGCTTTCTGAACAGTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCGTAAATACAGa                                                                                      <  1:458620/250‑1 (MQ=255)                                                                                                                                                                  acaAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGaa                                                                                     <  2:246092/250‑1 (MQ=255)                                                                                                                                                                       cacCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTg                                                                                <  2:517188/250‑1 (MQ=255)                                                                                                                                                                         ccAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCg                                                                              <  1:507761/250‑1 (MQ=255)                                                                                                                                                                         ccAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCg                                                                              <  2:321826/250‑1 (MQ=255)                                                                                                                                                                                 ttCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCaaaa                                                                      <  1:242953/250‑1 (MQ=255)                                                                                                                                                                                      aaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCTG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACg                                                                 <  1:447075/250‑1 (MQ=255)                                                                                                                                                                                             gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCAtttt                                                          <  1:477611/250‑1 (MQ=255)                                                                                                                                                                                              ttATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCAttttt                                                         <  1:291462/250‑1 (MQ=255)                                                                                                                                                                                                 taGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCt                                                       <  2:599749/249‑1 (MQ=255)                                                                                                                                                                                                  aGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTTTTTACTTTTTTAAAATTTTTTAAACCCAATATTGTTTTTTTTT‑TCTTGG‑CCCGCCCGTACAAAACTAAATCGTAATATTTGATTG‑CCAAAAC‑CACCTTCCTTTTGCTTTTTTTAAAA‑ACCATTTACTAATTTTATTAAA‑ATAGATTTAAAATTTTACTTTTTTAATACAAAATGTTGTGGCTCCAAATTACACCAttttttttt                                                     >  2:261025/1‑246 (MQ=255)                                                                                                                                                                                                  aGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑CTTTGATTATTTTTTATACAATCTAGAAACCCAAAAATGTTTTTTTTT‑TAATAG‑CCACAACGTACAGAAATAATAAGTCGTTATGGGTGAGCCACAAC‑CCTTTTCTTTTGTTTTTATTGAAGA‑ACCTTAATTTAATTTA‑TTTAG‑ATAGAATCCAAATTTTACTTTTTTTATACAGAATGTTTATGCCCAAAATTAGATCATTTTTCttt                                                     >  1:5146/1‑248 (MQ=255)                                                                                                                                                                                                        taaataaaAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAAttt                                               <  2:265947/250‑1 (MQ=255)                                                                                                                                                                                                        taaataaaAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAAttt                                               <  2:434940/250‑1 (MQ=255)                                                                                                                                                                                                            taaaAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTTTTTACTTTTTTTAAATTTTTTAAAACCAATAACGTTTTTTTTT‑TCATTG‑CCCGCACCGACAGAAATAAATCGTAGTAATTGTTGA‑CCAACAC‑CACCTTCTTTTTGCTTTTTTTGAAA‑AACAATACTGTATTTAGTTAAGAAAAAATTAAAAATTTAA‑TTCTTTAATAAAAAAAGTTTGGGGCCAAAATTACTTCCTTTTTTTTTCTATTTTTCa                                           >  2:321848/1‑250 (MQ=255)                                                                                                                                                                                                              aaaGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCAtt                                         <  1:112482/250‑1 (MQ=255)                                                                                                                                                                                                                 gATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAAc                                      <  1:42245/250‑1 (MQ=255)                                                                                                                                                                                                                                         gg‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGAc               >  2:365648/1‑249 (MQ=255)                                                                                                                                                                                                                                                 ttCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGt       <  1:202913/250‑1 (MQ=255)                                                                                                                                                                                                                                                       ttCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAg  >  1:473036/1‑249 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             TGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCT‑ACATGG‑CCAGCACGTACAGAACTAATGCGTAGTATTGGATGA‑CCAACAC‑CACTTTCATTTGGCTTTGATGAAAG‑ACCATTACGTAATTCAGTTAAG‑ATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAG  >  NZ_CP009273/4570210‑4570703

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