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%	195.0 / NA	62	*238* (TAG→TAA)	ybhM	BAX Inhibitor‑1 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/34);  total (28/34)

AGAGTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAGCCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGATAAAG  >  NC_000913/821371‑821632                                                                                                                                        |                                                                                                                               agagTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCg                                                                                                                             <  1:1327925/139‑1 (MQ=255)   agTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTa                                                                                                                           >  1:610107/1‑139 (MQ=255)         ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                     <  1:1019001/139‑1 (MQ=255)         ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                     <  2:786082/139‑1 (MQ=255)           aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGcaca                                                                                                                   <  1:769100/139‑1 (MQ=255)           aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGcaca                                                                                                                   <  1:220761/139‑1 (MQ=255)              ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATAc                                                                                                                >  2:1141863/1‑139 (MQ=255)                 gTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTaa                                                                                                             <  2:610107/139‑1 (MQ=255)                     gTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGaa                                                                                                         <  1:164891/139‑1 (MQ=255)                                        aTTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATc                                                                                      >  1:57252/1‑139 (MQ=255)                                                tctATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTc                                                                              >  2:260760/1‑139 (MQ=255)                                                   aTTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGcaca                                                                                                                   <  2:598358/99‑1 (MQ=255)                                                   aTTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGcaca                                                                                                                   >  1:598358/1‑99 (MQ=255)                                                    ttACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCAttt                                                                          >  2:1194934/1‑139 (MQ=255)                                                    ttACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCAttt                                                                          >  2:1257207/1‑139 (MQ=255)                                                          ttttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTa                                                                                                                           >  2:118569/1‑84 (MQ=255)                                                          ttttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTa                                                                                                                           <  1:118569/84‑1 (MQ=255)                                                          ttttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGcc                                                                    >  2:619266/1‑139 (MQ=255)                                                          ttttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGcc                                                                    >  2:1222328/1‑139 (MQ=255)                                                          ttttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGcc                                                                    >  2:28543/1‑139 (MQ=255)                                                          ttttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGcc                                                                    >  1:352131/1‑139 (MQ=255)                                                          ttttATCGGCTTCTTTGGGATACTAGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTaa                                                                                                    >  2:1454717/1‑107 (MQ=255)                                                          ttttATCGGCTTCTTTGGGATACTAGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTaa                                                                                                    <  1:1454717/107‑1 (MQ=255)                                                            ttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCg                                                                  >  2:258806/1‑139 (MQ=255)                                                            ttATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCg                                                                  >  1:814977/1‑139 (MQ=255)                                                              atgggCTTCTTCGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTa                                                                                                                           >  1:851374/4‑80 (MQ=255)                                                              aTCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTa                                                                                                                           <  2:851374/80‑1 (MQ=255)                                                              aTCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTa                                                                >  2:1178703/1‑139 (MQ=255)                                                                      cttTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGTCGCTTCATTTTTCGCCCGTAGCttcttc                                                        >  2:881180/1‑139 (MQ=255)                                                                      cttTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACGTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCttcttc                                                        >  2:877256/1‑139 (MQ=255)                                                                            gATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTgg                                                  <  1:260760/139‑1 (MQ=255)                                                                                cTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTaa                                                                                                    <  1:1137360/85‑1 (MQ=255)                                                                                cTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTaa                                                                                                    >  2:1137360/1‑85 (MQ=255)                                                                                cTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACAt                                                                                       >  2:96553/1‑98 (MQ=255)                                                                                cTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACAt                                                                                       <  1:96553/98‑1 (MQ=255)                                                                                cTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCttt                                              >  1:1201207/1‑139 (MQ=255)                                                                                 ttGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTc                                             >  1:1333950/1‑139 (MQ=255)                                                                                   gCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                     <  2:1353556/65‑1 (MQ=255)                                                                                   gCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGca                                                                                                                     >  1:1353556/1‑65 (MQ=255)                                                                                    ccGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCTCATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAg                                          <  2:885613/139‑1 (MQ=255)                                                                                    ccGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAg                                          <  2:352131/139‑1 (MQ=255)                                                                                      gcgATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCa                                        <  2:171511/139‑1 (MQ=255)                                                                                      gcgATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCa                                        <  2:1202068/139‑1 (MQ=255)                                                                                       cgATAGCTATATCGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCaa                                       <  2:662147/139‑1 (MQ=255)                                                                                       cgATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCaa                                       <  1:1194934/139‑1 (MQ=255)                                                                                        gATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAAc                                      <  2:814977/139‑1 (MQ=255)                                                                                         aTAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACa                                     <  1:881180/139‑1 (MQ=255)                                                                                                      ttGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCa                                        >  1:745157/1‑123 (MQ=255)                                                                                                      ttGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCa                                        <  2:745157/123‑1 (MQ=255)                                                                                                      ttGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGAt                           >  1:548100/1‑136 (MQ=255)                                                                                                      ttGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGAt                           <  2:548100/136‑1 (MQ=255)                                                                                                       tGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACa                                     <  2:956398/125‑1 (MQ=255)                                                                                                       tGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACa                                     >  1:956398/1‑125 (MQ=255)                                                                                                       tGTCTGGCAAAGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAgc                       <  2:252193/139‑1 (MQ=255)                                                                                                                 aGGCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCcagcagca             <  1:619266/139‑1 (MQ=255)                                                                                                                  ggCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCcagcagcag            <  2:57252/139‑1 (MQ=255)                                                                                                                  ggCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCcagcagcag            <  2:1080313/139‑1 (MQ=255)                                                                                                                   gCATACGCGTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGt           <  1:1141863/139‑1 (MQ=255)                                                                                                                         gcgTTTTTTTCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGATa     <  1:258806/139‑1 (MQ=255)                                                                                                                            tttttttCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGATAAAg  <  1:1222328/139‑1 (MQ=255)                                                                                                                            tttttttCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTCGGa       >  1:896649/1‑134 (MQ=255)                                                                                                                            tttttttCATTAACCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAAAAGATCAGCGCCAGCAGCAGTGGGa       <  2:896649/134‑1 (MQ=255)                                                                                                                                        |                                                                                                                               AGAGTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAGCCGTAAATGCACATACTAAAGAACTTAACTAAACTTCACATCGCCGCTTCATTTTTCGCCCGTAGCTTCTTCGCCTGGCTTTCCAGCAACAGGTAACAGATCAGCGCCAGCAGCAGTGGGATAAAG  >  NC_000913/821371‑821632

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: ‑