Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,057,170	C→T	A289T (GCA→ACA)	glnA ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,057,170	0	C	T	100.0%	127.3 / NA	40	A289T (GCA→ACA)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/20);  total (20/20)

CATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTGCGCTACGTTG  >  NC_000913/4057039‑4057303                                                                                                                                    |                                                                                                                                      cATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTg                                                                                                                                >  2:376878/1‑139 (MQ=255)      ccGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGcc                                                                                                                           >  2:392127/1‑139 (MQ=255)            cTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGaa                                                                                                                     >  1:415084/1‑139 (MQ=255)            cTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGaa                                                                                                                     <  2:303856/139‑1 (MQ=255)            cTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGaa                                                                                                                     <  1:765430/139‑1 (MQ=255)                   gCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTa                                                                                                              >  1:213950/1‑139 (MQ=255)                                   ttATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGaa                                                                                                                     <  2:613215/116‑1 (MQ=255)                                   ttATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGaa                                                                                                                     >  1:613215/1‑116 (MQ=255)                                    tataAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACag                                                                                             <  2:401553/139‑1 (MQ=255)                                           gTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACAtg                                                                                      <  2:527722/139‑1 (MQ=255)                                             tggtggTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCt                                                                                                                          <  1:485832/101‑1 (MQ=255)                                             tgatggtCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCACCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCt                                                                                                                          >  2:485832/4‑101 (MQ=255)                                              ggtggtCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTgg                                                                                   >  2:613603/1‑139 (MQ=255)                                               gtggtCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCt                                                                                                                          <  1:174557/99‑1 (MQ=255)                                               gtggtCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCt                                                                                                                          >  2:174557/1‑99 (MQ=255)                                                   tCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTg                                                                              <  1:376878/139‑1 (MQ=255)                                                     gggTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCa                                                                            <  1:613603/139‑1 (MQ=255)                                                                    ttAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTAt                                                             >  1:791923/1‑139 (MQ=255)                                                                                    gtTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGTGCCGTTATCACCGAACATCGGttt                                             >  1:468990/1‑139 (MQ=255)                                                                                              gccgccAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCggctccggtatgcac                                   <  2:791923/139‑14 (MQ=255)                                                                                               ccgccAATGTAGTACAGCGCCTGCTCAGACAGTCCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACa                                                                                        <  1:565991/85‑1 (MQ=255)                                                                                               ccgccAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACa                                                                                        >  2:565991/1‑85 (MQ=255)                                                                                                cgccAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGc                                                                                                          <  1:488381/66‑1 (MQ=255)                                                                                                cgccAATGTAGTACAGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGc                                                                                                          >  2:488381/1‑66 (MQ=255)                                                                                                              aGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCt                                                                                                                          <  2:37170/36‑1 (MQ=255)                                                                                                              aGCGCCTGCTCAGACAGACCTGTGTATTTGTCGCCt                                                                                                                          >  1:37170/1‑36 (MQ=255)                                                                                                                cgcCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTgg                                                                                   >  2:614946/1‑73 (MQ=255)                                                                                                                cgcCTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTgg                                                                                   <  1:614946/73‑1 (MQ=255)                                                                                                                  ccTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACag                                                                                             >  1:600167/1‑61 (MQ=255)                                                                                                                  ccTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACag                                                                                             <  2:600167/61‑1 (MQ=255)                                                                                                                  ccTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACa                                                                                              >  2:716093/1‑60 (MQ=255)                                                                                                                  ccTGCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACa                                                                                              <  1:716093/60‑1 (MQ=255)                                                                                                                     gCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCAcc                                                         >  2:675612/1‑94 (MQ=255)                                                                                                                     gCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCAcc                                                         <  1:675612/94‑1 (MQ=255)                                                                                                                     gCTCAGACAGACCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTg            >  2:158984/1‑139 (MQ=255)                                                                                                                            cagacCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTGCGCTACg     <  2:454471/139‑1 (MQ=255)                                                                                                                               acCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTGCGCTACGTtg  >  2:416538/1‑139 (MQ=255)                                                                                                                               acCTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTGCGCTACGTtg  <  1:158984/139‑1 (MQ=255)                                                                                                                                ccTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCg                                                                       <  2:474558/69‑1 (MQ=255)                                                                                                                                ccTGTGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCg                                                                       >  1:474558/1‑69 (MQ=255)                                                                                                                                    |                                                                                                                                      CATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTGCGCTACGTTG  >  NC_000913/4057039‑4057303

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