Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,912,061	T→C	E536E (GAA→GAG)	glmS ←	L‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,912,061	0	T	C	100.0%	142.7 / NA	42	E536E (GAA→GAG)	glmS	L‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/22);  total (20/22)

GATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAA  >  NC_000913/3911983‑3912137                                                                               |                                                                             gatgTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAAcctc                                                                            >  1:710003/1‑81 (MQ=255)    gtgCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAAcctctt                                                                          >  1:709995/1‑80 (MQ=255)       cATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCaa                                                                       <  2:709948/80‑1 (MQ=255)         tGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGt                                                                    <  1:709950/81‑1 (MQ=255)              tCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGAt                                                               <  1:709951/81‑1 (MQ=255)                gCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGAttt                                                             >  2:710000/1‑81 (MQ=255)                gCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGAttt                                                             >  1:710026/1‑81 (MQ=255)                    cTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGt                                                         <  2:709951/81‑1 (MQ=255)                     ttACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGtt                                                        >  1:710002/1‑81 (MQ=255)                       aCAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGtttt                                                      >  2:710014/1‑81 (MQ=255)                       aCAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGtt                                                        <  2:709952/79‑1 (MQ=255)                          aaaCCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTcc                                                   <  1:709954/81‑1 (MQ=255)                             cccGCATCCTGATCGGCGAAGACATACCACTGGCCGCCACGCGCGCGAACCTCTTCAATGTTGAATTTCAGTTTTTCCAGc                                                <  2:709954/81‑1 (MQ=255)                              ccGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCACTTCAATGTTGGATTTCAGTTTTTCCAGCa                                               <  2:709955/81‑1 (MQ=255)                                   tCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCg                                          <  2:709956/81‑1 (MQ=255)                                     cTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgtt                                        <  1:709958/81‑1 (MQ=255)                                       gATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgttgt                                      >  2:710006/1‑81 (MQ=255)                                       gATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgttgt                                      >  1:710052/1‑81 (MQ=255)                                         tCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTc                                    <  1:709960/81‑1 (MQ=255)                                         tCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTc                                    <  2:709958/81‑1 (MQ=255)                                           ggCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCgg                                  <  1:709961/81‑1 (MQ=255)                                            gCGAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGt                                 >  2:709990/1‑81 (MQ=255)                                              gAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTg                                >  1:710008/1‑80 (MQ=255)                                              gAAGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTg                                >  2:710015/1‑80 (MQ=255)                                                aGACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGcaa                             >  2:709961/1‑81 (MQ=255)                                                 gACATACAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGcaac                            <  1:709963/81‑1 (MQ=255)                                                      aCAACTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGcaacaacaa                       >  1:710014/1‑81 (MQ=255)                                                        aaCTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAAt                      >  2:710009/1‑80 (MQ=255)                                                          cTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAAc                   >  1:710022/1‑81 (MQ=255)                                                          cTGACCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAAc                   <  1:709964/81‑1 (MQ=255)                                                             aCCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAAcc                  <  1:709966/79‑1 (MQ=255)                                                             aCCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCgg                >  1:710023/1‑81 (MQ=255)                                                             aCCGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCg                 <  2:709964/80‑1 (MQ=255)                                                              ccGCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCgg                >  2:710023/1‑80 (MQ=255)                                                                gCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCat             <  2:709966/81‑1 (MQ=255)                                                                gCCACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCa              >  2:710005/1‑80 (MQ=255)                                                                  cACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCatat           <  1:709968/81‑1 (MQ=255)                                                                  cACGCGCGCGAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCatat           >  1:710025/1‑81 (MQ=255)                                                                    cgcgcgcgAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCg         <  1:709969/81‑1 (MQ=255)                                                                         gcgAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCAt     >  2:710012/1‑80 (MQ=255)                                                                         gcgAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATc    <  2:709969/81‑1 (MQ=255)                                                                          cgAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCa   <  2:709970/81‑1 (MQ=255)                                                                           gAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCaa  >  1:710004/1‑81 (MQ=255)                                                                           gAACCTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCaa  >  2:710017/1‑81 (MQ=255)                                                                               |                                                                             GATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAA  >  NC_000913/3911983‑3912137

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