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%	124.0 / NA	39	*127* (TAG→TAA)	yigG	PRK11371 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/21);  total (18/21)

TGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATATA  >  NC_000913/4002683‑4002947                                                                                                                                   |                                                                                                                                       tgctTTATATTCTGGTGTTGCCCCCCAGAATCCAGG‑AAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCATTCATATATTCCTTACTCATAAATACTCCATAATTATTGtttt                                                                                                                               >  2:736080/1‑139 (MQ=255)    tttATATTCTGTTGTTGCCCCCCAGAATCCAGG‑AAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATAGTTTtgat                                                                                                                            <  1:736080/139‑1 (MQ=255)             tGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGt                                                                                                                    <  1:806196/139‑1 (MQ=255)                       ccccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGc                                                                                                          >  2:974186/1‑139 (MQ=255)                                                   cATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAAc                                                                              >  1:112643/1‑139 (MQ=255)                                                   cATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAAc                                                                              >  1:511450/1‑139 (MQ=255)                                                     tAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGGACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCAtt                                                                                                       <  1:758199/112‑1 (MQ=255)                                                     tAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTGGATGAATCAGTAGGTGCAAGCAtt                                                                                                       >  2:758199/1‑112 (MQ=255)                                                               aTCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCaaaaaaa                                                                  >  1:610365/1‑139 (MQ=255)                                                                         ccATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAAcc                                                        <  1:469979/139‑1 (MQ=255)                                                                                         gAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCa                                        <  2:455228/139‑1 (MQ=255)                                                                                                    gatatatTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGc                             <  1:19061/138‑1 (MQ=255)                                                                                                    cATATATTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGc                             <  2:112643/139‑1 (MQ=255)                                                                                                     atatatTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGc                             >  1:361885/1‑138 (MQ=255)                                                                                                     atatatTCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGc                             <  2:361885/138‑1 (MQ=255)                                                                                                          ttCCTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCAt                       >  2:513016/1‑139 (MQ=255)                                                                                                            ccTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtg                     >  2:385283/1‑139 (MQ=255)                                                                                                             cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgt                    >  1:659501/1‑139 (MQ=255)                                                                                                             cTTACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgt                    <  2:573503/139‑1 (MQ=255)                                                                                                               tACTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGtga                  >  2:799526/1‑139 (MQ=255)                                                                                                                aCTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACAATTGCCTGAACAGGAAAAATCTTAGGCAATCATtgtaac                 >  2:19061/1‑136 (MQ=255)                                                                                                                 cTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgtg                   >  1:1111228/1‑136 (MQ=255)                                                                                                                 cTCATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATtgtg                   <  2:1111228/136‑1 (MQ=255)                                                                                                                   cATAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTgatgat              <  2:925002/139‑1 (MQ=255)                                                                                                                    aTAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATa             <  1:592974/139‑1 (MQ=255)                                                                                                                     tAAATGCTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATag            <  1:394588/139‑1 (MQ=255)                                                                                                                     tAAATACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATag            <  1:513016/139‑1 (MQ=255)                                                                                                                       aaTACTCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTCCCTGAACAGGCAAAACCTTCGGCTATCATTGTGATGATagag          <  1:740980/139‑1 (MQ=255)                                                                                                                       aaTACTCCATAATTATTGTTTCGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATACCAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTTTGATGagagag          <  2:659501/139‑1 (MQ=255)                                                                                                                            tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAgatgat     >  2:542085/1‑139 (MQ=255)                                                                                                                            tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATC‑AAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTc                                >  1:283998/1‑111 (MQ=255)                                                                                                                            tCCATAATTATTGTTTTGATGAATCAGTAGGTGCAAGAATTAGCATACTGAAAGTGGAAAAATAACAAATC‑AAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTc                                <  2:283998/111‑1 (MQ=255)                                                                                                                            tCCATAATTATTGTCTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag          <  2:477152/134‑1 (MQ=255)                                                                                                                            tCCATAATTATTGTCTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag          >  1:477152/1‑134 (MQ=255)                                                                                                                              cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGatat   >  1:117337/1‑139 (MQ=255)                                                                                                                              cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATC‑AAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGc                             <  1:908638/112‑1 (MQ=255)                                                                                                                              cATAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATC‑AAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGc                             >  2:908638/1‑112 (MQ=255)                                                                                                                               aTAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGatata  <  1:979812/139‑1 (MQ=255)                                                                                                                               aTAATTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGatata  <  2:1117960/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       TGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATATA  >  NC_000913/4002683‑4002947

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