Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,765,234	G→A	G129G (GGG→GGA)	yibJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,765,234	0	G	A	100.0%	128.6 / NA	47	G129G (GGG→GGA)	yibJ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/25);  total (22/25)

GCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTT  >  NC_000913/3765201‑3765306                                  |                                                                         gctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                                <  2:4417713/73‑1 (MQ=255)  ctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                               <  2:1876939/74‑1 (MQ=255)  ctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                               >  1:581762/3‑76 (MQ=255)   tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                              >  2:3284304/2‑76 (MQ=255)   tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                              <  1:2041333/75‑1 (MQ=255)   tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                              >  1:257649/2‑76 (MQ=255)   tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                              >  1:4849582/2‑76 (MQ=255)     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCa                            <  1:2018697/76‑1 (MQ=255)      ggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt                           >  2:3641620/1‑76 (MQ=255)      ggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt                           <  1:740953/76‑1 (MQ=255)      ggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt                           <  1:3271094/76‑1 (MQ=255)         ggTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                        >  1:1569721/1‑76 (MQ=255)          gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt                       >  1:2371302/1‑76 (MQ=255)               tCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                  >  2:1554173/1‑76 (MQ=255)               tCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                  >  2:1529569/1‑76 (MQ=255)               tCACGCGGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTGGAATCCGGTCATAAATGTAGa                  <  2:144519/76‑1 (MQ=17)               tCACGCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                  <  2:102475/76‑1 (MQ=255)                cACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt                 <  2:878657/76‑1 (MQ=255)                 aCTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                <  2:4234530/76‑1 (MQ=255)                 aCTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                >  2:2051096/1‑76 (MQ=255)                 aCTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCCGTCATCAATGCAGATc                >  2:2079098/1‑76 (MQ=11)                  cTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATcc               <  2:4331616/76‑1 (MQ=255)                  cTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATcc               >  2:1346848/1‑76 (MQ=255)                   tCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCg              <  1:3830746/76‑1 (MQ=255)                   tCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCg              <  2:4769006/76‑1 (MQ=255)                   tCAGGAGCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGTCATAAACGCAGATCCg              >  2:2151892/1‑76 (MQ=14)                    cAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGc             <  2:3770701/76‑1 (MQ=255)                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCa            >  1:4787494/1‑76 (MQ=255)                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCa            <  2:4787494/76‑1 (MQ=255)                      ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa           <  1:1051240/76‑1 (MQ=255)                      ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCACAAATGTAGAGTCGCaa           >  1:4828677/1‑76 (MQ=17)                        aTCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGTAAgg         >  2:1129829/1‑76 (MQ=255)                        aTCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAgg         >  1:4492510/1‑76 (MQ=255)                          ccGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGtt       <  2:1170655/76‑1 (MQ=255)                          ccGATTGGAATGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGtt       <  2:916682/76‑1 (MQ=21)                          ccGAGTGGACTGAAGGGGGGATGGAATTTTTATCAGGATCCGTTGAATCCGGTCATAAAGGTAGATCCGCAAGGtt       <  2:1217099/76‑1 (MQ=21)                           cGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGTAAGGttt      <  1:1129830/76‑1 (MQ=21)                           cGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt      <  1:3259831/76‑1 (MQ=255)                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg     >  2:3014377/1‑76 (MQ=255)                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg     >  1:2812667/1‑76 (MQ=255)                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg     >  1:2758909/1‑76 (MQ=255)                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg     >  1:1366147/1‑76 (MQ=255)                              ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATACATGTAGATCCGCAAGGTTTGGt   <  1:887520/76‑1 (MQ=255)                              ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt   <  2:3746896/76‑1 (MQ=255)                              ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt   <  2:4127268/76‑1 (MQ=255)                              ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt   >  1:1866694/1‑76 (MQ=255)                               tGGACTGAGGGGGGGATTGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGtt  <  1:2462140/76‑1 (MQ=25)                                  |                                                                         GCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTT  >  NC_000913/3765201‑3765306

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