Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,388,794 C→T A193T (GCG→ACG)  mscM ← miniconductance mechanosensitive channel MscM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,388,7940CT92.9% 30.4 / ‑5.9 14A193T (GCG→ACG) mscMminiconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/8);  minor base G (1/0);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGC  >  NC_000913/4388562‑4389001
                                                                                                                                                                                                                                        |                                                                                                                                                                                                               
gccgATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTCATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTTGCAGCCCGCTGCGTCAGCTCCCGTTCATCAAc                                                                                                                                                                                               >  1:234358/1‑251 (MQ=255)
gccgATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATAGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGAAGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGTCAGc                                                                                                                                                                                                             >  1:216513/1‑237 (MQ=255)
                                                aaaTCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGTCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGc                                                                                                                                               <  1:170744/251‑1 (MQ=255)
                                                                  gCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCGGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCTCCAGCTGTGAGCGTAAGCGCGCTAATGCCTGGCGGTTATTGGCAGACAGCTGCGTCAGCTCCAGTGCATCAACGAGCGCCTGAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGact                                                                                                                                <  1:296976/248‑3 (MQ=255)
                                                                             ggTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGTCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCg                                                                                                                  <  1:19076/251‑1 (MQ=255)
                                                                                                             gNNNTTGTCTNTNNAATGGNTTTNGCACGTCCGGCAAATCCGCATCCAGTTGCTGGCGCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGTCAGCTNCAGNTCATNAANGAGCGCCTTAATACGNGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGNNTNNNGNCGGTAAGc                                                                                                                   <  2:234358/218‑1 (MQ=255)
                                                                                                                           aaTTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGTCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTtgtn                                                                                                                                                <  1:39438/175‑2 (MQ=255)
                                                                                                                               gATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGTCAGCTCCAGTTCATCAACGAGCGCCGTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCg                                                                                                                                    <  1:134726/183‑1 (MQ=255)
                                                                                                                                           cnTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGGCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCCACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCCTTTAACTGACGGCGGGCGTCGGTTTGCTTTTGCGGccgg                                                    >  1:253951/1‑248 (MQ=255)
                                                                                                                                                           cAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGCCAGCTGCGTCAGCTCCCGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTAt                                                                                                                              >  1:72816/1‑161 (MQ=255)
                                                                                                                                                           cAGTTGCNNGCGCTCTTTNTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGNTGCGTCAGCNNCANTTCATNANNNANNNCCNTAAGACGTGCAGAGTCAGANTGNNACGCGANATTTTGTGCCTGATTGAGCGGAGTAt                                                                                                                              <  2:72816/161‑1 (MQ=255)
                                                                                                                                                             gTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAACGCCCTAATTCCTGGCGCTTATTGGAACACCGCTGAGTCAGCACCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGCCAGACTGAAACCCGGAATTTTGTGCCTGCTTGAGCGGAGTATTGCCGGTAAGCGGTCCCAGGGGGCGCTCGATCTCATTTAACTGACGGCGGGCGGCGGGTTGCTGGTGccgatcggaagagcaaccgcctga                                   >  1:212446/1‑227 (MQ=255)
                                                                                                                                                                   ggCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGTCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCTGCAGTTGATGCAGCGAATCGGCAAtctcn                            <  1:81137/251‑2 (MQ=255)
                                                                                                                                                                                             cGTAAGCGCGCTAATTCCTGGCGGTTAATGGCCGACAGCTGCGTCAGCTCCAGCTCATTAACCAGAGCCTTCAAACGTGCCGAGTCAGACTGCAACGCGCACATTTATGCCTGATTGAGCGGAGGATTGCCGGTAAGCATTCCCAGGCGGCGCTCGCTCTCATTTAACTGCCGGCGGGCGTAGGGTTGCTTGTGCGGCAGTTGATTCAACGAATTGGGAATGTCTCGGGCGGGCTCCTgtaaatcgggagg  >  1:236013/1‑239 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                               
GCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGC  >  NC_000913/4388562‑4389001

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