Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,360,951	C→T	intergenic (+335/‑225)	gltD → / → gltF	glutamate synthase subunit GltD/periplasmic protein GltF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,360,951	0	C	T	100.0%	50.6 / NA	16	intergenic (+335/‑225)	gltD/gltF	glutamate synthase subunit GltD/periplasmic protein GltF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/8);  total (8/8)

GAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAA  >  NC_000913/3360773‑3361186                                                                                                                                                                                   |                                                                                                                                                                                                                                            gAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTc                                                                                                                                                                     <  1:347023/251‑1 (MQ=255)     cacaTCAAGCAATTTCTCATTTACAATCCGTCTCAATTCCTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCCTGTAAATTGTCTCGATATTAATATCCAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACAt                                                                                                                                                                 >  2:45114/1‑251 (MQ=255)     cacaTCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACAt                                                                                                                                                                 >  1:123847/1‑251 (MQ=255)          cAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGAGTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATACATTATCCTTAATTATCTATATATCTTCTTATTTTTCAGCGCAAACATGCATGTAATGAACAATTATCTCCACAAACCACCTCTCCAGCACCTTATTATTCACATGATc                                                                                                                                                             <  2:174778/250‑1 (MQ=255)                      cATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAAt                                                                                                                                                <  2:360098/251‑1 (MQ=255)                       aTTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAAtt                                                                                                                                               >  1:303029/1‑251 (MQ=255)                                     tcAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATAc                                                                                                                                 >  1:257080/1‑251 (MQ=255)                                                                       tCTGCCAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACACCATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTACATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGttat                                                                                               >  2:181595/1‑251 (MQ=255)                                                                          gCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGttattt                                                                                             <  2:14167/250‑1 (MQ=255)                                                                                              aTTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGat                                                                         >  2:81006/1‑250 (MQ=255)                                                                                                 ccATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGatatat                                                                     <  2:257080/251‑1 (MQ=255)                                                                                                                       aaTATACAAAATATGAATATAAAAAACCAATATATTATCCCTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTttttatt                                                <  2:126025/250‑1 (MQ=255)                                                                                                                                                          tatCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCat             <  1:186155/250‑1 (MQ=255)                                                                                                                                                           atCCTTAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCatat           <  2:94196/251‑1 (MQ=255)                                                                                                                                                                tAATTATCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGtttt      >  1:291146/1‑251 (MQ=255)                                                                                                                                                                     aTCTATATATTTTCTTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaa  >  1:210521/1‑250 (MQ=255)                                                                                                                                                                                   |                                                                                                                                                                                                                                            GAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAA  >  NC_000913/3360773‑3361186

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