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%	195.7 / NA	60	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 (26/34);  total (26/34)

ATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCA  >  NC_000913/3855828‑3856093                                                                                                                                     |                                                                                                                                      aTGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  2:500026/139‑1 (MQ=255) aTGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  2:528439/139‑1 (MQ=255)  tGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                <  1:487456/139‑1 (MQ=255)  tGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCGGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                <  1:237908/139‑1 (MQ=255)     ccGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCAt                                                                                                                             >  1:515888/1‑139 (MQ=255)       gCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATac                                                                                                                           <  1:1034724/139‑1 (MQ=255)            gTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  1:549512/139‑1 (MQ=255)             taccAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGTGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     >  2:216937/3‑139 (MQ=255)             tGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     >  1:893480/1‑139 (MQ=255)                       gtTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATaa                                                                                                           >  1:911511/1‑139 (MQ=255)                           gATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                >  1:293436/1‑114 (MQ=255)                           gATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                <  2:293436/114‑1 (MQ=255)                            aTAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATc                                                                                                      >  1:1117810/1‑139 (MQ=255)                                ccAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCg                                                                                                  >  2:861687/1‑139 (MQ=255)                                   aCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                                                                                               <  1:194715/139‑1 (MQ=255)                                          cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  1:679613/109‑1 (MQ=255)                                          cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      >  2:679613/1‑109 (MQ=255)                                          cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  2:307537/139‑1 (MQ=255)                                          cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  2:1117810/139‑1 (MQ=255)                                                    tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAg                                                                              <  1:607461/139‑1 (MQ=255)                                                         cAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACAAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAAt                                                                         <  1:861687/139‑1 (MQ=255)                                                            cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAg                                                                                                                   <  1:944677/94‑1 (MQ=255)                                                            cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAg                                                                                                                   >  2:944677/1‑94 (MQ=255)                                                            cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGAt                                                                      >  2:751926/1‑139 (MQ=255)                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  2:21637/75‑1 (MQ=255)                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  2:384118/75‑1 (MQ=255)                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  2:534380/75‑1 (MQ=255)                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 >  1:534380/1‑75 (MQ=255)                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 >  1:384118/1‑75 (MQ=255)                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 >  1:21637/1‑75 (MQ=255)                                                                  aCGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaataa                                                                <  2:773561/139‑1 (MQ=255)                                                                       aCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  1:751926/139‑1 (MQ=255)                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     <  1:667010/77‑1 (MQ=255)                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     >  2:667010/1‑77 (MQ=255)                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        >  1:1080887/1‑106 (MQ=255)                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  2:1080887/106‑1 (MQ=255)                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgt                                                       >  1:650306/1‑139 (MQ=255)                                                                            gCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgtg                                                      <  2:220031/139‑1 (MQ=255)                                                                            gCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgtg                                                      <  1:476529/139‑1 (MQ=255)                                                                                 tGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATTCGGCATACAacac                                                                                                                      >  1:467081/1‑70 (MQ=255)                                                                                 tGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATTCGGCATACAacac                                                                                                                      <  2:467081/70‑1 (MQ=255)                                                                                   aTGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATa                                                                                                            <  1:590945/78‑1 (MQ=255)                                                                                   aTGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATa                                                                                                            >  2:590945/1‑78 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  2:1150858/64‑1 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             >  1:1150858/1‑64 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTa                                                                                >  1:699370/1‑77 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTa                                                                                <  2:699370/77‑1 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTaa                                                                          <  2:977776/83‑1 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTaa                                                                          >  1:977776/1‑83 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgtgt                                                     >  2:448039/1‑104 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgtgt                                                     <  1:448039/104‑1 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGt                                            >  1:681499/1‑113 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGt                                            <  2:681499/113‑1 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCct                                         >  2:568360/1‑116 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCct                                         <  1:568360/116‑1 (MQ=255)                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCa                  >  2:1189871/1‑139 (MQ=255)                                                                                                                   aTGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTg               <  1:1189871/139‑1 (MQ=255)                                                                                                                      cGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCttt                                     <  1:902960/114‑1 (MQ=255)                                                                                                                      cGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCttt                                     >  2:902960/1‑114 (MQ=255)                                                                                                                                tGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCa  <  2:893480/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      ATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCA  >  NC_000913/3855828‑3856093

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