Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,995,832	G→T	G84C (GGC→TGC)	yigA →	conserved protein YigA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,832	0	G	T	100.0%	123.8 / NA	36	G84C (GGC→TGC)	yigA	conserved protein YigA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (20/16);  total (20/16)

ATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTC  >  NC_000913/3995703‑3995961                                                                                                                                  |                                                                                                                                  aTACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTa                                                                                                                         >  2:93194/1‑140 (MQ=255) aTACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTa                                                                                                                         >  1:193045/1‑140 (MQ=255) aTACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTa                                                                                                                         >  2:47841/1‑140 (MQ=255) aTACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTa                                                                                                                         >  1:520993/1‑140 (MQ=255) aTACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTa                                                                                                                         >  1:426857/1‑140 (MQ=255)        tgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCg                                                                                                                      >  1:196860/1‑136 (MQ=255)        tgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCg                                                                                                                      <  2:196860/136‑1 (MQ=255)                   gTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCg                                                                                                                      >  2:104660/1‑125 (MQ=255)                   gTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCg                                                                                                                      <  1:104660/125‑1 (MQ=255)                           cACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGtc                                                                                               >  1:380659/1‑140 (MQ=255)                              cGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCa                                                                                            <  2:380659/140‑1 (MQ=255)                                     ttGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACcgccgc                                                                                     >  1:184495/1‑140 (MQ=255)                                                   caTGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTcgacg                                                                       <  2:334744/140‑1 (MQ=255)                                                   caTGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTcgacg                                                                       <  1:93194/140‑1 (MQ=255)                                                      ggCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGata                                                                    <  1:47841/140‑1 (MQ=255)                                                            cgcACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCCCCGCCGCCAGCAGTCTCGACGATATGCTGa                                                              <  2:520993/140‑1 (MQ=255)                                                            cgcACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGa                                                              >  1:568057/1‑140 (MQ=255)                                                                             aTGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTgg                                             >  2:585266/1‑140 (MQ=255)                                                                             aTGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTgg                                             >  2:24044/1‑140 (MQ=255)                                                                               gTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGgc                                           >  2:601870/1‑140 (MQ=255)                                                                               gTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGgc                                           >  2:551931/1‑140 (MQ=255)                                                                               gTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGgc                                           >  2:35704/1‑140 (MQ=255)                                                                                  cTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGgcgcg                                        <  1:601870/140‑1 (MQ=255)                                                                                  cTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGgcgcg                                        <  1:551931/140‑1 (MQ=255)                                                                                  cTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGgcgcg                                        <  1:266312/140‑1 (MQ=255)                                                                                  cTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGgcgcg                                        <  2:193045/140‑1 (MQ=255)                                                                                   tGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGgcgcgc                                       >  2:501190/1‑140 (MQ=255)                                                                                     gAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCTAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGa                                     >  2:552110/1‑140 (MQ=255)                                                                                       agagAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGAtc                                   <  2:184495/140‑1 (MQ=255)                                                                                           aaCATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCTGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCgg                               >  1:307620/1‑140 (MQ=255)                                                                                             cATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGcc                             >  2:291063/1‑140 (MQ=255)                                                                                                  cgcTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGAtc                                   <  2:213745/129‑1 (MQ=255)                                                                                                  cgcTGTTGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGAtc                                   >  1:213745/1‑129 (MQ=255)                                                                                                       ttGATGGAACAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGc                   <  1:501190/140‑1 (MQ=255)                                                                                                              aaCAGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTg            <  1:291063/140‑1 (MQ=255)                                                                                                                 aGGCTATCGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTgcgc         <  1:585266/140‑1 (MQ=255)                                                                                                                        cGCCAACGAATGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTc  >  2:360080/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                  ATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTC  >  NC_000913/3995703‑3995961

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