Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	24,091	C→T	R567R (CGC→CGT)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	24,091	0	C	T	96.2%	164.9 / NA	53	R567R (CGC→CGT)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/33);  minor base G (1/1);  total (19/34)
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.94e-01

TACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAATACCGTTTCGCCGCAGGATGTGATGAACAAACTGGGCGCGGATATTCTGCGTCTGTGGGTGGCATCAACCGACTACACCGGTGAAATGGCCGTTTCTGACGAGATCCTGAAACGTGCTGCC  >  NC_000913/23854‑24337                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       tACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCCGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATgg                                                                                                                                                                                                                                           >  2:11691/1‑250 (MQ=255)  aCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGt                                                                                                                                                                                                                                           <  1:371679/250‑1 (MQ=255)  aCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTc                                                                                                                                                                                                                                          >  1:7771/1‑251 (MQ=255)   cccTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTc                                                                                                                                                                                                                                          <  1:502111/250‑1 (MQ=255)   cccTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCt                                                                                                                                                                                                                                         >  1:366720/1‑251 (MQ=255)     cgcgAACTGATGGAAGAAGGGGCCAAACGCGTTGCAGGTTACGGCATCCAGGCGTGGTGGGATCTATATGCGAAAGAGATCCTCGGCGAAGAAGCTGACAAGTATGTGAAATTGCCGGATACTCTC‑GACGTATGGTTTGCCTCCGGATCTACCCACTCTGCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCTGCACCGTGGTTGGTTCATGTCtgc                                                                                                                                                                                                                                       <  2:373539/248‑3 (MQ=255)       cgAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTccc                                                                                                                                                                                                                                     >  1:8109/2‑251 (MQ=255)                   agaagTGGCAAAACGCGTTGCAGTTAACGGCATCCAGGCGTAGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGGGCCGGATACTCTC‑GACGTATGGATTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACAGTGGTTGGTTGATGTATTCCCTAATGATCTccc                                                                                                                                                                                                                         >  2:333189/1‑250 (MQ=255)                    gaagTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCAc                                                                                                                                                                                                                        <  1:73381/251‑1 (MQ=255)                             aaaCGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGaa                                                                                                                                                                                                               <  1:170297/251‑1 (MQ=255)                                cgcggtgaaattaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAATATATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCGTCTTTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTGCATGTCTTCCCTAATGATCTCCACCGCGATGaaagg                                                                                                                                                                                                            <  2:420144/236‑4 (MQ=255)                                   gttgaagttaacggCATCCAGGCGGGGTGGGATCTCGATGAGAAAGCGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGGATCTGGAAGGTTCTGACCCGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGCAGAAgggca                                                                                                                                                                                                          >  2:405826/13‑248 (MQ=255)                                         gttaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGcc                                                                                                                                                                                                   <  1:375991/245‑1 (MQ=255)                                            aacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGt                                                                                                                                                                                                 >  1:326704/4‑250 (MQ=255)                                               ggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCg                                                                                                                                                                                             <  1:426267/251‑1 (MQ=255)                                                gCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGt                                                                                                                                                                                            <  1:286426/251‑1 (MQ=255)                                                         cGTGGTGGGATCGCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTGCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTg                                                                                                                                                                                   <  1:293340/251‑1 (MQ=255)                                                               gggagcccGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCAc                                                                                                                                                                             <  2:158392/244‑1 (MQ=255)                                                                 gATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGGGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGGCTTCCCTAATGATCTCAACCGCGATGAAAGGAAAAGCGCCGTATCGTCAGGTACTTACCAACgg                                                                                                                                                                           >  2:368565/1‑251 (MQ=255)                                                                       gATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCtttcc                                                                                                                                                                     >  2:51790/1‑249 (MQ=255)                                                                             aaaGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGa                                                                                                                                                               <  1:370176/251‑1 (MQ=255)                                                                              aaGAGATCCTCGGCGACGAAGCTGCCGAGTATGTGAAAGTGCCGGATTCTCTC‑GACGTATGGTTTGACTCCGGACCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATTTATCTGGAAGGTTCTGACCAGCGCCGTGGTTGGTTCATGGCTTTCTTGATGATCGCCTCCGTGATGAAAGGCAAGGCGCCGTAGCGTCAGGATCTGATCCAGGGCTCTACCGTGGAt                                                                                                                                                              <  2:211957/251‑1 (MQ=255)                                                                               agagATCCTCGGCGACGAAGCTGACCAGTAAGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCTTCAGGTACTGATCCACGGCTTTACCGTGGATg                                                                                                                                                             <  2:289218/251‑1 (MQ=255)                                                                                 agaTCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGt                                                                                                                                                           <  1:481228/251‑1 (MQ=255)                                                                                            gacgaAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGAGGATGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGAACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCa                                                                                                                                                 >  2:384546/1‑250 (MQ=255)                                                                                             acgaAGCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCa                                                                                                                                                 >  1:259883/1‑249 (MQ=255)                                                                                                aaGCTGACCAGTATGTGAAAGTGCCGGATAC‑TCGCGACGGATGGTTTGACTCCGGATCTACCGACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGGCTTCCCTAATGATCTCCACCGCTATGAAAGGCAAAGCGCCGTATCGGCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATg                                                                                                                                            <  2:56897/251‑1 (MQ=255)                                                                                                  gCTGACCAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTc                                                                                                                                          <  1:427841/251‑1 (MQ=255)                                                                                                     gaccAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAGTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTTACCCACGGTTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTaa                                                                                                                                       <  1:191976/248‑1 (MQ=255)                                                                                                        cAGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACGCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTTTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCTGCAAGATGTCTAAATc                                                                                                                                    <  1:409289/251‑1 (MQ=255)                                                                                                         aGTATGTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATcc                                                                                                                                   <  2:358757/251‑1 (MQ=255)                                                                                                              gTGAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACAGGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAGCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGTCCCACGGCTTTACCGTGGATGGTTAGGGCCGCAAGATGTCTAAATCCATCgg                                                                                                                              <  1:444009/251‑1 (MQ=255)                                                                                                                gAAAGTGCCGGATACTCTC‑GACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGcc                                                                                                                            >  2:253612/1‑250 (MQ=255)                                                                                                                                        tATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGAAATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACAACGTTTCGCCGCAGGAtgtg                                                                                                     >  1:55187/1‑251 (MQ=255)                                                                                                                                         ttGGTTTGACTCCGGATCTACCCCCTCTTCTGTTGTTGACTTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAAGGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGAtgtg                                                                                                     <  2:81624/249‑1 (MQ=255)                                                                                                                                         aTGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGtga                                                                                                    <  2:85138/251‑1 (MQ=255)                                                                                                                                              tgtaCTCCGTATCGACCCACTCTTCTGTTGTTGATGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGGTCAGACCAGCACCGTGGTTTGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACTGCTTGACCGTGGATGTGCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGaa                                                                                                <  1:30227/247‑1 (MQ=255)                                                                                                                                                gACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACa                                                                                              <  1:420372/250‑1 (MQ=255)                                                                                                                                                 aCTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACaaa                                                                                            <  1:242477/251‑1 (MQ=255)                                                                                                                                                   tCCGGATCGTCCCACTCTTCTGTTGTTGACGTGCGTCCGGAAGTTGCCGGTCACGCAGCGGACATGTATCTGGAATGTTCTGACCAGCACCGGGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCTCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGTATTGTCAGGGCCGCAAGATGTCTAAAGCCATCGGCAACACCGTTTCGCCGTCGGATGTGATGATCAAAc                                                                                           <  2:324802/250‑1 (MQ=255)                                                                                                                                                    ccGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAGTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCCGCACCGTGGTTGGTTCATGGCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACAAACTg                                                                                         <  2:391510/251‑1 (MQ=255)                                                                                                                                                                aCTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACAAACTGGGCGCGGATAtt                                                                             <  1:415489/251‑1 (MQ=255)                                                                                                                                                                    ttctGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCCGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGGTGTGATGAACaca                                                                                            >  2:83466/1‑230 (MQ=255)                                                                                                                                                                    ttctGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACaaa                                                                                            <  1:83466/232‑1 (MQ=255)                                                                                                                                                                             tgACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCATCTCCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCATGGCCGCAAGATTTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACAAACTGGTCGCGGATATTCTGCGTCTgtgg                                                                 <  1:453039/250‑1 (MQ=255)                                                                                                                                                                                      tCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGGGGATGGTTCATGTCTTACCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCATTACCGTGGATGGACAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAAAAAACTGGGCGAGGATAATCTGCGTCTGGGGGGGGAATCa                                                        >  2:140976/1‑250 (MQ=255)                                                                                                                                                                                                 ccGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCCCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACAAACTGGGCGCGGATATTCTGCGTCTGTGGGTGGCATCAACCGACTACAcc                                            <  2:326704/251‑1 (MQ=255)                                                                                                                                                                                                          cAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACAAACTGGGCGCGGATATTCTGCGTCTGTGGGTGGCATCAACCGACTACACCGGTGAAATg                                   >  1:27385/1‑251 (MQ=255)                                                                                                                                                                                                           aGCGGACATGTAGCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACAAACTGGGCGCGGATATTCTGCGTCTGTGGGTGGCATCAACCGACTACACCGGTGAAATgg                                  <  2:119022/251‑1 (MQ=255)                                                                                                                                                                                                              ggACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCCCGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAACTCCATCGGCACCACCGTTTCCCCGCAGGATGTGATGAACAAACTGGGCGCGGCTATTCTGCGTTTGTGGGTGGCATCAACAGACTACACCGGTGAAATGGCCg                               >  2:238099/1‑251 (MQ=255)                                                                                                                                                                                                                             aGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACAAACTGGGCGCGGATATTCTGCGTCTGTGGGTGGCCTCAACCGACTACACCGGTGAAATGGCCGTTTCTGACTAGATcc                >  1:259952/1‑251 (MQ=255)                                                                                                                                                                                                                                    gACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACAAACTGGGCGCGGATATTCTGCGTCTGTGGGTGGCATCAACCGACTACACCGGTGAAATGGCCGTTTCTGACGAGATCCTGAAACg         <  1:131690/251‑1 (MQ=255)                                                                                                                                                                                                                                     aCCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTAACGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAACACCGTTTCGCCGCAGGATGTGATGAACAAACTGGGCGCGGATATTCTGCGTCTGGGGGTGGCATCAACCGACTACACCGGTGAAATGGCCGTTTCTGACGAGATCCTGAAACGt        >  1:279165/1‑251 (MQ=255)                                                                                                                                                                                                                                           aacgcgGTTGGTTCATGTCTTCCCTAATGATCTCCACAGCGATGTAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGAATGTCAGGGCCGCAAGCTGTCTAAATCCATCGGAAAAACCGTTTCGCCGCAGGTAGTGCTGAACAAACTGGGCGCGGATATTCTGCGTCTGTGGGTGGCATAAACCGACTACACCGGTGACATGACCGTTTCTggagcgatcttgacacaggtggcc  >  2:377586/3‑228 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       TACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTG‑GATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCAATACCGTTTCGCCGCAGGATGTGATGAACAAACTGGGCGCGGATATTCTGCGTCTGTGGGTGGCATCAACCGACTACACCGGTGAAATGGCCGTTTCTGACGAGATCCTGAAACGTGCTGCC  >  NC_000913/23854‑24337

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