Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,387,120 C→T E751K (GAA→AAA)  mscM ← miniconductance mechanosensitive channel MscM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,387,1200CT100.0% 58.4 / NA 20E751K (GAA→AAA) mscMminiconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/7);  total (13/7)

CTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGAGGTTTCAAGCCTCGCTAACAGCGCCTGTGCCGTTGCCAGATAACCCACCGCCGACGCCAGAATGGCAACCAACGGCGCGCCAATCATCATGTTCCACAGCAT  >  NC_000913/4386890‑4387363
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cTCCAGACTTTCTACGCCCTGTACCGTGGAGGGGACATCCCACAGCGAACTATTTTCCAGGAAGCCGAACGCGGAATGGATTTCTGACCACAGCACAATGACAGCCAGCCGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGtt                                                                                                                                                                                                                                 >  1:161111/1‑251 (MQ=255)
             aCGCCCTGTTCCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGTCTTCAATTGCTCCTTCCGGCCTACTGTGATGATTCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACAtctc                                                                                                                                                                                                                    <  1:331533/251‑1 (MQ=255)
                                 gACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGccc                                                                                                                                                                                                >  1:19990/1‑251 (MQ=255)
                                  aCATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACGAGCCGCAAGGATTGCGCACTGATGGCATCGAGAGCGACTGCGCTTTCATCGACTTCAATTGCTCCGTCCGGGCTACTTTGATGATGCGCTTTCGCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCg                                                                                                                                                                                               <  1:143918/251‑1 (MQ=255)
                                  aCATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAAGGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCTCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCg                                                                                                                                                                                               <  1:222413/251‑1 (MQ=255)
                                               aanTATNTNNGNNGAANCCGAAAGCGGAATGGATTTCTNNCCACAGCACNATGACAGAAAGCAGGGCGATCAACAGCAAAATTGAGCGCACCAGCCGCAAGGANTGCGCACTGNNGGNATCGANANNNANNNCGNTTTCATCGACTTCAATTGCNCCNNCCGGGCNACTGTGATGATGCGCGTTCTCTTCGc                                                                                                                                                                                                                                             <  2:69554/192‑1 (MQ=255)
                                                        ngAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGgcgc                                                                                                                                                                         >  1:84703/2‑251 (MQ=255)
                                                                                  tCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACAACAAAATTGAGCAAACCAGCAGCAAGGATTGAGCAATGATGGCATAGAGATCGACTTCGCTTTAATCGACTTCAATTGCTCATTCAGGGCTACTGTGAAGATGCGCTTTGTCTTCGCAACGCGCACGTTGCGATAACATCTCTGCCCGGCGATGCTTCGcacgctca                                                                                                                                                                                           >  1:208073/1‑200 (MQ=255)
                                                                                     gACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACGTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGt                                                                                                                                             <  1:45092/250‑1 (MQ=255)
                                                                                        cacaGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCCTCGAGATCGACTTCGCTTTCATCGCCTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTCACATCTCTGCCCGGCGATGCTTCGCCCGATCAACCGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGAt                                                                                                                                          >  1:143964/1‑250 (MQ=255)
                                                                                                    ncAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAACCACCAGt                                                                                                                               >  1:31514/2‑249 (MQ=255)
                                                                                                                gNGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGAt                                                                                                                                          >  1:254267/1‑226 (MQ=255)
                                                                                                                                       naCCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCg                                                                                                                                                                                                   >  1:80346/2‑146 (MQ=255)
                                                                                                                                                               gATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGCTGCTGCGCTTTCTCTTCGCCCCGCGCACGTTGCGCTAATATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGCTACACCACCAGTACCAGGAACCAGATGGCAACCGCGGTCTCAAGCCTCGCTAACAACGCCTGTGCCGTTGcc                                                                   >  1:173507/1‑250 (MQ=255)
                                                                                                                                                                        gagaTCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGAGGTTTCAAGCCTCGCTAACAGCGCCTGTGCCGTTGCCAGATAAccc                                                          >  1:103022/1‑250 (MQ=255)
                                                                                                                                                                                 ttCGCTTTCATCGACTTCAATTGCTCNTTCCGGNNTANTGTGATGATGCGCTTTCTCNTCNNNANNNGNCCGTTNCGNNAACATCTCTNCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCNAGCGGCGGGNNACGTGGTAAACCACCAGTAGCNGGANNCNNANGGCNAc                                                                                                           >  2:70349/1‑192 (MQ=255)
                                                                                                                                                                                                               cGGGCTACTGTGATGATGCGCTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATgg                                                                                                               <  1:123113/158‑1 (MQ=255)
                                                                                                                                                                                                               cGGGCNACTGTNANNATGCGCTTTCTCTTCGCCANGCGCACGTTGCGCTANCATCNCTGCCNgg                                                                                                                                                                                                             >  2:123113/1‑64 (MQ=255)
                                                                                                                                                                                                                 ggCTACTGCGATGATGCGCTTTCTCTTCGCCACTCGCACGTTGCGCTAACATCTCTGCCCGGCGATGc                                                                                                                                                                                                       <  1:284577/68‑1 (MQ=255)
                                                                                                                                                                                                                                gcgcTTTCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGCGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGGGGTTTCTAGCCTCGCTAACAGCGCCTGTGCCGTTGCCCGATAACCCACCGCCGACGCCCGAATGGCAACCAACGGCGCCCCCATCATCATGGTCCACCGCat  >  1:6829/1‑250 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCATCGAGATCGACTTCGCTTTCATCGACTTCAATTGCTCCTTCCGGGCTACTGTGATGATGCGCTTCCTCTTCGCCACGCGCACGTTGCGCTAACATCTCTGCCCGGCGATGCTTCGCCCGATCAAACGCCAGCCTGCGGCGCTGGATCAGCATCCAGCGGCGGATAACGTGATAAACCACCAGTAGCAGGAACCAGATGGCAACCGAGGTTTCAAGCCTCGCTAACAGCGCCTGTGCCGTTGCCAGATAACCCACCGCCGACGCCAGAATGGCAACCAACGGCGCGCCAATCATCATGTTCCACAGCAT  >  NC_000913/4386890‑4387363

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