Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,386,872	C→T	A833A (GCG→GCA)	mscM ←	mechanosensitive channel protein, miniconductance

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,386,872	0	C	T	100.0%	129.1 / NA	40	A833A (GCG→GCA)	mscM	mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/20);  total (20/20)

GTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATC  >  NC_000913/4386741‑4387007                                                                                                                                    |                                                                                                                                        gtgatggtggtgatggCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGa                                                                                                                                  >  2:193129/1‑139 (MQ=255)                   tAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTCATTGGCTCCAGAc                                                                                                                >  2:1116810/1‑139 (MQ=255)                   tAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGAc                                                                                                                <  1:193129/139‑1 (MQ=255)                        cGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                           <  2:337536/139‑1 (MQ=255)                          tACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTAc                                                                                                         >  1:794007/1‑139 (MQ=255)                          tACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTAc                                                                                                         >  2:166044/1‑139 (MQ=255)                           aCCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAATCACCAGAATGGCAATCCGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                        <  2:277907/139‑1 (MQ=255)                           aCCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                        <  2:794007/139‑1 (MQ=255)                               ggCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                           >  1:1070241/1‑132 (MQ=255)                               ggCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                           <  2:1070241/132‑1 (MQ=255)                               ggCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCt                                                                                                    >  1:814170/1‑139 (MQ=255)                                  gTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTa                                                                                                 <  1:166044/139‑1 (MQ=255)                                        tCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTgg                                                                                           <  1:437582/139‑1 (MQ=255)                                             gTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTg                                                                                      >  1:1063816/1‑139 (MQ=255)                                                ctgcAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                           >  2:944384/1‑115 (MQ=255)                                                ctgcAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                           <  1:944384/115‑1 (MQ=255)                                                ctgcAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa                                                                                   >  1:669728/1‑139 (MQ=255)                                                        tCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGa                                                                                          <  1:333007/124‑1 (MQ=255)                                                        tCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGa                                                                                          >  2:333007/1‑124 (MQ=255)                                                                          cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGaa                                                         >  1:677761/1‑139 (MQ=255)                                                                                    gcgcACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGcc                                                                                                      >  1:466154/1‑84 (MQ=255)                                                                                    gcgcACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGcc                                                                                                      <  2:466154/84‑1 (MQ=255)                                                                                           aGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGa                                        <  1:685523/139‑1 (MQ=255)                                                                                            gCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGa                                        <  1:946691/138‑1 (MQ=255)                                                                                            gCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGa                                        >  2:946691/1‑138 (MQ=255)                                                                                                   ggtgATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAc                                <  2:677761/139‑1 (MQ=255)                                                                                                   ggtgATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATCTTCGAGGAAGTCGAAAGCGGAATGGATTTCTGAc                                <  2:1063816/139‑1 (MQ=255)                                                                                                       atgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                           >  2:136909/1‑60 (MQ=255)                                                                                                       atgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                           <  1:136909/60‑1 (MQ=255)                                                                                                       atgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACcaca                            >  2:293837/1‑139 (MQ=255)                                                                                                          atAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGca                         <  1:293837/139‑1 (MQ=255)                                                                                                                 ccAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGAc                  <  2:814170/139‑1 (MQ=255)                                                                                                                  cAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                        <  2:1018132/52‑1 (MQ=255)                                                                                                                  cAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                        >  1:1018132/1‑52 (MQ=255)                                                                                                                   aGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGcc                                                                                                      <  1:416678/53‑1 (MQ=255)                                                                                                                   aGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGcc                                                                                                      >  2:416678/1‑53 (MQ=255)                                                                                                                          cAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAg         >  1:1010050/1‑139 (MQ=255)                                                                                                                               aGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                           <  2:1097578/36‑1 (MQ=255)                                                                                                                               aGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                           >  1:1097578/1‑36 (MQ=255)                                                                                                                               aGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGa    >  1:867712/1‑139 (MQ=255)                                                                                                                                 cactgcaccgAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATc  <  2:1010050/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        GTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATC  >  NC_000913/4386741‑4387007

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