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

CACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGTTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCA  >  NC_000913/586148‑586415                                                                                                                                     |                                                                                                                                        caccacATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   <  2:1539424/139‑1 (MQ=255)      caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATc                                                                                                                                     <  1:1595100/132‑1 (MQ=255)      caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATc                                                                                                                                     >  2:1595100/1‑132 (MQ=255)      caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              >  1:856430/1‑139 (MQ=255)           cTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATACTTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         <  2:771246/139‑1 (MQ=255)           cTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTAt                                                                                                                                      >  1:1011031/1‑126 (MQ=255)           cTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTAt                                                                                                                                      <  2:1011031/126‑1 (MQ=255)           cTTTTTGTCGCTGTTCGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         >  1:771246/1‑139 (MQ=255)                tGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    >  1:1464689/1‑139 (MQ=255)                       gTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                             <  2:439696/139‑1 (MQ=255)                        tttGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTa                                                                                                                                       >  2:680376/1‑112 (MQ=255)                        tttGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    <  1:541159/131‑1 (MQ=255)                        tttGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    >  2:541159/1‑131 (MQ=255)                                          cGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              >  1:1542961/1‑103 (MQ=255)                                          cGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              <  2:1542961/103‑1 (MQ=255)                                               tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     <  2:1464689/139‑1 (MQ=255)                                                    ccGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGc                                                                                <  2:1359315/139‑1 (MQ=255)                                                       tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCtat                                                                             >  2:1221556/1‑139 (MQ=255)                                                         aaaCGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              >  1:233981/1‑88 (MQ=255)                                                         aaaCGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              <  2:233981/88‑1 (MQ=255)                                                                            aTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTaa                                                                                                            <  2:1590554/87‑1 (MQ=255)                                                                            aTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTaa                                                                                                            >  1:1590554/1‑87 (MQ=255)                                                                                    cgacgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGaca                                                                                                 >  1:726046/1‑90 (MQ=255)                                                                                    cgacgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGaca                                                                                                 <  2:726046/90‑1 (MQ=255)                                                                                     gacgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCtttt                                               >  2:856769/1‑139 (MQ=255)                                                                                      acgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     <  2:249138/100‑1 (MQ=255)                                                                                      acgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     >  1:249138/1‑100 (MQ=255)                                                                                           ctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                             >  2:1513074/1‑71 (MQ=255)                                                                                           ctgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                             <  1:1513074/71‑1 (MQ=255)                                                                                            tgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGaca                                                                                                 >  1:1452281/1‑82 (MQ=255)                                                                                            tgctATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGaca                                                                                                 <  2:1452281/82‑1 (MQ=255)                                                                                               tataGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCt                                     <  2:856430/139‑1 (MQ=255)                                                                                                      acaTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    <  1:965346/53‑1 (MQ=255)                                                                                                      acaTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    >  2:965346/1‑53 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         <  2:418586/46‑1 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         <  1:511442/46‑1 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         >  2:511442/1‑46 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACAt                                                                                                                         >  1:418586/1‑46 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   <  2:268881/68‑1 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   >  1:268881/1‑68 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTaa                            >  1:1409236/1‑139 (MQ=255)                                                                                                        aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTaa                            >  1:1286001/1‑139 (MQ=255)                                                                                                           aTTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACa                         <  1:1221556/139‑1 (MQ=255)                                                                                                              gcgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAg                                           >  2:678617/1‑118 (MQ=255)                                                                                                              gcgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAg                                           <  1:678617/118‑1 (MQ=255)                                                                                                               cgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     <  2:79718/75‑1 (MQ=255)                                                                                                               cgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     >  1:79718/1‑75 (MQ=255)                                                                                                               cgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGc                     >  2:1321454/1‑139 (MQ=255)                                                                                                                 ccACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGTAGGCTgg                                   >  1:284907/1‑123 (MQ=255)                                                                                                                 ccACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGTAGGCTgg                                   <  2:284907/123‑1 (MQ=255)                                                                                                                 ccACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGa                   >  1:363945/1‑139 (MQ=255)                                                                                                                  cACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGaa                  >  2:128639/1‑139 (MQ=255)                                                                                                                    ccTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGa                   <  2:166893/136‑1 (MQ=255)                                                                                                                    ccTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGa                   >  1:166893/1‑136 (MQ=255)                                                                                                                    ccTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACt                >  1:1296090/1‑139 (MQ=255)                                                                                                                      tGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGTCAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTgg              >  2:1102051/1‑139 (MQ=255)                                                                                                                            atatTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     >  2:451095/1‑62 (MQ=255)                                                                                                                            atatTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     <  1:451095/62‑1 (MQ=255)                                                                                                                              atTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   >  2:1296366/1‑46 (MQ=255)                                                                                                                              atTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAAtt      >  1:865215/1‑139 (MQ=255)                                                                                                                              atTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAAtt      >  2:264968/1‑139 (MQ=255)                                                                                                                              atTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAAtt      >  2:123381/1‑139 (MQ=255)                                                                                                                              atTTTTGCGATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   <  1:1296366/46‑1 (MQ=255)                                                                                                                                ttttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCg    <  2:1296090/139‑1 (MQ=255)                                                                                                                                  ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGc                                                                          >  1:420362/1‑67 (MQ=255)                                                                                                                                  ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGc                                                                          <  2:420362/67‑1 (MQ=255)                                                                                                                                  ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCa  <  2:363945/139‑1 (MQ=255)                                                                                                                                  ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCa  <  1:1321454/139‑1 (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 < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑