Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,855,960	C→T	S4N (AGC→AAC)  *116* (TAG→TAA)	yidG ← yidH ←	inner membrane protein DUF202 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,855,960	0	C	T	100.0%	176.0 / NA	55	S4N (AGC→AAC)  *116* (TAG→TAA)	yidG yidH	inner membrane protein DUF202 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/32);  total (23/32)

CATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCA  >  NC_000913/3855827‑3856093                                                                                                                                      |                                                                                                                                      cATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATcc                                                                                                                                  >  2:800429/1‑139 (MQ=255)   tGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                <  2:75141/139‑1 (MQ=255)   tGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                <  2:284549/139‑1 (MQ=255)      ccGTCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAATGATGTGCGCTCCGGCTACAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCAt                                                                                                                             >  1:517573/1‑139 (MQ=255)      ccGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACTACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCAt                                                                                                                             >  2:499908/1‑139 (MQ=255)             gTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  1:199204/139‑1 (MQ=255)             gTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  2:77901/139‑1 (MQ=255)             gTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      >  1:413584/1‑139 (MQ=255)             gTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  2:413584/139‑1 (MQ=255)              tGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     <  1:499908/139‑1 (MQ=255)                     tgtgTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAgcac                                                                                                                      >  1:677302/1‑131 (MQ=255)                     tgtgTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAgcac                                                                                                                      <  2:677302/131‑1 (MQ=255)                         tttGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAAc                                                                                                          <  2:267507/139‑1 (MQ=255)                            gATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGAGGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGAt                                                                                                       <  2:661424/139‑1 (MQ=255)                               aGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTc                                                                                                                 >  2:404005/1‑126 (MQ=255)                               aGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTc                                                                                                                 <  1:404005/126‑1 (MQ=255)                                 ccAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCg                                                                                                  <  2:76658/139‑1 (MQ=255)                                    aCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                                                                                               <  2:598229/139‑1 (MQ=255)                                          tCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCa                                                                                         >  1:639168/1‑139 (MQ=255)                                           cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTTATGTGCGCTCCGGCTGCAGTCCGGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  1:797825/139‑1 (MQ=255)                                                 gcCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAatat                                                                                  <  2:142883/139‑1 (MQ=255)                                                      aGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     >  2:673804/1‑99 (MQ=255)                                                      aGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     <  1:673804/99‑1 (MQ=255)                                                         ccAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTaa                                                                          <  1:522051/139‑1 (MQ=255)                                                                  tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 <  2:639168/139‑1 (MQ=255)                                                                  tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 >  2:101590/1‑139 (MQ=255)                                                                      aaaCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCt                                                             <  2:34843/139‑1 (MQ=255)                                                                       aaCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTg                                                            <  1:101590/139‑1 (MQ=255)                                                                        aCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  1:800429/139‑1 (MQ=255)                                                                            cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                <  1:499401/66‑1 (MQ=255)                                                                            cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                >  2:499401/1‑66 (MQ=255)                                                                            cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgt                                                       >  2:724956/1‑139 (MQ=255)                                                                                       tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCa                                                                                                                              >  1:346562/1‑57 (MQ=255)                                                                                       tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCa                                                                                                                              <  2:346562/57‑1 (MQ=255)                                                                                              cGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCttt                                     >  1:47700/1‑139 (MQ=255)                                                                                                                 gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                                                                                               <  1:15024/62‑1 (MQ=255)                                                                                                                 gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                                                                                               >  2:15024/1‑62 (MQ=255)                                                                                                                 gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAAt                                                                         >  1:539211/1‑84 (MQ=255)                                                                                                                 gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAAt                                                                         <  2:539211/84‑1 (MQ=255)                                                                                                                 gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGc                                                              >  2:740856/1‑95 (MQ=255)                                                                                                                 gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGc                                                              <  1:740856/95‑1 (MQ=255)                                                                                                                 gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           >  1:623158/1‑98 (MQ=255)                                                                                                                 gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  2:623158/98‑1 (MQ=255)                                                                                                                 gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCa                  >  1:232856/1‑139 (MQ=255)                                                                                                                   aaTGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  2:88378/37‑1 (MQ=255)                                                                                                                   aaTGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      >  1:88378/1‑37 (MQ=255)                                                                                                                    aTGCGCCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  1:394051/95‑1 (MQ=255)                                                                                                                    aTGCGCCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           >  2:394051/1‑95 (MQ=255)                                                                                                                      gCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCg             >  1:289943/1‑139 (MQ=255)                                                                                                                              cTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAgc     <  1:771886/139‑1 (MQ=255)                                                                                                                                ttGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAgcgc   >  2:24079/1‑139 (MQ=255)                                                                                                                                 tGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCa  <  1:24079/139‑1 (MQ=255)                                                                                                                                 tGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCa  <  2:289943/139‑1 (MQ=255)                                                                                                                                  gCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCg                                                                                                  >  2:182094/1‑42 (MQ=255)                                                                                                                                  gCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCg                                                                                                  <  1:182094/42‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      CATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCA  >  NC_000913/3855827‑3856093

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