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

CTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGTTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCA  >  NC_000913/586145‑586415                                                                                                                                        |                                                                                                                                        cTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTAt                                                                                                                                      <  1:1044246/139‑1 (MQ=255)  ttCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATc                                                                                                                                     >  1:567104/1‑139 (MQ=255)   tCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATcc                                                                                                                                    <  1:11991/139‑1 (MQ=255)    caccacATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   <  1:541731/139‑1 (MQ=255)         caTCCCTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              <  1:578434/139‑1 (MQ=255)         caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGTCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATcaatcaat                                                                                                                              <  2:189746/139‑1 (MQ=255)         caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              <  1:1007677/139‑1 (MQ=255)            cACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTAt                                                                                                                                      >  2:464880/1‑128 (MQ=255)            cACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTAt                                                                                                                                      <  1:464880/128‑1 (MQ=255)                   tGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    <  2:1137717/139‑1 (MQ=255)                                    aTAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                             >  1:1021375/1‑129 (MQ=255)                                    aTAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                             <  2:1021375/129‑1 (MQ=255)                                    aTAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   <  1:356315/139‑1 (MQ=255)                                        ttCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACAct                                                                                               >  2:842569/1‑139 (MQ=255)                                        ttCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACAct                                                                                               <  2:567104/139‑1 (MQ=255)                                                  tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         <  2:268226/103‑1 (MQ=255)                                                  tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         >  1:268226/1‑103 (MQ=255)                                                  tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    >  1:276062/1‑108 (MQ=255)                                                  tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    <  2:276062/108‑1 (MQ=255)                                                        cGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCt                                                                               <  1:842569/139‑1 (MQ=255)                                                          tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATcc                                                                                                                                    >  2:625965/1‑84 (MQ=255)                                                          tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATcc                                                                                                                                    <  1:625965/84‑1 (MQ=255)                                                                             agaTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTg                                                                        >  1:1088709/1‑125 (MQ=255)                                                                             agaTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTg                                                                        <  2:1088709/125‑1 (MQ=255)                                                                               aTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTa                                                        >  1:283750/1‑139 (MQ=255)                                                                                  aaGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAAtt                                                     <  2:283750/139‑1 (MQ=255)                                                                                  aaGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAAtt                                                     >  1:1202417/1‑139 (MQ=255)                                                                                       cgatgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   <  1:42541/52‑1 (MQ=38)                                                                                       cgatgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   >  2:42541/5‑56 (MQ=255)                                                                                        gacgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   >  2:1152215/1‑55 (MQ=255)                                                                                        gacgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   <  1:1152215/55‑1 (MQ=255)                                                                                            tgctgcTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAg                                           <  2:864345/139‑1 (MQ=255)                                                                                             gctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGc                                          >  1:541026/1‑139 (MQ=255)                                                                                              ctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAg                                           <  1:448996/137‑1 (MQ=255)                                                                                              ctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAg                                           >  2:448996/1‑137 (MQ=255)                                                                                                  tataGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCt                                     <  2:1202417/139‑1 (MQ=255)                                                                                                    taGCCACATAATTGGGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATAGGACACTCTGTGACAATCTGGCTATAGGTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTgg                                   >  1:785311/1‑139 (MQ=255)                                                                                                         acaTAATTGCGCCACCTGAGTGATCTTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCtt                              <  1:721378/139‑1 (MQ=255)                                                                                                         acaTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCt                                     <  2:316548/132‑1 (MQ=255)                                                                                                         acaTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCt                                     >  1:316548/1‑132 (MQ=255)                                                                                                         acaTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCtt                              >  1:337449/1‑139 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATcc                                                                                                                                    <  1:774762/35‑1 (MQ=255)                                                                                                           aTAATTGCGCCACCTGAGTGATATTTTTGCTATcc                                                                                                                                    >  2:774762/1‑35 (MQ=255)                                                                                                              aTTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACa                         <  2:541026/139‑1 (MQ=255)                                                                                                                  cgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    <  2:1291814/44‑1 (MQ=255)                                                                                                                  cgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    >  1:1291814/1‑44 (MQ=255)                                                                                                                       ccTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACAct                                                                                               >  1:954949/1‑60 (MQ=255)                                                                                                                       ccTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACAct                                                                                               <  2:954949/60‑1 (MQ=255)                                                                                                                          gAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   >  2:462761/1‑53 (MQ=255)                                                                                                                          gAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGc             >  1:583877/1‑139 (MQ=255)                                                                                                                          gAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGGAACGCATACGa                                                                                                   <  1:462761/53‑1 (MQ=255)                                                                                                                                   ttttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCg    <  2:689439/139‑1 (MQ=255)                                                                                                                                     ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCa  <  2:482060/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        CTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGTTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCA  >  NC_000913/586145‑586415

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