Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,360,724	T→C	intergenic (+108/‑452)	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,724	0	T	C	100.0%	89.0 / NA	25	intergenic (+108/‑452)	gltD/gltF	glutamate synthase subunit GltD/periplasmic protein GltF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/16);  total (9/16)

GCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTC  >  NC_000913/3360651‑3360794                                                                          |                                                                       gCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAggg                                                                    <  1:637425/78‑1 (MQ=255) gCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAggg                                                                    <  1:637424/78‑1 (MQ=255)   gAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAAt                                                               <  1:637426/81‑1 (MQ=255)   gAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAAt                                                               >  2:637460/1‑81 (MQ=255)     aTCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAAttt                                                             <  2:637427/81‑1 (MQ=255)        aaaTTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTc                                                           >  2:637468/1‑80 (MQ=255)         aaTTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCaa                                                         >  2:637463/1‑81 (MQ=255)              cTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATaaa                                                     <  1:637428/80‑1 (MQ=255)                    aaGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCt                                              <  1:637429/81‑1 (MQ=255)                          gcgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTc                                        <  1:637430/81‑1 (MQ=255)                          gcgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTc                                        <  1:637431/81‑1 (MQ=255)                           cgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCa                                       <  2:637432/81‑1 (MQ=255)                           cgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTCCa                                       <  2:637433/81‑1 (MQ=255)                            gttgttTTCTGGCAAATTCGGCATGAGCGTTTGACATTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCAc                                      >  2:637474/1‑81 (MQ=255)                             ttgttTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCAcc                                     <  2:637434/81‑1 (MQ=255)                                ttttCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAg                                  >  2:637471/1‑81 (MQ=255)                                 tttCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGt                                 >  1:637472/1‑81 (MQ=255)                                       gCAAATGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtgat                            <  1:637435/80‑1 (MQ=255)                                         aaaTGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAg                         >  1:637470/1‑81 (MQ=255)                                         aaaTGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAg                         <  1:637436/81‑1 (MQ=255)                                           aTGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTg                       <  2:637437/81‑1 (MQ=255)                                            tGCGGCATGAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGa                      <  1:637438/81‑1 (MQ=255)                                                    gAGCGTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACAt               >  1:637466/1‑80 (MQ=255)                                                        gTTTGACCTTTCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAg           >  1:637477/1‑80 (MQ=255)                                                                  tCCTATAACAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTtctc  <  1:637439/79‑1 (MQ=255)                                                                          |                                                                       GCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTC  >  NC_000913/3360651‑3360794

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑