Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	821,506	G→A	*238* (TAG→TAA)	ybhM →	BAX Inhibitor‑1 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	821,506	0	G	A	100.0%	173.0 / NA	54	*238* (TAG→TAA)	ybhM	BAX Inhibitor‑1 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/28);  total (26/28)

TGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAGCCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGATAAAGT  >  NC_000913/821378‑821633                                                                                                                                 |                                                                                                                                tgGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGc                                                                                                                       <  1:1244971/139‑1 (MQ=255)  ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                      <  1:1087607/139‑1 (MQ=255)   gAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGcac                                                                                                                     <  2:1290571/139‑1 (MQ=255)    aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGcaca                                                                                                                    <  2:174130/139‑1 (MQ=255)           ttCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                      >  2:693656/1‑130 (MQ=255)           ttCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                      <  1:693656/130‑1 (MQ=255)                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTaa                                                                                                 >  1:121296/1‑139 (MQ=255)                                 aTTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGcaca                                                                                                                    >  1:1183281/1‑110 (MQ=255)                                 aTTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGcaca                                                                                                                    <  2:1183281/110‑1 (MQ=255)                                          ctATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTaa                                                                                                                           >  1:243865/1‑94 (MQ=255)                                          ctATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTaa                                                                                                                           <  2:243865/94‑1 (MQ=255)                                             ttACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCAttt                                                                           >  2:824820/1‑139 (MQ=255)                                             ttACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCAttt                                                                           >  2:407738/1‑139 (MQ=255)                                               aCTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCAttttt                                                                         <  2:510467/139‑1 (MQ=255)                                                cTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTc                                                                        >  1:586261/1‑139 (MQ=255)                                                cTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTc                                                                        >  1:293612/1‑139 (MQ=255)                                                     ttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTa                                                                                                                            >  2:744498/1‑82 (MQ=255)                                                     ttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTa                                                                                                                            <  1:744498/82‑1 (MQ=255)                                                     ttATCGGCTTCTCTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCg                                                                   >  1:92235/1‑139 (MQ=255)                                                           gCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                      <  2:609828/82‑1 (MQ=255)                                                           gCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                      >  1:609828/1‑82 (MQ=255)                                                                  tGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGcc                                                      >  2:746469/1‑139 (MQ=255)                                                                   gggATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCt                                                     <  2:293612/139‑1 (MQ=255)                                                                      aTACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGc                                                  <  2:520347/139‑1 (MQ=255)                                                                         cTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCttt                                               >  2:1275154/1‑139 (MQ=255)                                                                          ttGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTc                                              >  1:695207/1‑139 (MQ=255)                                                                          ttGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTc                                              >  2:189357/1‑139 (MQ=255)                                                                            gCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCa                                            >  2:124821/1‑139 (MQ=255)                                                                            gCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCa                                            >  1:179166/1‑139 (MQ=255)                                                                              cgcgATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGc                                          >  2:188051/1‑139 (MQ=255)                                                                               gcgATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                      <  2:258728/62‑1 (MQ=255)                                                                               gcgATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                      >  1:258728/1‑62 (MQ=255)                                                                                 gatagataTAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAAc                                       <  1:861858/139‑1 (MQ=255)                                                                                  aTAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACa                                      <  1:1275154/139‑1 (MQ=255)                                                                                    aGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAg                                           >  2:1079681/1‑132 (MQ=255)                                                                                    aGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAg                                           <  1:1079681/132‑1 (MQ=255)                                                                                           aCGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATAAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGa                             <  1:824820/139‑1 (MQ=255)                                                                                                tGTCTGGCAAAGGCATATGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAgc                        >  1:823498/1‑139 (MQ=255)                                                                                                tGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                      <  1:751831/45‑1 (MQ=255)                                                                                                tGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                      >  2:751831/1‑45 (MQ=255)                                                                                                tGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAgc                        <  1:635595/139‑1 (MQ=255)                                                                                                tGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAgc                        <  1:189357/139‑1 (MQ=255)                                                                                                    tGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACa                                      >  1:519321/1‑121 (MQ=255)                                                                                                    tGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACa                                      <  2:519321/121‑1 (MQ=255)                                                                                                       cAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGaa                                                                                                          >  2:935679/1‑50 (MQ=255)                                                                                                       cAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGaa                                                                                                          <  1:935679/50‑1 (MQ=255)                                                                                                        aaaGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTc                                                                               <  2:173432/76‑1 (MQ=255)                                                                                                        aaaGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTc                                                                               >  1:173432/1‑76 (MQ=255)                                                                                                           ggCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCcagcagcag             <  1:188051/139‑1 (MQ=255)                                                                                                                aCGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGa        <  2:121296/139‑1 (MQ=255)                                                                                                                aCGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGa        <  1:746469/139‑1 (MQ=255)                                                                                                                aCGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGa        <  2:823498/139‑1 (MQ=255)                                                                                                                aCGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGa        <  1:124821/139‑1 (MQ=255)                                                                                                                      ttttttCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGATAAAGt  >  2:234612/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                TGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAGCCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGATAAAGT  >  NC_000913/821378‑821633

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