Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	586,280	T→C	N210S (AAC→AGC)	envY ←	porin thermoregulatory transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	586,280	0	T	C	100.0%	231.5 / NA	64	N210S (AAC→AGC)	envY	porin thermoregulatory transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/34);  total (30/34)

CTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGTTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTC  >  NC_000913/586145‑586412                                                                                                                                        |                                                                                                                                     cTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTAt                                                                                                                                   <  1:1254273/139‑1 (MQ=255) cTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTAt                                                                                                                                   <  2:343630/139‑1 (MQ=255)    caccacATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                <  1:1061561/139‑1 (MQ=255)     accacATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCAt                                                                                                                               <  2:745453/139‑1 (MQ=255)      ccacATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATc                                                                                                                              <  1:422369/139‑1 (MQ=255)         caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           <  2:1073467/139‑1 (MQ=255)            cACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTACTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAAc                                                                                                                        <  2:346765/139‑1 (MQ=255)                  ttGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGt                                                                                                                  <  2:1268053/139‑1 (MQ=255)                   tGTCGCTGTTTGGCGAGATAGTTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                 >  1:946203/1‑139 (MQ=255)                            ttGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                        >  1:918341/1‑139 (MQ=255)                                   gaTAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACg                                                                                                 >  1:1425175/1‑139 (MQ=255)                                     tAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                          >  2:668869/1‑128 (MQ=255)                                     tAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                          <  1:668869/128‑1 (MQ=255)                                        ttCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACAct                                                                                            <  2:1425175/139‑1 (MQ=255)                                             cGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCtgtg                                                                                       >  1:520175/1‑139 (MQ=255)                                              ggTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                >  2:1197363/1‑129 (MQ=255)                                              ggTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                <  1:1197363/129‑1 (MQ=255)                                                tgtCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAAcggatagc                                                                                                                 >  1:1595670/1‑103 (MQ=255)                                                tgtCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAAcggatagc                                                                                                                 <  2:1595670/110‑8 (MQ=255)                                                tgtCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACa                                                                                    >  1:355946/1‑139 (MQ=255)                                                  tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           <  1:559436/98‑1 (MQ=255)                                                  tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           >  2:559436/1‑98 (MQ=255)                                                  tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                  <  2:918341/139‑1 (MQ=255)                                                            aaaCGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAg                                                                        >  1:172888/1‑139 (MQ=255)                                                                                  aaGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAAtt                                                  <  1:848995/139‑1 (MQ=255)                                                                                       cgacgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCttt                                             >  2:850756/1‑139 (MQ=255)                                                                                         acgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                  >  1:927983/1‑100 (MQ=255)                                                                                         acgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                  <  2:927983/100‑1 (MQ=255)                                                                                         acgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCg                                                            <  2:1318701/122‑1 (MQ=255)                                                                                         acgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCg                                                            >  1:1318701/1‑122 (MQ=255)                                                                                           gTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGa                                         <  2:172888/139‑1 (MQ=255)                                                                                            tgctgcTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                <  1:952947/51‑1 (MQ=255)                                                                                            tgctgcTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                >  2:952947/1‑51 (MQ=255)                                                                                             gctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATc                                                                                                                              <  1:145717/52‑1 (MQ=255)                                                                                             gctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATc                                                                                                                              >  2:145717/1‑52 (MQ=255)                                                                                              ctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                        <  1:1467660/73‑1 (MQ=255)                                                                                              ctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                        >  2:1467660/1‑73 (MQ=255)                                                                                               tgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGt                                                                                                           <  2:506966/69‑1 (MQ=255)                                                                                               tgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGt                                                                                                           >  1:506966/1‑69 (MQ=255)                                                                                                      gCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTgggg                              >  2:635073/1‑139 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           >  2:15402/1‑41 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           <  1:15402/41‑1 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           >  1:1448371/1‑41 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           <  2:1448371/41‑1 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           <  2:1372993/41‑1 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           >  1:1372993/1‑41 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           >  2:920803/1‑41 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                           <  1:920803/41‑1 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                        <  1:575242/60‑1 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                        >  2:575242/1‑60 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGAc                                                                                     >  2:204523/1‑79 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGAc                                                                                     <  1:204523/79‑1 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTaa                         >  2:393207/1‑139 (MQ=255)                                                                                                                   gcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                  <  2:945153/74‑1 (MQ=255)                                                                                                                   gcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                  >  1:945153/1‑74 (MQ=255)                                                                                                                     cACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                  <  1:1107457/72‑1 (MQ=255)                                                                                                                     cACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                  >  2:1107457/1‑72 (MQ=255)                                                                                                                       ccTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                        <  1:1294019/48‑1 (MQ=255)                                                                                                                       ccTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                        >  2:1294019/1‑48 (MQ=255)                                                                                                                           aGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                        <  2:83188/44‑1 (MQ=255)                                                                                                                           aGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                        >  1:83188/1‑44 (MQ=255)                                                                                                                           aGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCg         <  1:393207/139‑1 (MQ=255)                                                                                                                            gTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGa        <  1:635073/139‑1 (MQ=255)                                                                                                                                  tttttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTc  >  1:1254059/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                     CTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGTTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTC  >  NC_000913/586145‑586412

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