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%	165.2 / NA	51	A833A (GCG→GCA)	mscM	mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/22);  total (29/22)

ATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAG  >  NC_000913/4386744‑4387000                                                                                                                                 |                                                                                                                                 atggtggtgatggCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAggg                                                                                                                        >  1:1253456/1‑139 (MQ=255) atggtggtgatggCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAggg                                                                                                                        >  1:218242/1‑139 (MQ=255)       gtgatggCGTAACCCGTACCCGGCGTTAAATCCACGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCAATGCACCGAGGGTAATTg                                                                                                                  >  2:594056/1‑139 (MQ=255)              cGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAg                                                                                                           >  1:1517040/1‑139 (MQ=255)                     cGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                    <  2:1517040/139‑1 (MQ=255)                             gCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTg                                                                                            >  1:552311/1‑139 (MQ=255)                             gCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTCTCTACGCCCTg                                                                                            >  1:1123334/1‑139 (MQ=255)                              cGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGt                                                                                           >  2:374319/1‑139 (MQ=255)                                        aGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAgg                                                                                 <  2:552311/139‑1 (MQ=255)                                          gTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTg                                                                               >  2:261989/1‑139 (MQ=255)                                          gTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTg                                                                               <  1:374319/139‑1 (MQ=255)                                                    aTCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca                                                                     <  2:1123334/139‑1 (MQ=255)                                                          aGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGa                                                                              <  2:976367/124‑1 (MQ=255)                                                          aGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGa                                                                              >  1:976367/1‑124 (MQ=255)                                                                 gcagcGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGCATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTc                                                        >  2:1154724/1‑139 (MQ=255)                                                                    gcgcGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAg                                                     <  1:261989/139‑1 (MQ=255)                                                                     cgcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAgg                                                    >  1:1545958/1‑139 (MQ=255)                                                                      gcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGa                                                   >  2:950359/1‑139 (MQ=255)                                                                       cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCg                                                       <  2:125972/134‑1 (MQ=255)                                                                       cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCg                                                       >  1:125972/1‑134 (MQ=255)                                                                       cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGaa                                                  >  2:689061/1‑139 (MQ=255)                                                                       cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGaa                                                  <  1:689061/139‑1 (MQ=255)                                                                            aaGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGa                                             <  2:1545958/139‑1 (MQ=255)                                                                                    cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGaa                                     >  1:1617521/1‑139 (MQ=255)                                                                                        aGCTGCGTGGTGATGATAAACACCAGAATGGCA‑‑‑AGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGAttt                              <  2:894144/139‑1 (MQ=255)                                                                                             cGTGGTGATGATAAACACCAGCATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCt                            <  1:1154724/139‑1 (MQ=255)                                                                                                   gatgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACcac                      <  2:130487/139‑1 (MQ=255)                                                                                                     tgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                    >  1:653890/1‑59 (MQ=255)                                                                                                     tgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                    <  2:653890/59‑1 (MQ=255)                                                                                                     tgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                 >  1:1309318/1‑62 (MQ=255)                                                                                                     tgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                 <  2:1309318/62‑1 (MQ=255)                                                                                                       atAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAgg                                                                                 >  2:436376/1‑76 (MQ=255)                                                                                                       atAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAgg                                                                                 <  1:436376/76‑1 (MQ=255)                                                                                                         aaaCACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTCTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGcaca                >  1:474089/1‑139 (MQ=255)                                                                                                            cacCAGAATGGCAATCAGCACTGCACCGAGGGTAAt                                                                                                                    <  2:661493/36‑1 (MQ=255)                                                                                                            cacCAGAATGGCAATCAGCACTGCACCGAGGGTAAt                                                                                                                    >  1:661493/1‑36 (MQ=255)                                                                                                            cacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                    <  2:1056258/52‑1 (MQ=255)                                                                                                            cacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                    >  1:1056258/1‑52 (MQ=255)                                                                                                            cacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAgg                                                                                 >  2:176412/1‑71 (MQ=255)                                                                                                            cacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAgg                                                                                 <  1:176412/71‑1 (MQ=255)                                                                                                            cacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATg             >  1:1226369/1‑139 (MQ=255)                                                                                                                aGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAg         >  2:1723591/1‑139 (MQ=255)                                                                                                                 gAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCg                                                       <  2:641892/92‑1 (MQ=255)                                                                                                                 gAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCg                                                       >  1:641892/1‑92 (MQ=255)                                                                                                                  aaTGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCt                                                                                             <  2:1440487/53‑1 (MQ=255)                                                                                                                  aaTGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCt                                                                                             >  1:1440487/1‑53 (MQ=255)                                                                                                                      gCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGa                                 >  2:200324/1‑109 (MQ=255)                                                                                                                      gCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGa                                 <  1:200324/109‑1 (MQ=255)                                                                                                                       cAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAg                    >  2:109744/1‑121 (MQ=255)                                                                                                                       cAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAg                    <  1:109744/121‑1 (MQ=255)                                                                                                                       cAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAg  >  2:328848/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                 ATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAG  >  NC_000913/4386744‑4387000

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