Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,058,009	A→G	M9T (ATG→ACG)	glnA ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,058,009	0	A	G	100.0%	111.6 / NA	33	M9T (ATG→ACG)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/15);  total (18/15)

AGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCATCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGATAAAGCGAAATCTGTGCCAACTTTTAAATTG  >  NC_000913/4057875‑4058145                                                                                                                                       |                                                                                                                                         aggagCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCgtcgtc                                                                                                                                     >  2:76945/1‑140 (MQ=255)         tCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACgtg                                                                                                                             <  1:734756/140‑1 (MQ=255)         tCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCACAAATCAACAAAATTCGCTTCGTGCCCGTTCAGCGTCGTCAGTACgtg                                                                                                                             >  2:295580/1‑140 (MQ=255)            aaCATTTTGCCTTCTTCGAAGACTTCAGCGTTCACCTTATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTc                                                                                                                          <  1:295580/140‑1 (MQ=255)            aaCATTTTGCCTTCTTCGAAGACTTCAGCCTTCACCTTATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTc                                                                                                                          <  1:222607/140‑1 (MQ=255)                ttttGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCg                                                                                                                      <  1:156230/140‑1 (MQ=255)                ttttGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCg                                                                                                                      >  2:664822/1‑140 (MQ=255)                  ttGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGa                                                                                                                    <  2:537460/140‑1 (MQ=255)                    gCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACa                                                                                                                  <  1:176822/140‑1 (MQ=255)                    gCCTTCTTCGAAGAAATCAGCACTCACCTGATGAGCAAGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACa                                                                                                                  <  1:369538/140‑1 (MQ=255)                          ttcGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACttt                                                                                                            >  1:705060/1‑140 (MQ=255)                            cGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTaa                                                                                                          <  1:733888/140‑1 (MQ=255)                            cGAAGAATTCAGCATACACCTGATGCGCAAGGATAGTAACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTaa                                                                                                          <  1:356718/140‑1 (MQ=255)                                             aCCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCa                                                                                         >  1:421226/1‑140 (MQ=255)                                                  atgaGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCTGCGGACAAACTATAGCTCTCCGGGATTGGTCATGgtc                                                                                    >  1:250900/1‑140 (MQ=255)                                                   gagacAGGTATAGTGAGGTGCTATTCTTTACCTCTGGTAGCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGgtcg                                                                                   <  2:373284/136‑1 (MQ=255)                                                       cAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTcgtcgt                                                                               >  1:768074/1‑140 (MQ=255)                                                               gTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGa                                                                       >  1:6613/1‑140 (MQ=255)                                                               gTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTCAACTCTCATGGCTTCGTCATGCTCGTCGTGGTAACGa                                                                       >  1:376871/1‑140 (MQ=255)                                                                    gTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCt                                                                  >  1:139317/1‑140 (MQ=255)                                                                     tGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCACCTGGATCGGTCATTGTCGTCGTGGTAACgctgtctc                                                                 >  1:373284/1‑133 (MQ=255)                                                                              ttACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTAttt                                                        >  1:544362/1‑140 (MQ=255)                                                                                      aGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCgtg                                                <  1:664822/140‑1 (MQ=255)                                                                                       gTATCGGTGAAGCGCAAATCAACAAACTTCCCTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCgtgt                                               <  1:76945/140‑1 (MQ=255)                                                                                                aaGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGCTCGTCGTGCTAAGGAATTCTGCAATGATTTAGGCCGTGTCGCCGTaaa                                      >  2:223014/1‑140 (MQ=255)                                                                                                     cAAATCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGa                                 >  2:611895/1‑140 (MQ=255)                                                                                                       aaTCAACAAACTTCACTTCGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGATa                               >  2:680074/1‑140 (MQ=255)                                                                                                                acttcacttcGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGATAAAGCGAAAt                      <  1:611895/140‑1 (MQ=255)                                                                                                                  ttcacttcGTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGATAAAGCGAAATCt                    >  2:17728/1‑140 (MQ=255)                                                                                                                       ttcGTGCTCGTTCAGCGCCGTCTGTACGTGTTCAGGGGACATACTTTAGCTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGATAAAGCGAAATCTGTGcc               <  2:250900/140‑1 (MQ=255)                                                                                                                          gTGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGATAAAGCGAAATCTGTGCCAAc            >  2:332749/1‑140 (MQ=255)                                                                                                                           tGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTGTTTTGGCCTTGTCGCCTTAAAAAAGATAAACCGAAATCTGTGCCAACt           >  2:307491/1‑140 (MQ=255)                                                                                                                           tGCTCGTTCAGCGTCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGATAAAGCGAAATCTGTGCCAACt           >  2:765643/1‑140 (MQ=255)                                                                                                                                    aGCGTCGTCAGTGCGGGTTCAGGGGGCATGCTTTAACCCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGATAAAGCGAAATCTGTGCCAACTTTTAAATTg  <  1:307491/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         AGGAGCCGTCAAACATTTTGCCTTCTTCGAAGAATTCAGCATTCACCTGATGAGCAGGGATAGTGACGTGCTGTTCTTTACCTTTAGTATCGGTGAAGCGCAAATCAACAAACTTCACTTCGTGCTCGTTCAGCATCGTCAGTACGTGTTCAGCGGACATACTTTAACTCTCCTGGATTGGTCATGGTCGTCGTGGTAACGAAATCTGCAATTATTTTGGCCGTGTCGCCGTAAAAAAGATAAAGCGAAATCTGTGCCAACTTTTAAATTG  >  NC_000913/4057875‑4058145

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