Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,588,041	G→T	G275G (GGC→GGA)	yjiA ←	P‑loop guanosine triphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,588,041	0	G	T	100.0%	112.9 / NA	33	G275G (GGC→GGA)	yjiA	P‑loop guanosine triphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/19);  total (14/19)

CTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAATCCA  >  NC_000913/4587908‑4588174                                                                                                                                      |                                                                                                                                      cTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTgg                                                                                                                                 >  2:353465/1‑140 (MQ=255)       tCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAAcagc                                                                                                                           <  2:826795/140‑1 (MQ=255)          tCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCCCCCCACGGCCTGTCCCAGTCGGCGCCGTAGAGGCGCTGGACCCCCCGGAACAGCAgg                                                                                                                        <  2:635817/140‑1 (MQ=255)              ccccGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGt                                                                                                                    >  1:1548/1‑140 (MQ=255)              ccccGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGt                                                                                                                    >  1:543894/1‑140 (MQ=255)                ccGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTa                                                                                                                  >  2:197785/1‑140 (MQ=255)                    aaaCGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGtt                                                                                                              <  1:436729/140‑1 (MQ=255)                          aGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGt                                                                                                        <  1:218775/140‑1 (MQ=255)                              cGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAAt                                                                                                    <  2:68706/140‑1 (MQ=255)                               gAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATc                                                                                                   <  2:1548/140‑1 (MQ=255)                               gAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATc                                                                                                   <  1:495132/140‑1 (MQ=255)                                  tctcTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCcac                                                                                                >  1:238308/1‑140 (MQ=255)                                    tctTCCTCCGGCAGCTGGGTCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAg                                                                                              <  2:238308/140‑1 (MQ=255)                                    tctTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAg                                                                                              <  1:114207/140‑1 (MQ=255)                                    tctTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAg                                                                                              <  1:643983/140‑1 (MQ=255)                                    tctTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAg                                                                                              <  2:535024/140‑1 (MQ=255)                                         ctccGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATcc                                                                                         <  1:536888/140‑1 (MQ=255)                                         ctccGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATcc                                                                                         <  1:353465/140‑1 (MQ=255)                                              gCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCCGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTg                                                                                    >  2:5789/1‑140 (MQ=255)                                               cAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCCGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGt                                                                                   <  1:5789/140‑1 (MQ=255)                                                            aTAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGttt                                                                      >  2:640889/1‑140 (MQ=255)                                                                   ccATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATccgcc                                                               <  1:608656/140‑1 (MQ=255)                                                                     aTCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGTATCCCTTTGGAACGGAGCAGTTTATCCGCCGa                                                             >  1:800456/1‑140 (MQ=255)                                                                         tGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTcc                                                         >  1:218430/1‑140 (MQ=255)                                                                                             ccccACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCcagcagc                                                   >  1:817555/1‑126 (MQ=255)                                                                                             ccccACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCcagcagc                                                   <  2:817555/126‑1 (MQ=255)                                                                                                     ccTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGaaa                             <  2:218430/140‑1 (MQ=255)                                                                                                        gTCCCAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACtt                          >  2:463974/1‑140 (MQ=255)                                                                                                          cccAGTCGGCGCTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCg                        <  1:640889/140‑1 (MQ=255)                                                                                                                   cgcTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCt               >  2:538542/1‑140 (MQ=255)                                                                                                                    gcTGTAGAGGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTa              <  1:538542/140‑1 (MQ=255)                                                                                                                          gagGCGCTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGt        <  1:463974/140‑1 (MQ=255)                                                                                                                                cTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAATCCa  >  2:128546/1‑140 (MQ=255)                                                                                                                                cTGGACTCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAATCCa  >  2:379384/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                      CTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAATCCA  >  NC_000913/4587908‑4588174

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