Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,622	A→T	E274D (GAA→GAT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,622	0	A	T	97.7%	116.4 / NA	43	E274D (GAA→GAT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (7/35);  total (8/35)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTG  >  NC_000913/3378‑3853                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        cgcAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGcc                                                                                                                                                                                                                                   <  1:612735/251‑1 (MQ=255)        ccgccgcCAGGATTTCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGcgcg                                                                                                                                                                                                                            <  2:239528/250‑1 (MQ=255)               cAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGa                                                                                                                                                                                                                     <  1:199385/251‑1 (MQ=255)               cAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATAGGCGCGGTATCGAGCGGTATCTCCGGCTCCGGCCCGACCTTTTTCGCTCTTTGAGACAAACCGGATACCGCCCAGCGCGTTGCCgc                                                                                                                                                                                                                     >  1:24179/1‑250 (MQ=255)                        attgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGGTACCGCCCAGCGCGTTGCCGACTGGTTggg                                                                                                                                                                                                            >  1:222923/1‑251 (MQ=255)                         ttgcGCACGGGCGACAGTTGGATGGCTTCAGTCAAGCCTGCTATTCCCGACAGCCTGAGTTTGCCGCGAAGCTGATGAAAGATGG‑TATCGCTGCACCCTGCCGTGAACGGTTACTGCCTGGCTGCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTggg                                                                                                                                                                                                            <  2:65022/250‑1 (MQ=255)                           gcgcCCGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTACTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACATTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTaa                                                                                                                                                                                                         <  1:502385/251‑1 (MQ=255)                               aCGGGCGAGATCTGGTTGGGTTCATTCACGCCTGCTATCCCCGTCAGCCTGTGCTGTCCGCGAAGCTGATGAAAGATGT‑TAACGCGGAACCCGAGCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGGAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTTACAAGCCGGATACCGCCCAGCGCGTTGCCTACTGGTTGGGTAAGAAc                                                                                                                                                                                                     <  2:292297/251‑1 (MQ=255)                                            ggctggCTTCATTCACGCCTGCTATTCCCGTTAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATc                                                                                                                                                                                        <  1:236460/247‑1 (MQ=255)                                                            gCCTGCTATTCCCGTCAGCCTGAGCTGGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCa                                                                                                                                                                        <  1:139006/251‑1 (MQ=255)                                                               tGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGTAGCGAGCGGTATCTCCCGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTTGGTAAGAACTACCTGCAAAATCAGGAAGGGTTTGTTCatat                                                                                                                                                                     >  1:365009/1‑251 (MQ=255)                                                                  ttttCCCGGCAGCCTTGGCTTGCCGCGAATCTGATGAAAGA‑GTGTATCGCTGAACCCTTCCGTGAACGGTTACTGCCTTGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGTAATCGGCGCTGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTGCCTGCAAAATCAGGAAGGTTTTGTTCATAttt                                                                                                                                                                   <  2:409996/248‑1 (MQ=255)                                                                      cccGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCgg                                                                                                                                                              <  2:416242/251‑1 (MQ=255)                                                                       ccGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGc                                                                                                                                                             <  1:568262/251‑1 (MQ=255)                                                                                 gAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTACGGCTCCGGCCCGCCCTTGTTCGCTCTGTGTGGAAAGCAGGATACCTCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTAACTGCCAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATAcggc                                                                                                                                                   >  1:45209/1‑251 (MQ=255)                                                                                      tGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACggcgggc                                                                                                                                               <  1:218711/250‑1 (MQ=255)                                                                                      tGCCGCGAAGCTGATGAAAGATGG‑TGGCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCTGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCGCCGGCTCCGGCCCGACCTTGTTCGCTCTGGGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACggcgggc                                                                                                                                               <  2:573873/250‑1 (MQ=255)                                                                                            gAAGCTGATGAAAGATGT‑TATCGCTTAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAg                                                                                                                                        <  1:439399/251‑1 (MQ=255)                                                                                            gAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAACCAGCGGGAACTCCGGCTCCGGCCCGGCCTTGTTCGCTCTGTTGGACAAACCGGGTACCGCCCCGCGCGTTTGCGACCGGTTGGGTAAGAAATACCCGCCAAATCAGGGAGGGTTTGGTCATATTTGCCgggcgggtacgggggggggaagg                                                                                                                                         >  2:637205/1‑229 (MQ=255)                                                                                             aaGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCGCGAGt                                                                                                                                       <  1:93018/251‑1 (MQ=255)                                                                                              aGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATTGGCGCGGTAGCGAGAGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACCAGACGGAAACCGCCCCGCGCGGTGCCGCCTGGATGGGGAAGAACATCCTGCCAACACAGGAAGGGTTTGTTCATATTTTCCGGCTGGGAACggcggggagagcg                                                                                                                                         >  2:315958/1‑241 (MQ=255)                                                                                                  gatgaAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACt                                                                                                                                    <  1:303490/249‑1 (MQ=255)                                                                                                     gaAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGaaa                                                                                                                               <  1:330807/251‑1 (MQ=255)                                                                                                        aGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACt                                                                                                                            <  1:186893/251‑1 (MQ=255)                                                                                                             t‑tATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAACCAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGGCTGGGTGGGGAAGAACTACCTGCAAAAACAGGAAGGGTTTTTTCATATTTGCCGGCTGGATACGGGGGGGGCAcagggacatgaaaaatgaaa                                                                                                                        >  1:300655/1‑231 (MQ=255)                                                                                                                         cccTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACaa                                                                                                            <  1:192179/251‑1 (MQ=255)                                                                                                                          ccgaCCGTGTAAGGTTACTGCCTGGCGTCCGGCAGGCGCGGCCGTCGGTTGCGTAAATCGGCGCGGTAGCTAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTTGTTGTGTAAGACCTACCTGCAAAATCAGGAAGGTTTTGATCATATTTGCCGGCTGGATACGGCGGGCTCTCGAGTACTGGAAAACTAAATGAAACTCTACAAt                                                                                                           <  2:345240/248‑1 (MQ=255)                                                                                                                          ccTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAAt                                                                                                           <  1:248952/251‑1 (MQ=255)                                                                                                                             aCCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACTGCGGGCGCACGAGTACTGGAAAACTAAATGAAACGCTACAATCt                                                                                                         <  1:553685/250‑1 (MQ=255)                                                                                                                                           cggcctggCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACaatga                                                                                          <  2:222923/245‑4 (MQ=255)                                                                                                                                             gcctggCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACaatga                                                                                          <  1:143563/245‑4 (MQ=255)                                                                                                                                              cctggcatCCGGCTGGCGCGGCAGGCGGTTGCGTAAAGCTGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTGGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTTAGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGAGACGGCGGGCGCACGAGTACTGGAAAACGAAATGAAACTCGACAATCTGAAAGATcacgttgagca                                                                                       <  2:278036/244‑9 (MQ=255)                                                                                                                                                          caggcgcggcaggcgGTTGCGTAAATCGGCTCGTTAGCGGGCGGTATCTCCTGCTCCGGCCCGACCTTGTTCGCGCTGTGTTACAAGCCGGATACCGCCCAGCGCGTTGCCGTCTGGTTGGGTAAGATCTACCTGCAAAATCAGGAAGGTTTTGTGCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAACGAATCTTTCCAATCTGAAAGATCGCAATGAGCAgggaatctgtgc                                                                           <  2:281530/251‑11 (MQ=255)                                                                                                                                                                           tgcgGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAg                                                           <  1:31102/249‑1 (MQ=255)                                                                                                                                                                           ggggcaaaTCTGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGGTCGGCCTTTGGCGCTCGGTGTGACAAGCCGGGTACTGTCCAGGGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAAGCAGGAAGGTTGTGTTCACATGTGCCGGCTGGATACGGCGGGCACACGAGTACTGGAAAACTGAATGAAACTCTACAATCGGAAAGATCACAATGAGCAGGTCAGCTTGGCGTAAGCCGTAACCCAgg                                                          <  2:352715/246‑1 (MQ=255)                                                                                                                                                                                    cggcgctggagtgagcgGTATCTCCGGTCCCGGCCCGATCTTGTTCGGTCTTTGTGACAGGCGGTATCCCGCCCGGTGCGGTGCCTGCGGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTGTTGTTCATACTTGCCGGCTGGATACGGCGTGCGCACGAGTACTGGAAAACTAAATGAAACACTACAATCTGAAAGAGCGCAAGGAGCAGGGCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCa                                                 <  2:447401/239‑1 (MQ=255)                                                                                                                                                                                          ggTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATc                                           <  1:637205/251‑1 (MQ=255)                                                                                                                                                                                                     gTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGttt                                <  1:315958/251‑1 (MQ=255)                                                                                                                                                                                                           tctgcaccggcCCGTCCTTGTTCGCTCTTTGGGACAATCCGGATTCCGCCCAGCGCGTTGCCGACTGGTTGTGTAAGAACTACCTGCAAAAGCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAACCTCTGCAATCTTAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGttttt                              <  2:557738/241‑1 (MQ=255)                                                                                                                                                                                                                   ggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCaaaa                                                                                                                                                                                          <  1:440996/83‑1 (MQ=255)                                                                                                                                                                                                                   ggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCaaaa                                                                                                                                                                                          >  2:440996/1‑83 (MQ=255)                                                                                                                                                                                                                              gttggctatgtgtgACAAGCCGGATACAGCCTATCGTGGTGCCGACTGGGTGGGTAAGAACTACTTGCAAAACCAGGAAGGTTTTGATCATATTTGCCGGCTGGATACGGCGGGCGCAAGATTACGGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGCCAGCGTTGCGCAAGCCGGAACCCAGGGGGTGGGCAAAAATCATGGGCTGTGTTTTCCGCACGACCTGCCTGAATTc        <  2:190210/242‑1 (MQ=255)                                                                                                                                                                                                                                     ctGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctg  <  1:399093/249‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        CGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTG  >  NC_000913/3378‑3853

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