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%	182.6 / NA	56	*234* (TAG→TAA)	yhaK	redox‑sensitive bicupin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/26);  total (30/26)

TCATGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAGTTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGAAG  >  NC_000913/3254888‑3255154                                                                                                                                     |                                                                                                                                       tCATGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGttt                                                                                                                                  <  1:304893/139‑1 (MQ=255)      cgcTTACGCATCATGAAGATAAAG‑GGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                            <  1:842359/139‑1 (MQ=255)        gttACGCAACATGAAGAGAAAGCGGCGCTGACCTGCGGTGTTGGGGCGTTTATTCGTGACGAGGTTAACATTACGATGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTATTACAgg                                                                                                                           <  1:1180750/138‑1 (MQ=255)        cTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                             <  1:1430021/137‑1 (MQ=255)        cTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                             >  2:1430021/1‑137 (MQ=255)          tACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAggg                                                                                                                          >  2:845049/1‑138 (MQ=255)          tACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAggg                                                                                                                          <  1:845049/138‑1 (MQ=255)                 catGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatg                                                                                                                  >  2:1293976/1‑139 (MQ=255)                            aGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTaaaaa                                                                                                       >  1:774953/1‑139 (MQ=255)                                           cGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                            >  1:1406508/1‑103 (MQ=255)                                           cGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                            <  2:1406508/103‑1 (MQ=255)                                                     gCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAgg                                                                                                                           >  2:818816/1‑94 (MQ=255)                                                     gCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAgg                                                                                                                           <  1:818816/94‑1 (MQ=255)                                                     gCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTg                                                                              >  1:1319586/1‑139 (MQ=255)                                                     gCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTg                                                                              >  2:1173602/1‑139 (MQ=255)                                                          tATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGcc                                                                         >  1:1485335/1‑139 (MQ=255)                                                                      ggCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGaa                                                             >  2:681872/1‑139 (MQ=255)                                                                                 cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                                >  2:546544/1‑77 (MQ=255)                                                                                 cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                                <  1:546544/77‑1 (MQ=255)                                                                                  gCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGgag                                                                                                                        <  1:39585/68‑1 (MQ=255)                                                                                  gCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGgag                                                                                                                        >  2:39585/1‑68 (MQ=255)                                                                                      gttaccGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                            <  2:849428/56‑1 (MQ=255)                                                                                      gttaccGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                            >  1:849428/5‑60 (MQ=255)                                                                                      gTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatg                                                                                                                  >  1:1245027/1‑70 (MQ=255)                                                                                      gTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatg                                                                                                                  <  2:1245027/70‑1 (MQ=255)                                                                                              ttCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAggg                                                                                                                          <  2:1597520/54‑1 (MQ=255)                                                                                              ttCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAggg                                                                                                                          >  1:1597520/1‑54 (MQ=255)                                                                                                   cACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTg                                <  1:1173602/139‑1 (MQ=255)                                                                                                    aCTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAgcg                                                                                        <  2:87712/82‑1 (MQ=255)                                                                                                    aCTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAgcg                                                                                        >  1:87712/1‑82 (MQ=255)                                                                                                    aCTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTgg                               >  2:549276/1‑139 (MQ=255)                                                                                                     cTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCaaa                                                                                             >  1:845026/1‑76 (MQ=255)                                                                                                     cTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCaaa                                                                                             <  2:845026/76‑1 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGgaga                                                                                                                       >  2:1390430/1‑47 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGgaga                                                                                                                       <  1:1390430/47‑1 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAgcg                                                                                        <  2:574407/78‑1 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAgcg                                                                                        >  1:574407/1‑78 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCg                                                                                 <  1:1013212/85‑1 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCg                                                                                 >  2:1013212/1‑85 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTAt                           >  2:651305/1‑139 (MQ=255)                                                                                                        cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTAt                           >  1:468643/1‑139 (MQ=255)                                                                                                          cgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTa                                                                                                                               <  2:488402/37‑1 (MQ=255)                                                                                                          cgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTa                                                                                                                               >  1:488402/1‑37 (MQ=255)                                                                                                          cgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGa                         >  2:632928/1‑139 (MQ=255)                                                                                                           gcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGaa                        <  1:681872/139‑1 (MQ=255)                                                                                                           gcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGaa                        <  2:1319586/139‑1 (MQ=255)                                                                                                            cTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAag                       <  1:1293976/139‑1 (MQ=255)                                                                                                                gCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGAt                   >  1:1234652/1‑139 (MQ=255)                                                                                                                     tAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTga              >  2:523385/1‑139 (MQ=255)                                                                                                                      aGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTgag             <  2:774953/139‑1 (MQ=255)                                                                                                                         tttGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGaa                        <  1:612529/125‑1 (MQ=255)                                                                                                                         tttGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGaa                        >  2:612529/1‑125 (MQ=255)                                                                                                                                 tCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCTGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTgg                               >  1:1415494/1‑110 (MQ=255)                                                                                                                                 tCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCTGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTgg                               <  2:1415494/110‑1 (MQ=255)                                                                                                                                 tCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGAAg  <  2:1485335/139‑1 (MQ=255)                                                                                                                                  cTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTgg                               >  2:512673/1‑109 (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 < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑