Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	486,641	T→C	S36P (TCA→CCA)	mscK →	mechanosensitive channel protein, intermediate conductance, K+ regulated

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	486,641	0	T	C	100.0%	184.5 / NA	53	S36P (TCA→CCA)	mscK	mechanosensitive channel protein, intermediate conductance, K+ regulated
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/23);  total (30/23)

TTCTTCAGGTTCAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCA  >  NC_000913/486507‑486764                                                                                                                                       |                                                                                                                            ttcttcAGGTTCAGAAACCTTCATTCATCATGACTATGTTCAAATATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTACTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATc                                                                                                                         >  2:1192195/1‑139 (MQ=255)            cAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCa                                                                                                                           >  1:376163/1‑126 (MQ=255)            cAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCa                                                                                                                           <  2:376163/126‑1 (MQ=255)              gAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTg                                                                                                            >  1:1018524/1‑139 (MQ=255)                       attcatCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACaaaa                                                                                                   >  2:1150371/1‑139 (MQ=255)                                 aCTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGc                                                                                         >  1:517000/1‑139 (MQ=255)                                   tATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAg                                                                                       >  1:26894/1‑139 (MQ=255)                                    aTGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAgg                                                                                      <  2:517000/139‑1 (MQ=255)                                    aTGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAgg                                                                                      <  2:26894/139‑1 (MQ=255)                                     tGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCa                                                                                                                           >  1:333521/1‑101 (MQ=255)                                     tGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGgc                                                                                     >  2:26326/1‑139 (MQ=255)                                             tATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTg                                                                             >  1:374106/1‑139 (MQ=255)                                               ttACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGa                                                                                              >  2:405835/1‑120 (MQ=255)                                               ttACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGa                                                                                              <  1:405835/120‑1 (MQ=255)                                                 aCAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                                         <  1:26326/139‑1 (MQ=255)                                                 aCAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                                         >  1:417860/1‑139 (MQ=255)                                                 aCAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                                         >  1:483066/1‑139 (MQ=255)                                                 aCAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                                         >  1:495152/1‑139 (MQ=255)                                                            cGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAAc                                                              <  1:563695/139‑1 (MQ=255)                                                                 ttttGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaaaaa                                                         >  2:831772/1‑139 (MQ=255)                                                                                aTTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATc                                                                                                                         <  1:1114966/60‑1 (MQ=255)                                                                                aTTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATc                                                                                                                         >  2:1114966/1‑60 (MQ=255)                                                                                aTTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCAACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAg                                          >  2:22793/1‑139 (MQ=255)                                                                                       gCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAAc                                   <  1:1150371/139‑1 (MQ=255)                                                                                         ttttgtttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCa                                                                                  <  2:488786/90‑1 (MQ=255)                                                                                         ttttgtttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCa                                                                                  >  1:488786/1‑90 (MQ=255)                                                                                          tttgtttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGa                                                                            <  2:1053229/95‑1 (MQ=255)                                                                                          tttgtttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGa                                                                            >  1:1053229/1‑95 (MQ=255)                                                                                          tttgtttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTgg                                <  1:831772/139‑1 (MQ=255)                                                                                            tgtttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTg                              <  1:678470/139‑1 (MQ=255)                                                                                                  tgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGa                                                                            <  2:566642/87‑1 (MQ=255)                                                                                                  tgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGa                                                                            >  1:566642/1‑87 (MQ=255)                                                                                                       tGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                                         <  2:667470/85‑1 (MQ=255)                                                                                                       tGTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                                         >  1:667470/1‑85 (MQ=255)                                                                                                        gTTATGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTg                  <  2:483066/139‑1 (MQ=255)                                                                                                         gtATGACAGAACACTGGGTTTGCGCGGGCGCCACCGAATGGTGATCTGACGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAGACAAAAAGATCTTTCTGCTCAGTACGAACTGGTGCACCAGGATCTGa                 <  2:860914/138‑1 (MQ=255)                                                                                                            tGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAg              <  1:1192195/139‑1 (MQ=255)                                                                                                            tGCCAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCAACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAg              <  1:22793/139‑1 (MQ=255)                                                                                                              acAGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCGGCCGACAAAAGCGGACCTGCAGGCGCAACTTAACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGAAGCAGGATCTGACAGAt            >  1:794879/2‑139 (MQ=255)                                                                                                                aGAACACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATAc          >  2:792925/1‑139 (MQ=255)                                                                                                                  aaCACGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACAt        >  2:635503/1‑139 (MQ=255)                                                                                                                   acacGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAg              >  1:356692/1‑132 (MQ=255)                                                                                                                   acacGGCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAg              <  2:356692/132‑1 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTg                                                                                          <  1:1071950/52‑1 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTg                                                                                          >  2:1071950/1‑52 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCt                                             <  2:4062/97‑1 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCt                                             >  1:4062/1‑97 (MQ=255)                                                                                                                       ggCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGcc   >  1:234216/1‑139 (MQ=255)                                                                                                                        gCGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCa  >  1:926610/1‑139 (MQ=255)                                                                                                                         cGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATc                    >  2:431640/1‑120 (MQ=255)                                                                                                                         cGTTTGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATc                    <  1:431640/120‑1 (MQ=255)                                                                                                                            ttGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAg              >  1:1164099/1‑123 (MQ=255)                                                                                                                            ttGCGCGGGCGCCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAg              <  2:1164099/123‑1 (MQ=255)                                                                                                                                       |                                                                                                                            TTCTTCAGGTTCAGAAACCTTCATTCATCATGACTATGTTCCAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCA  >  NC_000913/486507‑486764

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