Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	377,407	0	C	T	22.8%	88.6 / 33.3	58	intergenic (+96/+128)	yaiL/frmB	DUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base T (5/8);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-02
Rejected as consensus: Frequency below/above cutoff threshold.

AAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTC  >  NC_000913/377157‑377640                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           aaGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTc                                                                                                                                                                                                                                            >  1:390703/1‑251 (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)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           AAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACG‑CTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTC  >  NC_000913/377157‑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: ‑