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

GCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGTTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCA  >  NC_000913/586142‑586415                                                                                                                                           |                                                                                                                                        caccTTCACCACATCACAGTTTGTCGGACGTTGGCGAGCTAATACAACGGTGTCAGGCCGTAAACCGCATCAACAACAGTTACGAAATACGACGGGCTGCAATAACCACATAATTGCGACACCTGAATGATATGTCTGt                                                                                                                                         >  1:247927/3‑139 (MQ=255)      tCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATcc                                                                                                                                    >  1:192808/1‑139 (MQ=255)            caTCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAAt                                                                                                                              >  1:722353/1‑139 (MQ=255)                 cTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   >  2:1199941/1‑129 (MQ=255)                 cTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCa                                                                                                                                   <  1:1199941/129‑1 (MQ=255)                     ttGTCGCTGTTTGGCGAGATAATTCATCGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGt                                                                                                                     <  1:361542/139‑1 (MQ=255)                      tGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTa                                                                                                                    <  2:192808/139‑1 (MQ=255)                         cGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACgg                                                                                                                 >  1:475251/1‑139 (MQ=255)                                      gaTAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                             >  2:496886/1‑130 (MQ=255)                                      gaTAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTa                                                                                                             <  1:496886/130‑1 (MQ=255)                                                     tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     <  1:983397/139‑1 (MQ=255)                                                     tCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     <  2:580628/139‑1 (MQ=255)                                                             tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCAt                                                                                                                                  >  1:442508/1‑86 (MQ=255)                                                             tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCAt                                                                                                                                  <  2:442508/86‑1 (MQ=255)                                                             tAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCtat                                                                             >  2:663379/1‑139 (MQ=255)                                                                               gagaTAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGtttt                                                           >  2:1218292/1‑139 (MQ=255)                                                                                          cgacgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTCTTTTAATTTCttt                                                >  2:480449/1‑139 (MQ=255)                                                                                            acgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                           >  1:482552/1‑78 (MQ=255)                                                                                            acgTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                           <  2:482552/78‑1 (MQ=255)                                                                                            acgTGCTGATATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCttttt                                              >  2:793135/1‑139 (MQ=255)                                                                                               tgctgcTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAg                                           <  2:722353/139‑1 (MQ=255)                                                                                                         gCCACATAATTGCGCCACCTGAGTGATATTTTTGCTa                                                                                                                                       >  1:509445/1‑37 (MQ=255)                                                                                                         gCCACATAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTgggg                                 >  2:562442/1‑139 (MQ=255)                                                                                                            acaTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAg                                           <  1:648685/126‑1 (MQ=255)                                                                                                            acaTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCtt                              >  1:591209/1‑139 (MQ=255)                                                                                                            acaTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAGCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAg                                           >  2:648685/1‑126 (MQ=255)                                                                                                              aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                           <  1:835962/60‑1 (MQ=255)                                                                                                              aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAAc                                                                                                           >  2:835962/1‑60 (MQ=255)                                                                                                              aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCt                                     <  1:76815/130‑1 (MQ=255)                                                                                                              aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCt                                     >  2:1156547/1‑130 (MQ=255)                                                                                                              aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCt                                     >  2:76815/1‑130 (MQ=255)                                                                                                              aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCt                                     <  1:1156547/130‑1 (MQ=255)                                                                                                              aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTgggg                                 >  2:1029025/1‑134 (MQ=255)                                                                                                              aTAATTGCGCCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTgggg                                 <  1:1029025/134‑1 (MQ=255)                                                                                                                    gcgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGAc                                                                                        >  2:951647/1‑73 (MQ=255)                                                                                                                    gcgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGAc                                                                                        <  1:951647/73‑1 (MQ=255)                                                                                                                    gcgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAg                      >  2:612445/1‑139 (MQ=255)                                                                                                                    gcgcCACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAg                      >  1:28064/1‑139 (MQ=255)                                                                                                                        cACCTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATCTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGaa                  <  1:663379/139‑1 (MQ=255)                                                                                                                          ccTGAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACt                >  2:237558/1‑139 (MQ=255)                                                                                                                             gAGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGc             >  1:931586/1‑139 (MQ=255)                                                                                                                              aGTGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCg            <  1:562442/139‑1 (MQ=255)                                                                                                                                tGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   >  2:493761/1‑50 (MQ=255)                                                                                                                                tGATATTTTTGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGa                                                                                                   <  1:493761/50‑1 (MQ=255)                                                                                                                                        ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACTCTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCa  <  1:793135/139‑1 (MQ=255)                                                                                                                                        ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAAt                                                                                     <  2:194113/56‑1 (MQ=255)                                                                                                                                        ttGCTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCa  <  1:1218292/139‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        GCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCTTAAAAACAGAGATAAAGTACGACGTGCTGCTATAGCCACATAATTGCGCCACCTGAGTGATATTTTTGTTATCCATCAATAACATCTGTACGGCGTAACGCATACGACACTCTGTGACAATCTGGCTATAGCTGGTATTTTCGTTTTTTAATTTCTTTTTGAGCAGGCTGGGGCTTAAACATAGCGAACTGGCGACAATTCGCA  >  NC_000913/586142‑586415

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