Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,389,226 G→A L49F (CTC→TTC)  mscM ← miniconductance mechanosensitive channel MscM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,389,2260GA96.8% 98.7 / NA 31L49F (CTC→TTC) mscMminiconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/11);  minor base C (0/1);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAA  >  NC_000913/4389110‑4389353
                                                                                                                    |                                                                                                                               
tAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTcc                                                                                                                 >  2:389597/1‑133 (MQ=255)
      tgCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCtgt                                                                                                           >  2:318442/1‑133 (MQ=255)
      tgCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCtgt                                                                                                           >  2:58708/1‑133 (MQ=255)
      tgCGCGCAGAGCAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCtgt                                                                                                           >  1:385110/1‑133 (MQ=255)
        cgcgcAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGc                                                                                                         >  2:599912/1‑133 (MQ=255)
        cgcgcAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGc                                                                                                         >  1:520105/1‑133 (MQ=255)
        cgcgcAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGc                                                                                                         >  1:519990/1‑133 (MQ=255)
        cgcgcAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGc                                                                                                         >  2:566799/1‑133 (MQ=255)
                         agTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTg                                                                                        >  2:127626/1‑133 (MQ=255)
                          gTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGc                                                                                       <  1:58708/133‑1 (MQ=255)
                           tttCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGcc                                                                                      <  2:520105/133‑1 (MQ=255)
                            ttCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCt                                                                                     <  2:488696/133‑1 (MQ=255)
                                gATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTc                                                                                 <  2:387838/133‑1 (MQ=255)
                                         cGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGTCTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCGGCTCCAGTTCCTg                                                                        <  2:385110/133‑1 (MQ=255)
                                                  gCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCG‑TCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGAtttgt                                                              >  2:333816/1‑133 (MQ=255)
                                                      atatTGTTTGATGCGATCAAGGGACCCTCTTCTTTCCTCAACCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGAtttgtttg                                                           <  2:251101/133‑1 (MQ=255)
                                                         ttGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTa                                                        >  2:313826/1‑133 (MQ=255)
                                                              tGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTTCAGTTCCGGAGGGACTTGTTTGCTATCggg                                                   >  1:186755/1‑133 (MQ=255)
                                                                 tGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGc                                                >  1:395706/1‑133 (MQ=255)
                                                                 tGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGc                                                >  2:8214/1‑133 (MQ=255)
                                                                             aaCCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTa                                    <  1:566799/133‑1 (MQ=255)
                                                                             aaCCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTa                                    <  1:318442/133‑1 (MQ=255)
                                                                                     cGTTCCTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGcccccc                            <  1:599912/133‑1 (MQ=255)
                                                                                          cTCAAGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTg                       >  1:75912/1‑133 (MQ=255)
                                                                                             aaGCGCATTTAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAgg                    >  1:96080/1‑133 (MQ=255)
                                                                                                     ttAAGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGc            <  1:389597/133‑1 (MQ=255)
                                                                                                        aGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCAt         <  2:519990/133‑1 (MQ=255)
                                                                                                        aGGCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCAt         <  1:127626/133‑1 (MQ=255)
                                                                                                         ggCAGACTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCACCAACTCAGGCACAAGGCCATc        >  1:298337/1‑133 (MQ=255)
                                                                                                            aGACTGGACCGCCTCTACGACTCCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGa     <  1:284597/133‑1 (MQ=255)
                                                                                                              aCTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGaaa   >  2:82823/1‑133 (MQ=255)
                                                                                                               cTGGAACGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCCCGGCGTACGCCCTCCAACTGAGGCACCAGGCAATCAGaaaa  >  2:284597/1‑133 (MQ=255)
                                                                                                                    |                                                                                                                               
TAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAA  >  NC_000913/4389110‑4389353

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