Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,520	A→T	P240P (CCA→CCT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,520	0	A	T	96.7%	191.8 / ‑6.3	63	P240P (CCA→CCT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  major base T (24/35);  minor base G (0/2);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGC‑GGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCA  >  NC_000913/3275‑3759                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   gACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGCAAGATGTTATCGCTGAACCCTACCGTGAACGGTTCCTGCCAGGCtt                                                                                                                                                                                                                                               >  2:124718/1‑251 (MQ=255)       atcaGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggc                                                                                                                                                                                                                                          >  1:198553/1‑250 (MQ=255)                cAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggcgcggca                                                                                                                                                                                                                                <  1:258088/251‑1 (MQ=255)                  aGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGCGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggcgcggcagg                                                                                                                                                                                                                              >  1:202828/1‑251 (MQ=255)                   gTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATAGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaagcgcggcaggc                                                                                                                                                                                                                             >  2:256240/1‑251 (MQ=255)                                 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)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   GACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGC‑GGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTG‑ACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGG‑TTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCA  >  NC_000913/3275‑3759

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