Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,177,613	A→G	T86A (ACT→GCT)	secE →	Sec translocon subunit SecE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,177,613	0	A	G	100.0%	132.3 / NA	36	T86A (ACT→GCT)	secE	Sec translocon subunit SecE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/19);  total (17/19)

CCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGACTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAG  >  NC_000913/4177481‑4177734                                                                                                                                     |                                                                                                                          ccGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGa                                                                                                                           <  2:351541/133‑1 (MQ=255) ccGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGa                                                                                                                           >  1:351541/1‑133 (MQ=255)    cTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGa                                                                                                                 >  1:207624/1‑140 (MQ=255)    cTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGa                                                                                                                 >  1:295576/1‑140 (MQ=255)    cTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGa                                                                                                                 >  2:144832/1‑140 (MQ=255)      gcgtgcgCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGaaa                                                                                                               >  2:222257/1‑140 (MQ=255)      gcgtgcgCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGaaa                                                                                                               <  1:59343/140‑1 (MQ=255)       cgtgcgCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGc                                                                                                                      <  1:346428/132‑1 (MQ=255)       cgtgcgCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGc                                                                                                                      >  2:346428/1‑132 (MQ=255)              tGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGcac                                                                                                       <  2:323107/140‑1 (MQ=255)               ggCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGcaca                                                                                                      >  2:170770/1‑140 (MQ=255)                              gATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTgg                                                                                       <  1:115851/140‑1 (MQ=255)                               aTTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGc                                                                                      >  2:285525/1‑140 (MQ=255)                                 tgctgcAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCt                                                                                     <  2:314917/139‑1 (MQ=255)                                 tgctgcAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCt                                                                                     >  1:314917/1‑139 (MQ=255)                                        gCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTAc                                                                             <  1:144832/140‑1 (MQ=255)                                        gCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTAc                                                                             <  1:285525/140‑1 (MQ=255)                                            gtggtgTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCa                                                                         <  2:207624/140‑1 (MQ=255)                                                   cgcgCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGt                                                                  <  2:239907/140‑1 (MQ=255)                                                       cTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACt                                                              >  2:245381/1‑140 (MQ=255)                                                          ttAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTAc                                                                             >  1:352535/1‑122 (MQ=255)                                                          ttAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTAc                                                                             <  2:352535/122‑1 (MQ=255)                                                                    aaGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGa                                                 >  2:2885/1‑140 (MQ=255)                                                                      ggTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACt                                               <  1:95601/140‑1 (MQ=255)                                                                         aaaGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGa                                            <  1:411712/140‑1 (MQ=255)                                                                           aGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATg                                          <  2:295576/140‑1 (MQ=255)                                                                                ccGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTAtt                                     >  1:94981/1‑140 (MQ=255)                                                                                       ttttGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTc                              <  1:222257/140‑1 (MQ=255)                                                                                           gCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCTGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCt                          >  1:86532/1‑140 (MQ=255)                                                                                            cccGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTg                         <  2:202954/140‑1 (MQ=255)                                                                                             ccGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTgg                        <  1:35624/140‑1 (MQ=255)                                                                                                       cgTACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTAt              >  1:9579/1‑140 (MQ=255)                                                                                                         tACCGAAGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCa            <  1:245381/140‑1 (MQ=255)                                                                                                              aaGTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGc       >  2:48130/1‑140 (MQ=255)                                                                                                                gTCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCt     <  2:163743/140‑1 (MQ=255)                                                                                                                 tCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGc       <  1:382599/137‑1 (MQ=255)                                                                                                                 tCCGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGc       >  2:382599/1‑137 (MQ=255)                                                                                                                   cGTAAGGTCATTTGGCCGGCTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAg  >  1:171797/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                          CCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCGCTGTTAACGACAAAAGGTAAAGCTACCGTTGCTTTTGCCCGTGAAGCGCGTACCGAAGTCCGTAAGGTCATTTGGCCGACTCGCCAGGAAACATTGCACACCACGCTGATTGTGGCTGCGGTTACCGCAGTAATGTCACTGATCCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAG  >  NC_000913/4177481‑4177734

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