Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,423,401	A→G	intergenic (+207/+22)	yhdZ → / ← rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,401	0	A	G	100.0%	161.9 / NA	47	intergenic (+207/+22)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/18);  total (29/18)

GTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATT  >  NC_000913/3423257‑3423546                                                                                                                                                 |                                                                                                                                                  gTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGc                                                                                                                                             <  1:619987/151‑1 (MQ=37)   cTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTTCTTCTCAACGTTTAAAATCGCCTCCCGCGGATTTGTTCTTCTCCGGGGGGGGTTTACCGCAAAACAAAAAGTAAAAAAAAAAACCCCGTCCTCCGGCttggcc                                                                                                                                            >  2:530450/1‑145 (MQ=2)      gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCtttt                                                                                                                                        <  2:294592/151‑1 (MQ=37)                      gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCt                                                                                                                        >  1:614036/1‑151 (MQ=255)                      gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCt                                                                                                                        >  2:72443/1‑151 (MQ=255)                         cTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGc                                                                                                                     >  2:389400/1‑151 (MQ=255)                          tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGc                                                                                                                     >  2:65002/1‑150 (MQ=255)                            cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGt                                                                                                                  <  1:187769/151‑1 (MQ=255)                                           gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCAt                                                                                                   >  2:421290/1‑151 (MQ=255)                                            ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATg                                                                                                  >  1:7803/1‑151 (MQ=255)                                                     cGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                         >  2:156394/1‑151 (MQ=255)                                                       cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAcccc                                                                                        >  1:465635/1‑150 (MQ=255)                                                        aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcac                                                                                      <  1:589380/151‑1 (MQ=255)                                                           gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACt                                                                                   >  2:344694/1‑151 (MQ=255)                                                               aaaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCa                                                                               <  2:197492/151‑1 (MQ=255)                                                               aaaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCa                                                                               <  1:718211/151‑1 (MQ=255)                                                                 aTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATc                                                                             <  2:170829/151‑1 (MQ=255)                                                                  tCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCg                                                                            >  2:522075/1‑151 (MQ=255)                                                                             cggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCg                                                                 >  1:591779/1‑151 (MQ=255)                                                                                 tttGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTc                                                             >  1:5158/1‑151 (MQ=255)                                                                                      ccTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                                         <  1:217171/150‑1 (MQ=255)                                                                                      ccTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                        >  1:206134/1‑151 (MQ=255)                                                                                       cTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCCTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                       >  2:513971/1‑151 (MQ=255)                                                                                               gagaGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTACGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                                 >  2:325112/1‑149 (MQ=255)                                                                                                    tgtTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                          >  2:183321/1‑151 (MQ=39)                                                                                                      ttCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgggg                                        >  1:269096/1‑151 (MQ=38)                                                                                                         aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCa                                     <  1:4387/151‑1 (MQ=38)                                                                                                          ccGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc                                      >  2:326873/1‑149 (MQ=38)                                                                                                           cGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAgg                                   >  2:679151/1‑151 (MQ=38)                                                                                                            gACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAgg                                   >  2:393729/1‑150 (MQ=38)                                                                                                              caaacaaCAGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg                                <  1:38605/151‑1 (MQ=38)                                                                                                                 acaacaGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac                             <  2:640961/151‑1 (MQ=38)                                                                                                                 acaacaGATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac                             >  2:95780/1‑151 (MQ=38)                                                                                                                       gATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcgc                       <  1:641741/151‑1 (MQ=38)                                                                                                                       gATAAAACAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcgc                       <  2:473647/151‑1 (MQ=38)                                                                                                                             aCAAAAGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAc                  >  1:400994/1‑150 (MQ=38)                                                                                                                                aaaGGCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgcc              >  2:756413/1‑151 (MQ=25)                                                                                                                                  agtctcAGTATTCCGGCTGAGCATTTTGTTATATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgccgc            <  2:651724/146‑1 (MQ=11)                                                                                                                                    gCCCAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCa          <  2:102612/151‑1 (MQ=25)                                                                                                                                     cccAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATCTCTGACAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCGTGGGGTCAGGTGGGACCTCCGCGCTACTGCCGCCa          <  2:15198/150‑1 (MQ=255)                                                                                                                                      ccAGTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTATGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgccgcc           >  2:508376/1‑148 (MQ=255)                                                                                                                                         gTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTg                >  1:24826/1‑140 (MQ=37)                                                                                                                                         gTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTg                <  2:24826/140‑1 (MQ=37)                                                                                                                                         gTCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAggc       >  2:44326/1‑147 (MQ=19)                                                                                                                                          tCTTCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACaa     >  1:325426/1‑150 (MQ=38)                                                                                                                                            ttCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAAtt  >  1:466465/1‑151 (MQ=38)                                                                                                                                             tCCGGCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAAtt  <  1:414756/150‑1 (MQ=38)                                                                                                                                                 |                                                                                                                                                  GTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATT  >  NC_000913/3423257‑3423546

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