Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,255,020	G→A	*234* (TAG→TAA)	yhaK →	redox‑sensitive bicupin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,255,020	0	G	A	100.0%	193.6 / NA	60	*234* (TAG→TAA)	yhaK	redox‑sensitive bicupin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/31);  total (29/31)

TCATGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAGTTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGAAG  >  NC_000913/3254888‑3255154                                                                                                                                     |                                                                                                                                       tCATGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGttt                                                                                                                                  <  2:180776/139‑1 (MQ=255)           aCGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGgag                                                                                                                        <  2:1540566/139‑1 (MQ=255)                   tGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                                >  2:545687/1‑139 (MQ=255)                                  gcTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGc                                                                                                 >  1:1342437/1‑139 (MQ=255)                                                aTGGGGCGTTTATTCGTGGCGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGc                                                                                   >  1:926796/1‑139 (MQ=255)                                                aTGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGc                                                                                   <  2:926796/139‑1 (MQ=255)                                                aTGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGc                                                                                   <  2:1407036/139‑1 (MQ=255)                                                aTGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGc                                                                                   <  2:1342437/139‑1 (MQ=255)                                                aTGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCACACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGc                                                                                   >  1:1311782/1‑139 (MQ=255)                                                  ggggCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCg                                                                                 >  2:1358047/1‑139 (MQ=255)                                                     gCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGt                                                                               >  2:471162/1‑138 (MQ=255)                                                     gCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGt                                                                               <  1:471162/138‑1 (MQ=255)                                                     gCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTg                                                                              >  2:1558609/1‑139 (MQ=255)                                                        tttATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                             <  1:683036/89‑1 (MQ=255)                                                        tttATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                             >  2:683036/1‑89 (MQ=255)                                                        tttATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                            >  1:1100630/1‑90 (MQ=255)                                                        tttATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                            <  2:1100630/90‑1 (MQ=255)                                                        tttATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGcc                                                                                      <  1:515630/128‑1 (MQ=255)                                                        tttATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGcc                                                                                      >  2:515630/1‑128 (MQ=255)                                                         ttATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGc                                                                          <  1:625315/139‑1 (MQ=255)                                                             tCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCgt                                                                      <  1:1558609/139‑1 (MQ=255)                                                              cGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCgtg                                                                     <  1:40107/139‑1 (MQ=255)                                                                      ggCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTa                                                                                                                               <  2:335036/73‑1 (MQ=255)                                                                      ggCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTa                                                                                                                               >  1:335036/1‑73 (MQ=255)                                                                                 cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGTCTGTCTAATTGTTTTTa                                                                                                                               >  2:576908/1‑62 (MQ=255)                                                                                 cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGTCTGTCTAATTGTTTTTa                                                                                                                               <  1:576908/62‑1 (MQ=255)                                                                                 cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTa                                                                                                                               <  1:1265981/62‑1 (MQ=255)                                                                                 cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTa                                                                                                                               >  2:1265981/1‑62 (MQ=255)                                                                                 cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                             <  1:1016926/64‑1 (MQ=255)                                                                                 cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                             >  2:1016926/1‑64 (MQ=255)                                                                                 cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTc                                                  >  1:442231/1‑139 (MQ=255)                                                                                      gTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGaa                                                         <  1:1608353/127‑1 (MQ=255)                                                                                      gTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGaa                                                         >  2:1608353/1‑127 (MQ=255)                                                                                       ttGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGatgatgatg                                                                                                               <  1:809990/72‑1 (MQ=255)                                                                                       ttGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGatgatgatg                                                                                                               >  2:809990/1‑72 (MQ=255)                                                                                       ttGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGaa                                                         <  2:1588443/126‑1 (MQ=255)                                                                                       ttGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGaa                                                         >  1:1588443/1‑126 (MQ=255)                                                                                          ccGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCa                                         <  1:545687/139‑1 (MQ=255)                                                                                                ccccACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGttt                                                                                                                                  <  1:1135386/44‑1 (MQ=255)                                                                                                ccccACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAAtt                                   >  1:934404/1‑139 (MQ=255)                                                                                                cccaACTGCGCGCTTTGCTGGTAGATTTGCCTGTCTAATTGttt                                                                                                                                  >  2:1135386/1‑44 (MQ=255)                                                                                                    aCTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTgg                               <  1:1358047/139‑1 (MQ=255)                                                                                                      tGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                            <  1:327772/44‑1 (MQ=255)                                                                                                      tGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                            >  2:327772/1‑44 (MQ=255)                                                                                                        cgcgcTTTGCTGGTAGATTTGGCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTg                                                                              <  1:408811/88‑1 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                                <  1:296051/54‑1 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                                >  2:296051/1‑54 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTg                                                                              >  1:365938/1‑88 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTg                                                                              <  2:365938/88‑1 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTg                                                                              >  2:408811/1‑88 (MQ=255)                                                                                                          cgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatg                                                                                                                  <  2:1142626/50‑1 (MQ=255)                                                                                                          cgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatg                                                                                                                  >  1:1142626/1‑50 (MQ=255)                                                                                                            cTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGcc                                                                                      >  2:193381/1‑76 (MQ=255)                                                                                                            cTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGcc                                                                                      <  1:193381/76‑1 (MQ=255)                                                                                                                       gATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTgagc            <  2:442231/139‑1 (MQ=255)                                                                                                                                gTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGaa   >  1:502242/1‑139 (MQ=255)                                                                                                                                 tCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGAAg  <  2:1311782/139‑1 (MQ=255)                                                                                                                                  cTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGcc                                                                                      >  1:1081585/1‑54 (MQ=255)                                                                                                                                  cTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCa                                               >  2:316057/1‑93 (MQ=255)                                                                                                                                  cTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCa                                               <  1:316057/93‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       TCATGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAGTTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGAAG  >  NC_000913/3254888‑3255154

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