Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,358,782	C→A	S1353Y (TCC→TAC)	gltB →	glutamate synthase subunit GltB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,358,782	0	C	A	100.0%	162.0 / NA	47	S1353Y (TCC→TAC)	gltB	glutamate synthase subunit GltB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (25/22);  total (25/22)

CGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTCCGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCT  >  NC_000913/3358645‑3358914                                                                                                                                          |                                                                                                                                     cGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACg                                                                                                                                    <  2:1123474/140‑1 (MQ=255)  gCCATTCGTCCCCCGGTTGGTTCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACgg                                                                                                                                   <  2:1170171/140‑1 (MQ=255)      ttCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCt                                                                                                                               <  1:1199382/140‑1 (MQ=255)                 gttggttCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGt                                                                                                                    >  1:638455/1‑140 (MQ=255)                    ggttCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGa                                                                                                                 >  1:1229968/1‑140 (MQ=255)                    ggttCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGa                                                                                                                 <  2:1069630/140‑1 (MQ=255)                         cGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCa                                                                                                            <  1:1194067/140‑1 (MQ=255)                                   aGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACaa                                                                                                  <  1:1067278/140‑1 (MQ=255)                                   aGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACaa                                                                                                  <  2:68623/140‑1 (MQ=255)                                   aGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACaa                                                                                                  <  2:1229968/140‑1 (MQ=255)                                        tGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGtt                                                                                             >  2:563861/1‑140 (MQ=255)                                        tGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGtt                                                                                             >  1:977937/1‑140 (MQ=255)                                        tGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGtt                                                                                             >  1:1471984/1‑140 (MQ=255)                                          aagcaagcATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTtgt                                                                                           <  1:153271/140‑1 (MQ=255)                                            gcaagcaTTATCGGCAACTCCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGa                                                                                         <  1:891508/140‑1 (MQ=255)                                                 caTTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAatata                                                                                    >  2:672038/1‑140 (MQ=255)                                                          cAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGggtg                                                                           <  1:830226/140‑1 (MQ=255)                                                          cAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGggtg                                                                           <  1:563861/140‑1 (MQ=255)                                                            acacCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGggtggt                                                                         >  1:1055872/1‑140 (MQ=255)                                                                   cctgTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGGAAGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGGCt                                                                  >  2:339746/1‑140 (MQ=255)                                                                         tGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTc                                                            <  1:1226022/140‑1 (MQ=255)                                                                           gcgcGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGggtggt                                                                         >  2:142411/1‑125 (MQ=255)                                                                           gcgcGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGggtggt                                                                         <  1:142411/125‑1 (MQ=255)                                                                            cgcgACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTgg                                                         >  2:11285/1‑140 (MQ=255)                                                                                 ccGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTaaa                                                    >  1:1173025/1‑140 (MQ=255)                                                                                  cGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTaaaa                                                   <  2:1055872/140‑1 (MQ=255)                                                                                         cgtcTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGt                                            <  2:1169979/140‑1 (MQ=255)                                                                                           tcTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTa                                          >  1:290947/1‑140 (MQ=255)                                                                                           tcTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTa                                          >  1:21443/1‑140 (MQ=255)                                                                                                   ccGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGt                                  <  1:339746/140‑1 (MQ=255)                                                                                                    cGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGGGTTAACTTCGGTg                                 <  1:11285/140‑1 (MQ=255)                                                                                                    cGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGt                                            <  2:318633/129‑1 (MQ=255)                                                                                                    cGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGt                                            >  1:318633/1‑129 (MQ=255)                                                                                                      cAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCg                               >  2:635483/1‑140 (MQ=255)                                                                                                      cAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCg                               >  1:963721/1‑140 (MQ=255)                                                                                                      cAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCg                               >  1:445363/1‑140 (MQ=255)                                                                                                      cAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTATGGGTAAAACCGGCGTTAACTTCGGTGCg                               >  2:1104310/1‑140 (MQ=255)                                                                                                         gCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGc                            >  1:422880/1‑140 (MQ=255)                                                                                                             cgcgGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGa                        >  1:992822/1‑140 (MQ=255)                                                                                                             cgcgGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGa                        >  2:393836/1‑140 (MQ=255)                                                                                                              gcgGGTGAACGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGAc                       >  2:1250787/1‑140 (MQ=255)                                                                                                                       cGTTTCGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGttt              <  1:1250787/140‑1 (MQ=255)                                                                                                                            cGGCGTGCGTAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCtt         >  2:794555/1‑140 (MQ=255)                                                                                                                               cgtgcgtAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACg      <  1:992623/140‑1 (MQ=255)                                                                                                                                 tgcgtAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCt          <  2:1103536/134‑1 (MQ=255)                                                                                                                                 tgcgtAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCt          >  1:1103536/1‑134 (MQ=255)                                                                                                                                  gcgtAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtc   >  1:1430013/1‑140 (MQ=255)                                                                                                                                   cgtAACTACGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTtct  >  1:488729/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                     CGCCATTCGTCCTCCGGTTGGTTCCGCCTTCCGCAGCCATGAAGCAAGCATTATCGGCAACACCTGCCTGTATGGCGCGACCGGTGGTCGTCTGTATGCCGCAGGCCGCGCGGGTGAACGTTTCGGCGTGCGTAACTCCGGTGCTATCACCGTGGTAGAAGGCATTGGCGACAACGGTTGTGAATATATGACGGGTGGTATCGTCTGCATTCTGGGTAAAACCGGCGTTAACTTCGGTGCGGGCATGACCGGCGGTTTCGCTTACGTTCT  >  NC_000913/3358645‑3358914

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑