Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,547	C→T	V249V (GTC→GTT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,547	0	C	T	98.3%	182.6 / ‑4.7	60	V249V (GTC→GTT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/37);  total (21/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

TGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGC‑GGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC  >  NC_000913/3307‑3792                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        tgaGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCCGGC‑GGGTGCGGGAATCg                                                                                                                                                                                                                                               >  1:4161/1‑251 (MQ=255)     gggCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGgcgc                                                                                                                                                                                                                                           <  1:237920/250‑1 (MQ=255)          gtgGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAATATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGTCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTag                                                                                                                                                                                                                                      <  1:155158/251‑1 (MQ=255)          gtgGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTag                                                                                                                                                                                                                                      <  1:102511/251‑1 (MQ=255)           tgGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGGagc                                                                                                                                                                                                                                     >  1:93988/1‑251 (MQ=255)                            ggggATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATctccggct                                                                                                                                                                                                                    <  1:407959/251‑1 (MQ=255)                              ggATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGTGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATctccggctcc                                                                                                                                                                                                                  <  1:169126/251‑1 (MQ=255)                               gATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATctccggctccg                                                                                                                                                                                                                 <  2:107765/251‑1 (MQ=255)                                           gACGGCAGAAGCCAGGGCAATTTTACCGGCGCAGTAGCGCCTCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGGCATCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGt                                                                                                                                                                                                     <  2:230664/251‑1 (MQ=255)                                                    agcccgggcttgtttACCGGCGCCGTATCGCCGCCAGGATTGAATGGCGCACGGGCGACATCTTGCTGGCTGTATTGACGCCTGATAATCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATTAAAGATGTTATCGCTGAACCCTACCGGGAACTATTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCTGGT‑GCGGAAATCGGCGCGGTAGCGAGCGTTATCGCCGGCTCCGGGCCGACCTGGTTCGCTCtgt                                                                                                                                                                                            <  2:437424/239‑1 (MQ=255)                                                       cAGGGATATTTTGTCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCAAGGGCGACATCTGGCTGGCTGCATTCACGACTGCTATTCCCGTCAGCCTGAGCTTGCCGAGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCtgtgtg                                                                                                                                                                                         <  2:171171/251‑1 (MQ=255)                                                             tATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGGAATCGGCGCGGGAGCGAGCGGTATCTCCGGCCCCGGCCCGACCTTGTTCGCTCTGTGGG‑GCAAGc                                                                                                                                                                                  >  1:469198/1‑251 (MQ=255)                                                              aTTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGGGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGGCCCTGTTTGCCCTGTGTG‑GCAAGcc                                                                                                                                                                                 >  1:114657/1‑251 (MQ=255)                                                                    ccGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGGTCCGGCCCGACCTTGTTCGCTCTGTGTG‑AAAAGCCgggtac                                                                                                                                                                           >  1:316394/1‑247 (MQ=255)                                                                      ggCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCgggtaccg                                                                                                                                                                         >  1:114963/1‑245 (MQ=255)                                                                          cAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑GCAAGCCGGATACCGCCCa                                                                                                                                                                     <  1:452226/251‑1 (MQ=255)                                                                              aTCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAgcg                                                                                                                                                                  <  1:36261/250‑1 (MQ=255)                                                                                  ccgccAGGATTGCATTGCGCACGGGCGTCTTCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGGTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTg                                                                                                                                                             <  2:229175/251‑1 (MQ=255)                                                                                    gacaggagtgcatttctcaatggCGACATCTGGCTGGCTTCATTCCCGCCTGCTAGTCCCGTCAGCCTGATCTTGCCGCGAATCTGATGAAAGATGTTATCGCGGAACCCTACCGTGAACGGTTACTTCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCTGTGCGTTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGcc                                                                                                                                                           <  1:248046/230‑1 (MQ=255)                                                                                           ttgcattgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGACAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCACGACCTTGTTCGCTCTGGGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTggg                                                                                                                                                    >  1:273976/1‑250 (MQ=255)                                                                                               attgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTggg                                                                                                                                                <  2:93988/251‑1 (MQ=255)                                                                                                  gcgcACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGGAATTGGCGCGGGAGCGAGCGGTGTCTCCGGCTCCGGGCCGACCCTTTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCACGCGCGTTGCCGACCGGtttgggt                                                                                                                                              >  2:248164/1‑245 (MQ=255)                                                                                                      ccGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGC‑GGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAAc                                                                                                                                         <  1:41443/250‑1 (MQ=255)                                                                                                       cGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGGAATCCGCGCGGGGGCGAGCGGTATCTCCGGCCCCGGCCCGCCCTTGTTCGCTCTGTGGG‑ACAAGCCGGGTACCGCCCCGCGCGCTGCCCGCTGGGTGGGGAAGAAc                                                                                                                                         >  2:396522/1‑250 (MQ=255)                                                                                                                 ctggctggCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCCGCCCCACCTTTGTCGCTCTGTGGGGACA‑GCCGGATAACGCCCAGCGCGTTGTCGACTGGGTGGGGTAGAACTACCTCCaaa                                                                                                                               >  2:89776/1‑250 (MQ=255)                                                                                                                  tggctggcTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGCAATCGGCGCGGTAGCGAGCGGTATCTCCCGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGGCCAGCGCGTTGCCGTCTGGTTGGGgaagaaataactggcaaat                                                                                                                             >  1:59038/1‑238 (MQ=255)                                                                                                                   ggctggCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCaaa                                                                                                                               <  1:193852/244‑1 (MQ=255)                                                                                                                      tggCTTCATTCACGCCTGATATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGCTGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCCGGC‑GGTTGAGGAAATCGGCGCGGTAGCGAGCGGTATATCCGGCTCCGGGCCGACCCTGTTCGCTCTTTGTG‑ACAAGCCGGTACCCGCCCCGCGCGTTGGCGCCCGGTTGGGGAAGAACTAACTCCAAACTCacg                                                                                                                         >  2:338859/2‑249 (MQ=255)                                                                                                                       ggCTTCATTCACGCCGGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAg                                                                                                                          <  1:333563/249‑1 (MQ=255)                                                                                                                            cattcaCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGACAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGGAATCGGCCCGGGAGGGAGGGGGATCTCCGGCTCCGGCCCGCACTTGTTCTCTCTGTGGG‑ACAAACCAGATAACGCCCAGGGCGTTGCCGCATGGTTGGGGCAGAACTACCAGCCAACTCAGGGAGG‑tt                                                                                                                  >  2:221367/1‑251 (MQ=255)                                                                                                                            cattcaCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTAGCGCTGAGCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGTCCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑tt                                                                                                                  <  1:47213/251‑1 (MQ=255)                                                                                                                                caCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTg                                                                                                               <  1:248164/250‑1 (MQ=255)                                                                                                                                    ccTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTGATCGCGGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCat                                                                                                          <  2:129051/251‑1 (MQ=255)                                                                                                                                          aTTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACAACCTGCAAAATCAGGAAGG‑TTTTGTTCATAttt                                                                                                      <  1:30623/249‑1 (MQ=255)                                                                                                                                          aTTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCCGCGCGGGAGCGAGCGGGATCTCCCGCTCCGGCCCGACCCTGTTCGCTCTGTGGG‑ACAAGCCGGATACCGCCCCGCGCGTTGCCGGCCGGGTGGGGAAGAAACACCCGCCAAATCAGGAAGGGTTTTTTCCATATTTg                                                                                                     >  2:395888/1‑251 (MQ=255)                                                                                                                                             cccGGCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCgg                                                                                                 <  2:469198/251‑1 (MQ=255)                                                                                                                                                     ccTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGGAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCGGCCCGCCCTTGTTCGCCCTGTGGG‑ACAAGCCGGATACCGCCCAGCGCGGTGCCCGCCGGGTGGGGAAAAACTTCCCGCCAAAACAGGGAGG‑GTTTGTTCATATTTTCCggacgggatc                                                                                         >  2:108617/1‑243 (MQ=255)                                                                                                                                                      cTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGCCCGCCCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGGTGCCGACTGGTTGGGGAAGAACTACCTGCAAAATCAGGAAGG‑GTTTGTTCATATTTGCCGGCTGGATAcg                                                                                        >  1:282959/1‑251 (MQ=255)                                                                                                                                                       tGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCTGCGCGGTAGCGAGCGGTATCTCCGTCTCCGGcc                                                                                                                                                                                                              <  2:473950/134‑1 (MQ=255)                                                                                                                                                       tGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGcc                                                                                                                                                                                                              >  1:473950/1‑134 (MQ=255)                                                                                                                                                            ttGCCGCGAAGCTGATGAAAGATGTTAGCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACggcggg                                                                                   <  2:10184/250‑1 (MQ=255)                                                                                                                                                             tGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACggcggg                                                                                   <  1:34373/249‑1 (MQ=255)                                                                                                                                                             tGCCGCGAAGCGGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGGG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACggcgggc                                                                                  <  2:4161/250‑1 (MQ=255)                                                                                                                                                                  cgAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGa                                                                            <  1:383077/251‑1 (MQ=255)                                                                                                                                                                   gAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGGTGCCGAATGGGTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑GTTTGTTCATATTTGCCGGCTGGATACggcgggaggcccgg                                                                           >  1:264586/1‑243 (MQ=255)                                                                                                                                                                     aGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTa                                                                         <  1:232557/251‑1 (MQ=255)                                                                                                                                                                      gCTGATGAAAGATGTTATCGGTGAACCCTACCGTGAACGGTTACGGCCTGGCGTCCGGCAGGCGCGGCTGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTAc                                                                        <  2:368541/251‑1 (MQ=255)                                                                                                                                                                      gCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGGCTGGGTGGGGAAGAAATAACTGCAAAATAAGGAAGG‑GTTTGTTCATATTTGCCGGCTGGATACggcgggggccccgagac                                                                        >  1:197126/1‑240 (MQ=255)                                                                                                                                                                         gatgaAAGATGTTAGCGCTGAACCCTTCCGTGAACGGTTACGGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACt                                                                       <  2:235354/249‑1 (MQ=255)                                                                                                                                                                                    ttATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑GTTTGTTCATATTTGCCGGCTGGATACggcgggggggcggggagctggaaacacaaa                                                           >  1:98969/1‑226 (MQ=255)                                                                                                                                                                                      aTCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACggcgggcggtcggcgtattggaaaataaaaca                                                         >  1:187649/1‑226 (MQ=255)                                                                                                                                                                                        cGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGaaa                                                       <  1:396522/250‑1 (MQ=255)                                                                                                                                                                                             atcccgaccttGAACGGTTACTGCCTGGCGTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAGCTCCGGCTCCGGCCCTACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTGCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTAc                                                 <  2:181350/241‑1 (MQ=255)                                                                                                                                                                                                ccTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAAt                                              <  1:40595/251‑1 (MQ=255)                                                                                                                                                                                                   aCCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATc                                             <  1:395888/249‑1 (MQ=255)                                                                                                                                                                                                    ccGTGAACGGTTACTGCCGGGCTGCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTTCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCCGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGa                                          <  2:469846/251‑1 (MQ=255)                                                                                                                                                                                                      gTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGaaa                                        <  1:130505/251‑1 (MQ=255)                                                                                                                                                                                                             gTTACTGCCGGGCTTCCGTCAGGCGCGGCAGTC‑GGTTGCGGAAATCGGCTCGTGAGCGAGCGTTAGCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTTGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCTCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATca                                   <  1:285796/249‑1 (MQ=255)                                                                                                                                                                                                                           tCCGGCAGGCGCGGCAGGC‑GGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAg                    <  1:89776/250‑1 (MQ=255)                                                                                                                                                                                                                                            gc‑gGTTGCGGAAATCGGCGCGGTAGCGGGCGGTACCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAAc  <  2:344302/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        TGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGC‑GGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC  >  NC_000913/3307‑3792

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: ‑