Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,312,830	C→T	D21N (GAC→AAC)	yciB ←	IspA family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,830	0	C	T	100.0%	185.4 / NA	57	D21N (GAC→AAC)	yciB	IspA family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/28);  total (29/28)

GCGTCAAGCCACCGAAGACGACCACCAGAACAAAAGTGATCAGGGCCATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTCATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAATCAACATATACAGGCGGAACAGA  >  NC_000913/1312695‑1312953                                                                                                                                        |                                                                                                                            gCGTCAAGCCACCGAAGACGACCACCAGAACAAAAGTGATCAGGGCCATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTtata                                                                                                                          >  1:874504/1‑139 (MQ=255)  cGTCAAGCCACCGAAGACGACCACCAGAACAAAAGTGATCAGAGCCATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATaa                                                                                                                         >  1:1105872/1‑139 (MQ=255)    tCAAGCCACCGAAGACGACCACCAGAACAAAAGTGATCAGGGCCATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAAt                                                                                                                       <  2:1474448/139‑1 (MQ=255)                    gacCACCAGAACAAAAGTGATCAGGGCCATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGaaa                                                                                                       <  2:1105872/139‑1 (MQ=255)                            gAACAAAAGTGATCAGGGCCATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAAc                                                                                               <  2:933981/139‑1 (MQ=255)                               cAAAAGTGATCAGGGCCATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCACATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAg                                                                                            <  2:173788/139‑1 (MQ=255)                                            ggCCATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCa                                                                               >  1:1275620/1‑139 (MQ=255)                                              ccATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAAc                                                                                               >  1:109772/1‑121 (MQ=255)                                              ccATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAAc                                                                                               <  2:109772/121‑1 (MQ=255)                                                 tcttctCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGaaa                                                                          >  2:1567070/1‑139 (MQ=255)                                                     ctcAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGc                                                                      <  2:1275620/139‑1 (MQ=255)                                                            ttACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCAttt                                                               >  1:382290/1‑139 (MQ=255)                                                             tACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCt                                                                     >  1:1189081/1‑132 (MQ=255)                                                             tACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCt                                                                     <  2:1189081/132‑1 (MQ=255)                                                                gAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACg                                                           <  2:382290/139‑1 (MQ=255)                                                                              cTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCt                                                                     <  2:610739/115‑1 (MQ=255)                                                                              cTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCt                                                                     >  1:610739/1‑115 (MQ=255)                                                                                      aaGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAAcgcg                                                                                                          <  1:1164945/70‑1 (MQ=255)                                                                                      aaGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAAcgcg                                                                                                          >  2:1164945/1‑70 (MQ=255)                                                                                      aaGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAg                                     >  1:216550/1‑139 (MQ=255)                                                                                        gCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCg                                                     <  2:768628/121‑1 (MQ=255)                                                                                        gCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCg                                                     >  1:768628/1‑121 (MQ=255)                                                                                        gCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGt                                   >  2:60087/1‑139 (MQ=255)                                                                                           caATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTAt                                >  1:1080282/1‑139 (MQ=255)                                                                                              tCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCt                                                                     <  1:21257/99‑1 (MQ=255)                                                                                              tCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCt                                                                     >  2:21257/1‑99 (MQ=255)                                                                                              tCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCa                                               >  1:430631/1‑121 (MQ=255)                                                                                              tCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCa                                               <  2:430631/121‑1 (MQ=255)                                                                                              tCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTgg                             >  2:1429411/1‑139 (MQ=255)                                                                                                  cGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGaa                         <  1:60087/139‑1 (MQ=255)                                                                                                  cGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGaa                         <  2:874504/139‑1 (MQ=255)                                                                                                    tGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAAt                        >  1:806321/1‑138 (MQ=255)                                                                                                    tGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAAt                        <  2:806321/138‑1 (MQ=255)                                                                                                      gCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTaaaa                                                                                    <  2:54031/76‑1 (MQ=255)                                                                                                      gCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTaaaa                                                                                    >  1:54031/1‑76 (MQ=255)                                                                                                      gCGACGATCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAATCaa                     >  1:1218450/1‑139 (MQ=255)                                                                                                         acgaTCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCaa                                              <  1:1469490/111‑1 (MQ=255)                                                                                                         acgaTCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCaa                                              >  2:1469490/1‑111 (MQ=255)                                                                                                         acgaTCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCG‑AAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCa                                               <  2:1040951/109‑1 (MQ=255)                                                                                                         acgaTCAGCGCCGCAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCG‑AAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCa                                               >  1:1040951/1‑109 (MQ=255)                                                                                                                      cAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACg                                                           <  1:1473796/85‑1 (MQ=255)                                                                                                                      cAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACg                                                           >  2:1473796/1‑85 (MQ=255)                                                                                                                      cAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTgg                             >  1:1079563/1‑115 (MQ=255)                                                                                                                      cAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAATCAACATATACAGGCGGAAc     >  2:1334709/1‑139 (MQ=255)                                                                                                                      cAGTAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAATCAACATATACAGGCGGAAc     >  1:287583/1‑139 (MQ=255)                                                                                                                      cAGTAGCCGCATACATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTgg                             <  2:1079563/115‑1 (MQ=255)                                                                                                                         tAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAATCAACATATACAGGCGGAACaga  <  2:287583/139‑1 (MQ=255)                                                                                                                         tAGCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAATCAACATATACAGGCGGAACaga  <  1:1429411/139‑1 (MQ=255)                                                                                                                           gCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAg                                                                                            <  1:918235/47‑1 (MQ=255)                                                                                                                           gCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAg                                                                                            >  2:918235/1‑47 (MQ=255)                                                                                                                           gCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAATCAACATATACAGGCgg        >  1:659846/1‑131 (MQ=255)                                                                                                                           gCCGCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAATCAACATATACAGGCgg        <  2:659846/131‑1 (MQ=255)                                                                                                                              gCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCg                                                     <  2:1167077/83‑1 (MQ=255)                                                                                                                              gCATAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCg                                                     >  1:1167077/1‑83 (MQ=255)                                                                                                                                aTAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTg                                                                       >  2:13580/1‑63 (MQ=255)                                                                                                                                aTAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTg                                                                       <  1:13580/63‑1 (MQ=255)                                                                                                                                aTAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGaa                         <  2:1527716/109‑1 (MQ=255)                                                                                                                                aTAGATGTTATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGaa                         >  1:1527716/1‑109 (MQ=255)                                                                                                                                        |                                                                                                                            GCGTCAAGCCACCGAAGACGACCACCAGAACAAAAGTGATCAGGGCCATCTTCTCAACCTTACGAAAGCGAACCCAGCTATAGATAAGCACAATCGCCGTGGCGACGATCAGCGCCGCAGTAGCCGCATAGATGTCATAAATCTTGTAAAACGCGAAAAAGACAACCAGCGGTAAAAAATCAAGAAACTGCTTCATTTTACGATTCCGTAATCAAATGCTTAAGGTTATTGGCGAATCAACATATACAGGCGGAACAGA  >  NC_000913/1312695‑1312953

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