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%	234.5 / NA	66	N210S (AAC→AGC)	envY	porin thermoregulatory transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/33);  total (33/33)

CACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGTTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCA  >  NC_000913/586148‑586415                                                                                                                                     |                                                                                                                                        caccacATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   <  2:926265/139‑1 (MQ=255) caccacATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   >  2:1790298/1‑139 (MQ=255) caccacATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCACCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   <  2:1165620/139‑1 (MQ=255)    cacaTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCa                                                                                                                                <  1:1973955/139‑1 (MQ=255)      caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              <  1:1475890/139‑1 (MQ=255)      caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              <  2:1860794/139‑1 (MQ=255)      caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              <  2:1798321/139‑1 (MQ=255)         cACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTa                                                                                                                                       >  1:813423/1‑127 (MQ=255)           cTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         <  2:1644370/139‑1 (MQ=255)           cTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         >  2:1734403/1‑139 (MQ=255)                       gTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                             <  1:90593/139‑1 (MQ=255)                          tGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACg                                                                                                          >  1:234044/1‑139 (MQ=255)                                  tAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATc                                                                                                                                     >  2:1731304/1‑104 (MQ=255)                                  tAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATc                                                                                                                                     <  1:1731304/104‑1 (MQ=255)                                   aaTTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGaca                                                                                                 <  1:382531/139‑1 (MQ=255)                                          cGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCtgtg                                                                                          <  2:1846200/139‑1 (MQ=255)                                               tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCTCCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     <  1:1843610/139‑1 (MQ=255)                                                       tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         <  1:494961/95‑1 (MQ=255)                                                       tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         >  2:494961/1‑95 (MQ=255)                                                       tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   <  1:429725/117‑1 (MQ=255)                                                       tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   >  2:429725/1‑117 (MQ=255)                                                       tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCtat                                                                             >  1:905358/1‑139 (MQ=255)                                                             gCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         <  2:1842287/89‑1 (MQ=255)                                                             gCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         >  1:1842287/1‑89 (MQ=255)                                                                  aaaaaCAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTa                                                                                                                                       >  1:1303715/1‑70 (MQ=255)                                                                  aaaaaCAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTa                                                                                                                                       <  2:1303715/70‑1 (MQ=255)                                                                          agaTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         <  2:587225/76‑1 (MQ=255)                                                                          agaTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         >  1:587225/1‑76 (MQ=255)                                                                            aTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCtgtg                                                                                          >  1:1855944/1‑105 (MQ=255)                                                                            aTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCtgtg                                                                                          <  2:1855944/105‑1 (MQ=255)                                                                                    cgacgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCttt                                                >  2:562551/1‑139 (MQ=255)                                                                                         tgctgcTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAg                                           <  1:1734403/139‑1 (MQ=255)                                                                                           ctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTa                                                                                                                                       <  1:189355/45‑1 (MQ=255)                                                                                           ctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              <  1:170604/54‑1 (MQ=255)                                                                                           ctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              >  2:170604/1‑54 (MQ=255)                                                                                            tgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATATGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAg                                        >  2:960142/1‑139 (MQ=255)                                                                                                   gCCACATAATTGCGCCACCTGAGTGATATTTTTGCTa                                                                                                                                       >  2:859431/1‑37 (MQ=255)                                                                                                   gCCACATAATTGCGCCACCTGAGTGATATTTTTGCTa                                                                                                                                       <  1:859431/37‑1 (MQ=255)                                                                                                      acaTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCtt                              >  2:1742395/1‑139 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              >  2:159531/1‑41 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              <  1:159531/41‑1 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         <  2:1671437/46‑1 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         >  1:1671437/1‑46 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGaca                                                                                                 <  2:1361232/70‑1 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGaca                                                                                                 >  1:1361232/1‑70 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCtgt                                                                                           >  1:2001862/1‑76 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCtgt                                                                                           <  2:2001862/76‑1 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTaa                            <  2:234044/139‑1 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGAGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAACAGGCAGGGGCTTaa                            >  1:334424/1‑139 (MQ=255)                                                                                                             tGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATa                       >  1:395661/1‑139 (MQ=255)                                                                                                              gcgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAg                      >  2:1868139/1‑139 (MQ=255)                                                                                                               cgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   >  2:1410344/1‑61 (MQ=255)                                                                                                               cgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   <  1:1410344/61‑1 (MQ=255)                                                                                                               cgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     >  2:1167658/1‑75 (MQ=255)                                                                                                               cgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     <  1:1167658/75‑1 (MQ=255)                                                                                                                 ccACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCtt                              <  1:431553/128‑1 (MQ=255)                                                                                                                 ccACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCtt                              >  2:431553/1‑128 (MQ=255)                                                                                                                 ccACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGa                   >  2:628586/1‑139 (MQ=255)                                                                                                                 ccACCTGAGTGATATTTTTGCTATCCATCAATAACACCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGa                   >  2:342742/1‑139 (MQ=255)                                                                                                                  cACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGaa                  >  1:1390798/1‑139 (MQ=255)                                                                                                                            atatTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                           >  1:53170/1‑40 (MQ=255)                                                                                                                            atatTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                           <  2:53170/40‑1 (MQ=255)                                                                                                                              atTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAAtt      >  1:179358/1‑139 (MQ=255)                                                                                                                                  ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   <  2:776349/42‑1 (MQ=255)                                                                                                                                  ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCa  <  1:1717374/139‑1 (MQ=255)                                                                                                                                  ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCa  >  2:1717374/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                        CACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGTTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCA  >  NC_000913/586148‑586415

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