Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,002,813	C→T	*127* (TAG→TAA)	yigG ←	PRK11371 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,002,813	0	C	T	100.0%	186.2 / NA	58	*127* (TAG→TAA)	yigG	PRK11371 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/29);  total (29/29)

CTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATAT  >  NC_000913/4002694‑4002946                                                                                                                        |                                                                                                                                      cTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAg                                                                                                                    <  2:622663/139‑1 (MQ=255)  tGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                   >  2:611438/1‑139 (MQ=255)    gtgtTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAg                                                                                                                 >  1:1305023/1‑139 (MQ=255)       ttGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                              <  1:1069087/139‑1 (MQ=255)        tGCCCTCCAAAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCAGCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGc                                                                                                             <  2:1409974/139‑1 (MQ=255)         gCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTtgatga                                                                                                                         <  2:466564/126‑1 (MQ=255)         gCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTtgatga                                                                                                                         >  1:466564/1‑126 (MQ=255)                         gAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGt                                                                                                               <  2:637295/120‑1 (MQ=255)                         gAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGt                                                                                                               >  1:637295/1‑120 (MQ=255)                             aCCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGt                                                                                        <  1:611438/139‑1 (MQ=255)                                aTGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGa                                                                                     >  1:1023031/1‑139 (MQ=255)                                        cATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTtgatg                                                                                                                          >  2:151418/1‑94 (MQ=255)                                        cATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTtgatg                                                                                                                          <  1:151418/94‑1 (MQ=255)                                        cATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                              <  1:21320/106‑1 (MQ=255)                                        cATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTg                                                                                                              >  2:21320/1‑106 (MQ=255)                                        cATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAAc                                                                             >  1:1423001/1‑139 (MQ=255)                                           aaaaaTTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAAt                                                                                                                       <  2:53816/94‑1 (MQ=255)                                           aaaaaTTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAAt                                                                                                                       >  1:53816/1‑94 (MQ=255)                                            aaaaTTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAAt                                                                         >  1:386179/1‑139 (MQ=255)                                                                   tttCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGc                                                  <  2:386179/139‑1 (MQ=255)                                                                                   aTCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTg                                                   >  2:654800/1‑122 (MQ=255)                                                                                   aTCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTg                                                   <  1:654800/122‑1 (MQ=255)                                                                                   aTCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATc                                  >  2:222487/1‑139 (MQ=255)                                                                                    tCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCt                                 <  1:222487/139‑1 (MQ=255)                                                                                             tatTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATc                        <  2:1423001/139‑1 (MQ=255)                                                                                                 ccTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt                      <  2:307867/137‑1 (MQ=255)                                                                                                 ccTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt                      >  1:307867/1‑137 (MQ=255)                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAg                                                                                                                    >  1:1682632/1‑42 (MQ=255)                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAg                                                                                                                    <  2:1682632/42‑1 (MQ=255)                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                   >  1:1317655/1‑43 (MQ=255)                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                   <  2:1317655/43‑1 (MQ=255)                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgt                   >  1:1692510/1‑139 (MQ=255)                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgt                   <  1:1244356/139‑1 (MQ=255)                                                                                                  cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgt                   >  2:1244356/1‑139 (MQ=255)                                                                                                    tACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGtga                 <  2:1692510/139‑1 (MQ=255)                                                                                                     aCTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCaaaaaa                                                                  >  2:1183858/1‑89 (MQ=255)                                                                                                     aCTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCaaaaaa                                                                  <  1:1183858/89‑1 (MQ=255)                                                                                                         aTAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATa            >  2:532687/1‑139 (MQ=255)                                                                                                            aaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag         <  2:1305023/139‑1 (MQ=255)                                                                                                               aCTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAgatg      <  1:532687/139‑1 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAatcatc                                                            <  2:757940/83‑1 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAatcatc                                                            >  1:757940/1‑83 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAatcatc                                                            >  1:533348/1‑83 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAatcatc                                                            <  2:533348/83‑1 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAAcc                                                       >  2:704322/1‑88 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAAcc                                                       <  1:704322/88‑1 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCt                                 >  2:262403/1‑110 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCt                                 <  1:262403/110‑1 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt                      <  2:1609227/121‑1 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt                      >  1:1609227/1‑121 (MQ=255)                                                                                                                 tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAgatgat    >  1:1522141/1‑139 (MQ=255)                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGaa                                             <  2:83999/96‑1 (MQ=255)                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGaa                                             >  1:83999/1‑96 (MQ=255)                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt                      >  1:445845/1‑119 (MQ=255)                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt                      <  2:445845/119‑1 (MQ=255)                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGatat  >  2:1530247/1‑139 (MQ=255)                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTACCATACt                                                                                              >  1:1667203/1‑47 (MQ=255)                                                                                                                   cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTACCATACt                                                                                              <  2:1667203/47‑1 (MQ=255)                                                                                                                        |                                                                                                                                      CTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATAT  >  NC_000913/4002694‑4002946

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