Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,244	C→T	L148L (CTC→CTT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,244	0	C	T	97.2%	114.6 / NA	37	L148L (CTC→CTT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/14);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAA  >  NC_000913/3006‑3480                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             gggAGCGTTTTGGCCAGGAGCTAGGCAAGCGAATTCCAGTGGTGATGACTCTGGAAAAGAATATGCCAATCGGTGCGGGCTTAGGCTCCAGCGCCTGGCCGGTGGTCGCGGCGCTGATGGCGATGAATGAACGCTGTGGCAAGCCGCGTAATGACACTCGTTTGCTGGCGTTGATGGTCGAGCTGGAAGGACGAAGCTCCGGCAGTATTCATTGCGACATCGCGGGACCTTGTTTGCAGGGTGGTATGc                                                                                                                                                                                                                                    <  2:244478/249‑1 (MQ=255)                aggagcttggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGCAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCgaaga                                                                                                                                                                                                                    >  1:143743/12‑250 (MQ=255)                  gagcttggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGaaa                                                                                                                                                                                                                  >  1:28358/10‑250 (MQ=255)                   agcttggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAAGATGCCAATCGGGTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGGGGCACCGTGTTTTCTTGGTGGTATGTAGTTGATGATCGAAGAAAAc                                                                                                                                                                                                                <  1:104750/243‑1 (MQ=255)                              aaaaTTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGGCATCATcag                                                                                                                                                                                                     <  1:211390/250‑1 (MQ=255)                                           cGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCa                                                                                                                                                                                        <  1:32890/251‑1 (MQ=255)                                            gATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGCTGATCGAAGAAAACGACATCATCAGCCAGCAAGTGcc                                                                                                                                                                                         >  1:22579/1‑249 (MQ=255)                                                 ctcTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTg                                                                                                                                                                                           >  2:115267/3‑242 (MQ=255)                                                 ctcTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTg                                                                                                                                                                                           <  1:115267/240‑1 (MQ=255)                                                 ctcTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGttt                                                                                                                                                                                  >  2:65777/3‑251 (MQ=255)                                                       aaaaGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTtgatg                                                                                                                                                                              >  1:234938/1‑249 (MQ=255)                                                          aGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACCCTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCCGCAAGTGCCAGGGTTTGATGAGTgg                                                                                                                                                                         >  1:186588/1‑251 (MQ=255)                                                          aGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTgg                                                                                                                                                                         >  1:94956/1‑251 (MQ=255)                                                                 gccaaTCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTgg                                                                                                                                                                   <  1:190412/246‑1 (MQ=255)                                                                              ggCTTAGGCTCCAGCGCCTGATCGGTGGTCGCGGCGCTGTTGGCGATGAATGAACACTGTGGCAAGCCGATTAATGACACTCGTTTGCTGGCTTTGATGGGCGCGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGATTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGGGCTGGCGTATcc                                                                                                                                                     >  2:182071/1‑251 (MQ=255)                                                                              ggCTGAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCTCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATcc                                                                                                                                                     <  2:88449/251‑1 (MQ=255)                                                                                       tCCAGCGCCTGTTCGGTGGTCGCGGCGCCGATGGCGATGAATGACCACTGTGGCACGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACACCGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTa                                                                                                                                             >  2:13705/1‑250 (MQ=255)                                                                                                            gcggcgCTGATGGCGATGCATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGcc                                                                                                                         >  1:219285/1‑249 (MQ=255)                                                                                                                     aTGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTAtttt                                                                                                              >  1:168264/1‑251 (MQ=255)                                                                                                                                tgaaCACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGt                                                                                                   >  2:209602/1‑251 (MQ=255)                                                                                                                                    cacTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATcg                                                                                               >  1:67934/1‑251 (MQ=255)                                                                                                                                          ggCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCAAGGGCTATTTTACCGGCGCAGTATCGCCGCCa                                                                                         <  2:234938/251‑1 (MQ=255)                                                                                                                                               gCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATGTGACCGGCGCAGTATCGCCGCCATGatt                                                                                    <  2:235958/251‑1 (MQ=255)                                                                                                                                                 cGCTTAATGACTCTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGattgc                                                                                  <  2:215568/251‑1 (MQ=255)                                                                                                                                                         gACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCa                                                                          <  1:20841/251‑1 (MQ=255)                                                                                                                                                               cGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATATCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAACACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGCCGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCAAGGGCg                                                                    >  2:130786/1‑251 (MQ=255)                                                                                                                                                                    gCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATc                                                               <  2:168264/251‑1 (MQ=255)                                                                                                                                                                      tGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCCGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTCTCGCCGCCAGGGTTGCATTGCGGCCGGGCGACATCTg                                                             >  2:231210/1‑251 (MQ=255)                                                                                                                                                                           ttGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGctgg                                                        >  1:33582/1‑248 (MQ=255)                                                                                                                                                                           ttGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGct                                                          <  1:65777/249‑2 (MQ=255)                                                                                                                                                                               tGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGATTTCTTGTTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCGGGGCTATTTTACCGGAGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGTCATCTGGCTGGCttc                                                    >  2:188885/1‑251 (MQ=255)                                                                                                                                                                                   cgagctggaacttcgactcaccGGCAGCATTCCTTACGACAACGTGGTACCGTGTTTTCTTGGTGGTATGCACTTGATGATCGAAGAAAACGCCATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGAGTACCCGGGGATTAAAGTCTCGACGGCCGAAGCCATGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCGGGCttcattc                                                >  2:224754/21‑251 (MQ=255)                                                                                                                                                                                   cGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCttcattc                                                >  2:39537/1‑251 (MQ=255)                                                                                                                                                                                                                      aCGACAACGTGGCACCGTGTTTGCTTGGTGGTATGCAGTTGATGAGCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGTCATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCcgcg             <  2:186588/251‑1 (MQ=255)                                                                                                                                                                                                                         ccAACGTGGTAACGTGTTTTCTCGGTGGTAAGCAGCTGATGAGCGAAGAAAACGACATCACCAGCCAGTAAGTGTCAGGGTTTCATGAGTGTCTGTGGGTGCTGGCGTATCCGGGTATTAAAGTCTCGACGGCAGAAGATAGGGCTATGTTACCGGCGCAGTATCGTCTCCAGGATTGCATTGCGCACGGGCGACATCTGGATGTTTTCATTCACGCCTGCTATTCCCTTCAGCCTGAGCTTGCCGCGAAg          <  1:175843/250‑1 (MQ=255)                                                                                                                                                                                                                           aaCGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCt        <  1:110781/251‑1 (MQ=255)                                                                                                                                                                                                                                   aCCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCCTCAGCCAGCAAGTGACAGGGTTTGATGAGTGGATGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGAGCAGAATCGACGCCAGGATTGCCTTGCGCACGGGCGACATCTGACTGGCTTCATTCACGCCTGCTATACGCGTCAGCCTGAGCATGCCGCGAAGATGATGaa  >  1:113189/1‑249 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             GGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAA  >  NC_000913/3006‑3480

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