Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,290,669	G→A	W143* (TGG→TAG)	rssB →	regulator of RpoS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,290,669	0	G	A	98.0%	170.7 / ‑6.6	50	W143* (TGG→TAG)	rssB	regulator of RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (42/7);  total (42/8)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTGGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATAATGGCGTACTTGCTGCCTTGTTATTACGCGCATTGT  >  NC_000913/1290426‑1290908                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 cACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGc                                                                                                                                                                                                                                            >  1:880443/1‑250 (MQ=255)      aTGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGt                                                                                                                                                                                                                                      <  2:695625/251‑1 (MQ=255)        gAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTg                                                                                                                                                                                                                                    >  1:29121/1‑251 (MQ=255)             ggCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAAc                                                                                                                                                                                                                               >  2:714754/1‑251 (MQ=255)                  aaaCTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGGGACTAGGATGCAATGGTTGATAGCCCTGc                                                                                                                                                                                                                          <  1:419703/251‑1 (MQ=255)                           gAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAgcggc                                                                                                                                                                                                                <  1:1055817/251‑1 (MQ=255)                           gAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAgcggc                                                                                                                                                                                                                <  2:721581/251‑1 (MQ=255)                           gAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAgcggc                                                                                                                                                                                                                >  1:1210656/1‑251 (MQ=255)                             gCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAgcggcg                                                                                                                                                                                                               >  2:494310/1‑250 (MQ=255)                                 ataCGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAAGGGTTGATAACCCTGCC‑GCAGCGGCGAAat                                                                                                                                                                                                           <  1:970075/250‑1 (MQ=255)                                    cGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAattatt                                                                                                                                                                                                       >  2:183176/1‑251 (MQ=255)                                       aaaagagGCGACCAGACCCCAGTTCTGGGGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGAGAGATGGTTTTTGCCTGTTTCTATCCCACCATGTTTAATTCGTGCGGTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGAGGCAATGGTTGATTACCCTACC‑GCAGCGGCGAACTTATTAc                                                                                                                                                                                                     >  2:1168567/4‑250 (MQ=255)                                          agagGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAgg                                                                                                                                                                                                  >  2:488360/1‑250 (MQ=255)                                                 aCCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACa                                                                                                                                                                                           >  1:395967/1‑250 (MQ=255)                                                  ccAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAAc                                                                                                                                                                                         >  2:1020102/1‑251 (MQ=255)                                                   cAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAAcc                                                                                                                                                                                        >  1:1034455/1‑251 (MQ=255)                                                          cAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGtt                                                                                                                                                                                 >  1:1110563/1‑251 (MQ=255)                                                                tGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTg                                                                                                                                                                                                                                    <  1:1189720/195‑1 (MQ=255)                                                                tGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCATGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTg                                                                                                                                                                                                                                    >  2:1189720/1‑195 (MQ=255)                                                                   tGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTg                                                                                                                                                                        >  2:436940/1‑251 (MQ=255)                                                                      tatCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGAtt                                                                                                                                                                     >  2:400816/1‑251 (MQ=255)                                                                              aCTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCCGGTGATTTCCCATTg                                                                                                                                                             >  2:1238746/1‑251 (MQ=255)                                                                                tGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGcc                                                                                                                                                           >  1:598923/1‑251 (MQ=255)                                                                                             gATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAgg                                                                                                                                                                          >  1:789282/1‑223 (MQ=255)                                                                                              atatTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTg                                                                                                                                                             >  2:575/1‑235 (MQ=255)                                                                                                     cAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCACTTgg                                                                                                                                      >  1:1164123/1‑251 (MQ=255)                                                                                                      aaaGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGCTAATTATCGTCAATTGGt                                                                                                                                     >  1:562747/1‑251 (MQ=255)                                                                                                                 tCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCCTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGAc                                                                                                                           >  1:266192/1‑250 (MQ=255)                                                                                                                           gAAGATGTTTTGCTGACACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCt                                                                                                                >  2:756497/1‑251 (MQ=255)                                                                                                                            aaGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTg                                                                                                               >  1:239988/1‑251 (MQ=255)                                                                                                                                     tgtagaaccaagtaaaaaagcTGAATCGCTTGCGGGAGTTGTTTTTTCCCTGTCCCTACCCCAGTTTTTTTAATTCGCGGGTGGAGGAAGAG‑AAAGG‑TTTTTTTCGGAACTGGAATGCAAGGGTTGA‑AACCCCGCCCGCGGGGGCAAAATTATTACAGAAACAAAAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGa                                                                                                       <  1:494310/230‑1 (MQ=255)                                                                                                                                        tGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATAtt                                                                                                   >  2:933338/1‑251 (MQ=255)                                                                                                                                           aaCCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGcc                                                                                                >  2:421802/1‑251 (MQ=255)                                                                                                                                               aGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCAc                                                                                            >  1:490145/1‑251 (MQ=255)                                                                                                                                                  tAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACttt                                                                                         >  2:327847/1‑251 (MQ=255)                                                                                                                                                        tCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGaa                                                                                    >  2:253892/1‑250 (MQ=255)                                                                                                                                                        tCTGAATCGCTTGCGCGAGATGGTTTTTACCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGaa                                                                                    >  2:759783/1‑250 (MQ=255)                                                                                                                                                                             ggTTTTTGCCTGTCTCGATCCCAGCATGTTTAATTCGCGCGGTGCGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCACTTGGTTGCCGCGGACCAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGTAATACGATATGGCATTTTGTTGcc                                                              >  2:266417/1‑251 (MQ=255)                                                                                                                                                                                 tttGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTg                                                           >  2:1001105/1‑250 (MQ=255)                                                                                                                                                                                  ttGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGAt                                                         >  1:162872/1‑251 (MQ=255)                                                                                                                                                                                     ccTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTc                                                      >  2:743168/1‑251 (MQ=255)                                                                                                                                                                                        ttctcTACCCCAGCATTTTTAATTCGCGCGGTG‑GGAAAGGGAAAGGGTTTTTT‑CGCGAGTAGGATGAAATGTTTAACACCCCTGCC‑GCAGGGGCAAATTTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCAc                                                    <  1:253892/249‑1 (MQ=255)                                                                                                                                                                                                     cATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTg                                                                                                  >  2:1114668/1‑191 (MQ=255)                                                                                                                                                                                                     cATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGCGCTGGACAt                                       >  2:23598/1‑250 (MQ=255)                                                                                                                                                                                                             aTTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGACGGTTCAGCAGGGGCTTTACCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCTCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCCCCCTAGCTGGTAATAATGGCGt                               >  2:716097/1‑250 (MQ=255)                                                                                                                                                                                                                 gcgcgTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTa                                                                   >  1:1243903/1‑210 (MQ=255)                                                                                                                                                                                                                    cgTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGt                                                                                                             >  1:902175/1‑165 (MQ=255)                                                                                                                                                                                                                           gaagaggaaAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATAATGGCGTACTTGCTGCCTTGtt                >  2:880403/1‑251 (MQ=255)                                                                                                                                                                                                                             agaggaaAGG‑CTTTTT‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATAATGGCGTACTTGCTGCCTTGttat              >  2:44656/1‑251 (MQ=255)                                                                                                                                                                                                                                          tttt‑CGCGACTAGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATAATGGCGTACTTGCTGCCTTGTTATTACGCGCATTGt  >  2:783272/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 CACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGG‑CTTTTT‑CGCGACTGGGATGCAATGGTTGATAACCCTGCC‑GCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGATGTCACCCGAGCTGGACATAATGGCGTACTTGCTGCCTTGTTATTACGCGCATTGT  >  NC_000913/1290426‑1290908

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