Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	377,402	0	G	A	24.2%	88.1 / 35.9	62	intergenic (+91/+133)	yaiL/frmB	DUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB
Reads supporting (aligned to +/- strand):  ref base G (23/24);  new base A (6/9);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-02
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTC  >  NC_000913/377153‑377640                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                aGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGc                                                                                                                                                                                                                                                >  2:21980/1‑251 (MQ=255)  gCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCg                                                                                                                                                                                                                                               <  1:344143/251‑1 (MQ=255)     aaGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTc                                                                                                                                                                                                                                            >  1:390703/1‑251 (MQ=255)     aaGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTAATGGGCGAAGTGTAGACCGGATAAGACGttt                                                                                                                                                                                                                                            >  1:179562/1‑250 (MQ=255)         tGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTAcgc                                                                                                                                                                                                                                        >  1:360623/1‑251 (MQ=255)          gaaatcGCCTTTATTCCCGCCAGCGTTGACGATAAAATTGCTCAGCGCGATGCCAGCAGTATATTTTTGCACAGCGCGCTAAGCGCCGAAGAGCAGGATGTAGTCGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACTTAGATAAGTCGTGCGCAGCTGCCTGATGCGACG‑CTATGCGCGTCTTAGCAGGTCTACTTTTTATGGGCGAAGTGTATGCCGGATACGGCGTCCAcgcc                                                                                                                                                                                                                                       <  2:429497/246‑1 (MQ=255)          aaaTACGCCAGTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTAcgcc                                                                                                                                                                                                                                       <  2:187715/250‑1 (MQ=255)              aCGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGCCTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCAt                                                                                                                                                                                                                                   >  1:80210/1‑251 (MQ=255)                gCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGCGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                                                                                                                 >  1:472513/1‑251 (MQ=255)                 ccATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCg                                                                                                                                                                                                                                <  1:146933/251‑1 (MQ=255)                        ccccGCCAGCGTTGCCGATAAAATTGCTCAGCGCGAGGCCAGCAGTATAGTGTTGCACAGCGCGCGCAGCGCCGAAGAGCAGGAGGAAGACGACCCGTAGGCCGACTTCAAAGTGCCTGATTATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCTTACGGTTTATGGGCGTAGTGTAGACCGGATAAGGCTTTCACGCCGCATGCTGCAGTCg                                                                                                                                                                                                                         <  2:472513/250‑1 (MQ=255)                          ccGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCATATGCTTGATGCGACG‑CGAGGCTCGTCTTAGCAGGCCTACGGTTTATGGGGTAAGTGTAGACCGGAGAAGACGTTGCCGCCGGATCCgccatttgtg                                                                                                                                                                                                                       <  2:88731/251‑10 (MQ=255)                          ccGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTg                                                                                                                                                                                                                       <  2:258020/251‑1 (MQ=255)                           cGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCatttgtgc                                                                                                                                                                                                                      <  1:299660/251‑8 (MQ=255)                           cGCCAGCGTGGCCGATAAAATTGCTCAGCGCGATGCAAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTAGGCGGGTCTTATCAGGCCGACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCAGCCTGCAGGCGTgc                                                                                                                                                                                                                      <  2:338595/251‑1 (MQ=255)                                     gCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTg                                                                                                                                                                                                           >  1:282364/1‑251 (MQ=255)                                           aaaaTTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCTCAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                                                                                                                                                                                                      <  2:264004/251‑1 (MQ=255)                                                  cTCAGCGCGATGCCAGAAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACAACCCGTATGCCGACTTCAAAGAGACTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGAACGGGTAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGAGGCCTGATGCGACGCtg                                                                                                                                                                                                >  2:153413/1‑249 (MQ=255)                                                                agcagTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAgg                                                                                                                                                                                >  1:28219/1‑251 (MQ=255)                                                                  cagTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGcc                                                                                                                                                                              <  1:258884/251‑1 (MQ=37)                                                                     tataGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTAc                                                                                                                                                                           >  1:28437/1‑251 (MQ=255)                                                                      ataGTGTTGCTCAGCGCTGTCAGCGCGGAAGAGCAGGATGAGGACGCCTCGTAGGCCGACTCCAAAGTGCCTGGTGATTGGATGTGGTAAACATATATAAGGCGTGCGCAGAAGCCTGAGGCGACG‑GTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATACGACGTTTACGCCGCATCCGGCATTTGTGGGCTGATGCCTGATGCGACGCTTT‑CGCGTCTTATCATGCCTAca                                                                                                                                                                          <  2:292570/251‑2 (MQ=255)                                                                           gtTGTACAGCGCTCTCAGCGCCTAAGAGCAGGATGAATACGCCCCGTATGCCGACTTCAAAGTGCCTGATTATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCTTTATGCGACG‑CTATGCGCGTCTTATCAGGCTTACGTTTTATGGGCTAAGTGTAGACCGGTTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCtt                                                                                                                                                                      <  1:20810/250‑1 (MQ=38)                                                                           gtTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGCTAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTAcaaatct                                                                                                                                                                     >  1:332211/1‑245 (MQ=25)                                                                               cacaGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTAcaaatctgtac                                                                                                                                                                 <  1:425220/251‑11 (MQ=25)                                                                                               ccGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTGCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGa                                                                                                                                                 <  1:254469/251‑1 (MQ=255)                                                                                                       aaGGATGAAGACGACCCTTATGCCGACTACCATTGGCCTGAGGATTGGAGGTGTTAAACATAGCGAAGTGCTGCTCTTATGCCTGA‑GACACGGCTATGCGAGTCTTATCAGACCGACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCTCGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGt                                                                                                                                         <  2:353225/250‑1 (MQ=255)                                                                                                             gaagtcgacgagTATGCCGTCATAAAAGTGCCGGATGATTTGATGTGGTAAACATAGATAAGTAGTGCGCAGACGCACGATGCGACT‑CTATGAGCGTTTTATCAGGCCTATGTTTTATGGGCGAAGTGTAGACCGGATAAGGGGTTCGCGTCGCATTCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                    <  2:179052/239‑1 (MQ=255)                                                                                                             gAAGACGACCCTTATGCCGACTTCAAAGTGCCTGATGATTTGATGGGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                   <  1:255723/251‑1 (MQ=255)                                                                                                             gAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCACAGGTGCCTGCTGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCATGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                   >  1:316222/1‑251 (MQ=255)                                                                                                                         tATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                         >  1:390372/1‑249 (MQ=255)                                                                                                                                  ttCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCtat                                                                                                              >  1:35961/1‑251 (MQ=255)                                                                                                                                       ggTGCCTGATGATTGGATGTGGTAAACCTAGATATGTCGTGCGACGATGCGTGATGCGACG‑CTATGCGCGTCTTACAAGGCCTACGTTTTAGGGGCGAAGTGTAGCCCGGATAAGGCGGTCGCGCCGCATCCTGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGTCAGTCGTGCTATTATCa                                                                                                         <  2:503918/250‑1 (MQ=37)                                                                                                                                        gTGCCTGATGATTTGATGTGGTAAACATCGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                                                                                                                                                                                                     >  2:86093/1‑143 (MQ=38)                                                                                                                                        gTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                                                                                                                                                                                                     <  1:86093/143‑1 (MQ=38)                                                                                                                                         ggCCTGAGGATTTGATGTGGTAAACATAGATAATTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTAGCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGTCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAAc                                                                                                       <  1:486635/250‑1 (MQ=255)                                                                                                                                           ccTGATGAGTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGc                                                                                                     <  1:21980/251‑1 (MQ=255)                                                                                                                                                     tGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGtt                                                                                           >  2:167201/1‑251 (MQ=255)                                                                                                                                                             tAAACATAGATAAGTCGTGCGCAGATGCCTGATGAGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCCTATTCAGTTTATTGGCg                                                                                   >  1:348821/1‑251 (MQ=255)                                                                                                                                                                 cATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCCGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCCCACCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGCCGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTCTTGGCGTGGt                                                                               >  1:305415/1‑251 (MQ=255)                                                                                                                                                                        aaGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAat                                                                        >  1:483144/1‑251 (MQ=255)                                                                                                                                                                          gTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGCTGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAatat                                                                      >  1:337817/1‑251 (MQ=255)                                                                                                                                                                              tGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTc                                                                  >  1:343517/1‑251 (MQ=255)                                                                                                                                                                               gggcAGATGCCTGATGCGGCG‑CTATGCGCGTCTGATCAGGCCTACGGTTTATGGGCGAAGGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTc                                                                  <  2:28437/248‑1 (MQ=255)                                                                                                                                                                                     aTGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGGGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAAAATGCTCGcccat                                                            >  2:445008/1‑247 (MQ=255)                                                                                                                                                                                         cTGATGAGACG‑CTATGCGCGTCTTATCATGCCTACGGTTTATGGGAGAAGAGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATTCCTGATGCGACCATTTGCGCGCCTTGTATGGCGAACGGCTTCCGGGTGACGTGTCTGACGGATAAGGAGATTACGCCGCCACCTGCGGTCGTGCTGTTATATATGCATATCCAGTTTATCGGCCTggatggtactaaactcgcaaataaaac                                                       >  2:407912/1‑226 (MQ=9)                                                                                                                                                                                           gaggcggcG‑CTATGCCGGTCTTATGAGGCCTACGGTTTTTGGGCGAAGTGTAGACCTGATAGGGCGTTCACGCCGCATCCTGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAATTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAc                                                       <  2:85788/242‑1 (MQ=37)                                                                                                                                                                                                           gcgtttgaTCAGGCCTACGGTTTATGGGCGAAGGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAg                                      <  2:391011/244‑1 (MQ=255)                                                                                                                                                                                                           gcgTCTTATCAGGCCTACGGTTTATGGGCGAAGGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATGATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAg                                      <  1:305035/251‑1 (MQ=255)                                                                                                                                                                                                            ggtctattagggCCTACGGTTAATGGGCGAAGTGTGGACGGGATAAGGCGTTCACGCCGCATCAGGAAGTCGTGCGC‑CATGCCTGAGGCGACGTTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTGTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAg                                      <  2:451427/240‑1 (MQ=37)                                                                                                                                                                                                                ttATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCtgtg                                 >  1:352432/1‑251 (MQ=255)                                                                                                                                                                                                                  atcgggcttACGGTTTAGGGGCGAAGTGTAGACCGTATATGGCGGTCACGCCGCAGCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCATACGGCTTACGGGCGAAGTGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGAt                               <  2:282364/243‑1 (MQ=38)                                                                                                                                                                                                                  aTCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGGTGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGGATTGCTATTATCAACGCATATTCAGTTTATTTGCGTGGTAGGCAATATTCTCGCACATAACACTGGAGACACATTAATATCtgtact                               >  2:279409/1‑248 (MQ=37)                                                                                                                                                                                                                  aTCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGAt                               <  1:280591/251‑1 (MQ=255)                                                                                                                                                                                                                    atgcccTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCCCGACGAATCCGGCAGTCGGGCGC‑CATGCCTGATGCGACTCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGTCGAAGTTTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATc                              <  2:149050/246‑1 (MQ=38)                                                                                                                                                                                                                          tACGGATTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGCTGCGCCGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                       >  1:301236/1‑251 (MQ=255)                                                                                                                                                                                                                           aCGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                       <  2:158832/250‑1 (MQ=255)                                                                                                                                                                                                                               tttATGGGCGAAGTGTAGACCGGCTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAAGCAATATACTTGCCAATAAAACTGGAGAAAAAATAATAGGTGTGCTCATAACCCTCGTg                  >  2:354826/1‑251 (MQ=38)                                                                                                                                                                                                                                 taggggCGAAGTGTATACCGGAAAAGGCGTTCACGCCGCGTCCGTCAGTCGTGCGC‑CATTCCTGAGTCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTAGTATCAACGCATATTCAGTTTATTGGCTTGGTAGGCAATAGGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGa                 <  2:323210/247‑1 (MQ=38)                                                                                                                                                                                                                                     ggCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATTATAACCCTCGTGATAACg            >  1:46975/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGAtt        <  1:480523/251‑1 (MQ=255)                                                                                                                                                                                                                                               aGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTc  >  1:203508/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                AGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTC  >  NC_000913/377153‑377640

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