Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	487,023	C→T	T163M (ACG→ATG)	mscK →	potassium dependent, small conductance mechanosensitive channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	487,023	0	C	T	100.0%	96.6 / NA	29	T163M (ACG→ATG)	mscK	potassium dependent, small conductance mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/14);  total (15/14)

CGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCAAAAGCAACGTGATTACGTGACG  >  NC_000913/486778‑487255                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         gaacacataaaaaaaAGAGACA‑TTCAACCACCGCCAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGcc                                                                                                                                                                                                                                    <  2:335771/242‑1 (MQ=255) cgatcgCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGc                                                                                                                                                                                                                                     >  1:88053/1‑251 (MQ=255)                    aCAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCCGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGc                                                                                                                                                                                                                  >  2:269792/1‑251 (MQ=255)                               aCGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAAc                                                                                                                                                                                                        <  2:126469/250‑1 (MQ=255)                               aCGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACg                                                                                                                                                                                                       >  2:114281/1‑251 (MQ=255)                               aCGGCAAAAAGGCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCTGGCGCTGGACGATTTGCAAAACGCACAAAACTATCTGGCGTCTTATAACATCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATaa                                                                                                                                                                                                         <  1:249554/249‑1 (MQ=255)                                                      cGGAAAAAATGCGCCAGGCGACCGCGGCGTTACCAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCCACAAAtt                                                                                                                                                                                >  2:138141/1‑251 (MQ=255)                                                       ggAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTc                                                                                                                                                                               <  1:147197/251‑1 (MQ=255)                                                            aaaTGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAgt                                                                                                                                                                          <  2:122551/251‑1 (MQ=255)                                                                              cggcgTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGAt                                                                                                                                                        <  2:359443/251‑1 (MQ=255)                                                                                          cACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGAc                                                                                                                                             >  2:269386/1‑250 (MQ=255)                                                                                                          tAACGACGACGACACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCccagcc                                                                                                                            >  1:301973/1‑251 (MQ=255)                                                                                                                         gcgcAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGc                                                                                                             <  2:81042/251‑1 (MQ=255)                                                                                                                            cAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTCATGCAgg                                                                                                          >  2:301605/1‑251 (MQ=255)                                                                                                                                  tCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCcag                                                                                                     <  2:52391/250‑1 (MQ=255)                                                                                                                                     gAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAgg                                                                                                 <  2:11198/251‑1 (MQ=255)                                                                                                                                              gTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTg                                                                                         >  1:298499/1‑250 (MQ=255)                                                                                                                                                           cTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTg                                                                            >  1:5376/1‑250 (MQ=255)                                                                                                                                                            tGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAgt                                                                                                                     >  1:288206/1‑208 (MQ=255)                                                                                                                                                            tGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAgt                                                                                                                     <  2:288206/208‑1 (MQ=255)                                                                                                                                                                          gCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTaaaa                                                             <  1:240945/250‑1 (MQ=255)                                                                                                                                                                              aGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTg                                                        <  1:114281/251‑1 (MQ=255)                                                                                                                                                                                       aCGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAac                                               >  2:166575/1‑251 (MQ=255)                                                                                                                                                                                            ttGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGACTGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGt                                          >  1:129363/1‑251 (MQ=255)                                                                                                                                                                                             tGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGt                                          <  2:88053/250‑1 (MQ=255)                                                                                                                                                                                                         aaaaCGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATAcc                             >  2:302171/1‑251 (MQ=255)                                                                                                                                                                                                                   ggCGTCTTATAACAGCCAGCTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGTAACACCGTCTTGCAGGATACCTTGCAAAAg                    <  1:287674/250‑1 (MQ=255)                                                                                                                                                                                                                                agccagcTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGTGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCTTCTTGCAGGATACCTTGCAAAAGCAACGTGATTACGt      >  2:93157/1‑251 (MQ=255)                                                                                                                                                                                                                                    agcTGGTTTCGTTACAGATGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCAAAAGCAACGTGATTACGTGACg  >  1:101437/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         CGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCAAAAGCAACGTGATTACGTGACG  >  NC_000913/486778‑487255

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