Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	109,394	G→T	L372L (CTG→CTT)	secA →	protein translocation ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	109,394	0	G	T	97.1%	222.3 / NA	69	L372L (CTG→CTT)	secA	protein translocation ATPase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (26/41);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTGGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGT  >  NC_000913/109148‑109639                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      catcatGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGCAGCGAAAGAAGGGGTGCAGATCCAGAACGAAAACCAAACgcttgctt                                                                                                                                                                                                                                                   >  2:435498/1‑251 (MQ=255)     atGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGAt                                                                                                                                                                                                                                               <  1:562315/251‑1 (MQ=255)     atGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGAt                                                                                                                                                                                                                                               >  1:227335/1‑251 (MQ=255)       ggATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCCTGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGGAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATc                                                                                                                                                                                                                                              >  1:453437/1‑250 (MQ=255)                  aGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAAAGGCGGCGCTGCGCGCTCCTGCGCTGTTTTCCCGTGACCTCGACTACATCGTTTAAGATGGGGAAGTTATCCTCGTTTACCAACACACCGGTCGTACCATGAAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCCGATCCAAAAAGAAAACCAAACGATGGCTTCGATCACATTCCCGAac                                                                                                                                                                                                                                  >  1:504109/1‑251 (MQ=255)                  aGTCTCTGTACTCTCCGGACAAAATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGGGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAac                                                                                                                                                                                                                                  >  1:161911/1‑251 (MQ=255)                    tctcTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTAACCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGGGCGCAGGTCCCGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAact                                                                                                                                                                                                                                 >  1:274934/1‑250 (MQ=255)                                ccGGCCAACATCATGCTGATGCACCACGTAACTGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTTAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGGCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTGGCTTCGATCATCTTCCAGAACTACTTCCGTCTGTa                                                                                                                                                                                                                    <  2:330622/251‑1 (MQ=255)                                ccGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCt                                                                                                                                                                                                                       <  1:360375/248‑1 (MQ=255)                                ccGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTg                                                                                                                                                                                                                      <  1:492913/249‑1 (MQ=255)                                ccGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGGGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCGTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTa                                                                                                                                                                                                                    <  2:187733/251‑1 (MQ=255)                                 cGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTAt                                                                                                                                                                                                                   <  1:401546/251‑1 (MQ=255)                                   gCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGGCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCCAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATg                                                                                                                                                                                                                  >  1:177270/1‑250 (MQ=255)                                   gCCAACAGCATGATGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGGGACGTCGACTACATCGTTAAAGATGGTGAAGTGATCATCGTTGACGAACACACCGGGCGTACCATGCAGGTCCGTCGCTGGTCCGATGGTCGGCACCAGGCGGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCAACTTCCAGAACTACTTCCGTCTGTATGa                                                                                                                                                                                                                 <  2:528319/251‑1 (MQ=255)                                      aaCATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGAGCTGTTTACCCCTGCCGTCGACTACATCGTTAAAGATGGGGAAGTTGTCATCGTTGACCAACACACCGGGCGGACCCTGCCAGGCCGTCGCTGGGCCGGTAGGCTGCCCCCGGCTGTTGGAGCGGAAACAGGTGTGCAGATCCCGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAAATACTTCCGTCTGTATGaca                                                                                                                                                                                                               >  2:164393/1‑248 (MQ=255)                                                 gATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGGCGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTcc                                                                                                                                                                                                                           >  1:245709/1‑227 (MQ=255)                                                 gATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTcc                                                                                                                                                                                                                           <  2:245709/227‑1 (MQ=255)                                                    gCACCACGTAACCGCGGCGCTGGGCGCTCATGCGCTGTTTACCCGGTACGTCGACTACAGCGTTAAAGATGGTGAAGTTATCATCGTTGCCGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGTTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGGCTGTATGAAAAACTGGCGGGGATGa                                                                                                                                                                                                <  2:555737/251‑1 (MQ=255)                                                      cccacGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGAcc                                                                                                                                                                                              <  1:201479/250‑1 (MQ=255)                                                         acGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGt                                                                                                                                                                                           <  1:168972/251‑1 (MQ=255)                                                             aaCGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCACATCCAGAACGAAAACCAAACGCTTGCTGCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTActg                                                                                                                                                                                       <  2:215128/251‑1 (MQ=255)                                                             aaCGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGGTGACGAACACACCGGGCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACTAAAACCAAACGCTTGCTTCGATCACCTTCCATATCTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTAct                                                                                                                                                                                        <  1:274666/250‑1 (MQ=255)                                                              aCGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTActgc                                                                                                                                                                                      <  2:123889/251‑1 (MQ=255)                                                                ggcggcGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAACCTGGCGGGGATGACCGGTActgctg                                                                                                                                                                                    >  1:365432/1‑251 (MQ=255)                                                                   ggcgCTGCGCGATCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGAt                                                                                                                                                                                  >  1:124585/1‑250 (MQ=255)                                                                    gcgcTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTTCATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATAc                                                                                                                                                                                <  1:69924/251‑1 (MQ=255)                                                                     cgcTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCCTCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTCTGAAAAACTGGCGGGGGTGACCGGTACTGCTGCTAcc                                                                                                                                                                               >  2:511235/1‑251 (MQ=255)                                                                     cgcTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATAcc                                                                                                                                                                               <  1:253078/251‑1 (MQ=255)                                                                       cTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGCAAGAAGGTGAGAAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGAGGACCGGTACTGCTGATACCgc                                                                                                                                                                             >  1:263265/1‑250 (MQ=255)                                                                             gcTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACAAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAACACTGGCGGGGATGCCCGGTAATGCTGAGAACCAAGCttt                                                                                                                                                                       >  2:23962/1‑251 (MQ=255)                                                                                       tGTTTACCCGTGACGGCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAAAACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCAAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGGTTACCGGTACTGGTGATACCGAAGCCTTCGAATGTAGc                                                                                                                                                             >  2:43234/1‑251 (MQ=255)                                                                                                   tcggcGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCGTGCAGGGCCGTCGCTGGTCCGATTGTCTGCACCAGGCTGTGTAAGCGAAAGAATGTGTGCATATCCAGAACGAAAACCAAACGCTTGCGTCGATCATCTTCCATAATTACTTCCGTCTGTATGAAAAACTGGCGGTGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAAGCTACaa                                                                                                                                                  <  2:520053/246‑1 (MQ=255)                                                                                                    cGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCTGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAg                                                                                                                                                 <  1:273621/250‑1 (MQ=255)                                                                                                            aCATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTTCCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTTCGTCTGTATGAAAAACTGGCGGTGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTGAGCTCAATCTACAAGCTGGATa                                                                                                                                          <  2:378318/249‑1 (MQ=255)                                                                                                                        aTGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGGACTGCTGATACCGAAGCTTTCGAATTTAGCTCAAACTACAAGCTGGATACCGTCGTTGTTCCg                                                                                                                            >  1:306902/1‑251 (MQ=255)                                                                                                                              aaGTTCTCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAGAACTGGCGGGGGTGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGGTCCGACCAAc                                                                                                                      >  1:400713/1‑251 (MQ=255)                                                                                                                                gTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCg                                                                                                                    <  2:419947/251‑1 (MQ=255)                                                                                                                                     cagcGTTGACCAACACACCGGACTTCCCCGGCAGGGACGTCGCTGGTCCGATGGGCGGACCCAGGCTGTGGACGCGAACGAAGGTGTGCAGATTCAGAACGATAACCATACGTTTGTGCCGACTACCTTCCAGAACCACTACTGTCGGTATTAAAAACTAGCGGGGATGGCCGGTACAGCTGATACCGAATCGTTCGAATTTAGCTCAATCTACAAGCTGGATACCGGCGTGGTTCCGACCAACCGTCCa                                                                                                                <  2:485394/247‑1 (MQ=255)                                                                                                                                           tGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTACGATGGTCTGCACCAGGCTGTGGAAGAGAAAGAAGGTGCGCAGCTCCAGAACGCACACCAAACGCTTGCTTCGATCAACTTCCAGCAATACTTACGTCTGTATGACAAACTGGCGGGGGTGCACGGTACGGCTCATCACGACGCTTTCGAATTTAGCGCAAACTCCAAGCTGGATACCGTCGTTGATCCGACCAAcccgccaactattt                                                                                                         >  2:195290/1‑239 (MQ=255)                                                                                                                                             aCGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGt                                                                                                       <  2:227335/251‑1 (MQ=255)                                                                                                                                              cGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGt                                                                                                       >  1:440520/1‑250 (MQ=255)                                                                                                                                                aaCACACCGGTCGCACCATGCAGGGCCATCGCTGGGCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCGTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAAGCTACTAGCTGGATACCGTCGTTGTCCCGACCAACCGTCCAATGATTCGTaaa                                                                                                    <  2:297342/251‑1 (MQ=255)                                                                                                                                                 acacacCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATGCGTAAAg                                                                                                   <  1:128307/251‑1 (MQ=255)                                                                                                                                                      ccGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTAACGCTGATACCGCAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGTCCAACCGGCCAATGATTCGTAAAGAGCTg                                                                                              >  2:118448/1‑251 (MQ=255)                                                                                                                                                               ccATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAGAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGAGACCGTCGTGGTTCCGACCACCCGTCCAATGATTCGTAAAGATCTGCCGTACCTg                                                                                     <  2:400713/251‑1 (MQ=255)                                                                                                                                                                 aTGCAGGGCCGTCGCTGGGCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGt                                                                                   <  1:511235/251‑1 (MQ=255)                                                                                                                                                                       ggCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTAGATACCGTCGTTGTTCCGACCAACCGTCCACTGATTCGTAAAGATCTGCCTGCCCTGGGCTACAt                                                                             >  1:75079/1‑251 (MQ=255)                                                                                                                                                                          cGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGGTACCGTCGTTTTTCCGCCAAACCGTCCAATGATTCGTAAAGAGCTGCCGGACCTGGTCTACATGAc                                                                          >  2:328319/1‑251 (MQ=255)                                                                                                                                                                                     cGATGGTCTGCCACAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGATTTCGAATTTAGCTCAAGCTACAAGCTGGATGCCGTCGTGGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGaaa                                                               <  2:504109/251‑1 (MQ=255)                                                                                                                                                                                      gATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGGGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGAGACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCTGACCTGGTCTACATGACTGAAGCGGaaaa                                                              <  2:474298/251‑1 (MQ=255)                                                                                                                                                                                         ggTCTGCCCCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATTACCGGGACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGaaaaaa                                                            <  2:365432/250‑1 (MQ=255)                                                                                                                                                                                            cTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATAACGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTTTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTc                                                         >  2:377711/1‑250 (MQ=255)                                                                                                                                                                                               cACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGc                                                     <  1:131630/251‑1 (MQ=255)                                                                                                                                                                                               cACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGc                                                     <  1:126426/251‑1 (MQ=255)                                                                                                                                                                                                            aaGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGCCTTGGTCTACATGACTGAAGCGGAAAAAAGTCAGGCGAGAATTGAAGat                                        <  1:164393/251‑1 (MQ=255)                                                                                                                                                                                                                  aagaagGTGTGCATATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGTACCTGGTGTACATGACTGAAGCGGAAAAAATTCAGGCGATCATGGAAGATATCaaa                                  <  2:474272/251‑1 (MQ=255)                                                                                                                                                                                                                  aagaagGTGTACAGATACAGAACGAAAACCAAACGCGTTCTTCGATCCCCTTCCAGAACTACCTCCGTCTGTATGAAAAACTGTCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCTAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAAACGTACAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCaaa                                  <  2:307304/251‑1 (MQ=255)                                                                                                                                                                                                                        gtgtGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGt                            <  2:430662/251‑1 (MQ=255)                                                                                                                                                                                                                           tgCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACt                         >  1:18492/1‑251 (MQ=255)                                                                                                                                                                                                                           ggCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTTCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGTATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCAGTGAAGATATCAAAGAACGTACt                         <  2:483902/250‑1 (MQ=255)                                                                                                                                                                                                                            gCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATAGCAAAGAACTTACTg                        <  1:459590/251‑1 (MQ=255)                                                                                                                                                                                                                             cAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACcgtc                                                                                                                                     <  2:289131/141‑1 (MQ=255)                                                                                                                                                                                                                             cAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACcgtc                                                                                                                                     >  1:289131/1‑141 (MQ=255)                                                                                                                                                                                                                              aGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCg                      <  1:435498/251‑1 (MQ=255)                                                                                                                                                                                                                                aTCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACa                                                                              >  1:303973/1‑193 (MQ=255)                                                                                                                                                                                                                                aTCCAGAACGAAAACCAAACGCTTGCTGCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACa                                                                              <  2:303973/193‑1 (MQ=255)                                                                                                                                                                                                                                     gAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCCACCGTCCAATAATTCGTAAAGATCTGCCGGACCCGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCCATGAAGATATAAAAGAACGTACTGCGAAAggac               >  1:403204/1‑249 (MQ=255)                                                                                                                                                                                                                                     aaaCGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGcc               <  1:284449/250‑1 (MQ=255)                                                                                                                                                                                                                                         gAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGcc           >  1:488355/1‑251 (MQ=255)                                                                                                                                                                                                                                                  aCGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTGTCGAATTTAGCTCAAGCTACAAGCTGGATACCGGCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAAGTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCtggt  <  2:453437/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CATCATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTGGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTTGTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCTGGT  >  NC_000913/109148‑109639

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