Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,165,671 C→T 18.9% G159R (GGG→AGG)  PP_3666 ← MHS family metabolite MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,165,6710CT18.9% 71.9 / 11.9 37G159R (GGG→AGG) PP_3666MHS family metabolite MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (16/14);  new base T (5/2);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

CAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGG  >  NC_002947/4165531‑4165817
                                                                                                                                            |                                                                                                                                                  
cagGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGgcgc                                                                                                                                           <  2:208066/150‑1 (MQ=255)
 agGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCGCGCCGCCTGGgcggcgtg                                                                                                                                          <  1:398652/150‑3 (MQ=255)
                                 ggCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCa                                                                                                                                                   >  2:515801/1‑109 (MQ=255)
                                 ggCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCa                                                                                                                                                   <  1:515801/109‑1 (MQ=255)
                                 ggCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCGCGGCGTa                                                                                                          >  1:281567/1‑150 (MQ=255)
                                    ttCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGc                                                                                                                         >  1:280477/1‑132 (MQ=255)
                                    ttCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGc                                                                                                                         <  2:280477/132‑1 (MQ=255)
                                    ttCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTc                                                                                                       >  2:340038/1‑150 (MQ=255)
                                                    tCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGTCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAc                                                                                       >  1:498233/1‑150 (MQ=255)
                                                     cATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCGCGCCGCCTGGgcggcgtgg                                                                                                                                         >  1:247726/1‑96 (MQ=255)
                                                     cATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCGCGCCGCCTGGgcggcgtgg                                                                                                                                         <  2:247726/99‑4 (MQ=255)
                                                      aTGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAccc                                                                                     <  1:168241/150‑1 (MQ=255)
                                                       tGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCACCGGCGTACTCACCCGAGGCGGAAAAcccc                                                                                    >  2:363066/1‑150 (MQ=255)
                                                                          ccgcACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagca                                                                 >  2:139848/1‑150 (MQ=255)
                                                                              aCCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacag                                                             >  2:135830/1‑150 (MQ=255)
                                                                               ccGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacagc                                                            >  2:136058/1‑150 (MQ=255)
                                                                                cGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacagca                                                           >  2:30814/1‑150 (MQ=255)
                                                                                    caacaacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGc                                                       <  1:30814/150‑1 (MQ=255)
                                                                                      acaacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCg                                                     >  2:81722/1‑150 (MQ=255)
                                                                                       caacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCgg                                                    <  2:281567/150‑1 (MQ=255)
                                                                                           cccGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCg                                                                                                  >  1:441355/1‑100 (MQ=255)
                                                                                           cccGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCg                                                                                                  <  2:441355/100‑1 (MQ=255)
                                                                                                 gccgTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGc                                          <  2:498233/150‑1 (MQ=255)
                                                                                                 atcgTGCCGGCCAGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGc                                          >  2:777492/3‑150 (MQ=255)
                                                                                                       cTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATc                                    <  2:687076/150‑1 (MQ=255)
                                                                                                          gCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTgc                                 <  1:340038/150‑1 (MQ=255)
                                                                                                          gCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACGCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTgc                                 <  1:571095/150‑1 (MQ=255)
                                                                                                                aCGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAc                                                                                       >  1:504210/1‑90 (MQ=255)
                                                                                                                aCGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAc                                                                                       <  2:504210/90‑1 (MQ=255)
                                                                                                                            tACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAg               >  2:464306/1‑150 (MQ=255)
                                                                                                                             aCAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGTTCATCGGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGc              <  2:707005/150‑1 (MQ=255)
                                                                                                                              gagtacgcgccgcCTGGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGcc             >  1:206007/8‑150 (MQ=255)
                                                                                                                                 tacgcgccgcCTGGGCGGCGTGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGAt          >  1:728207/5‑150 (MQ=255)
                                                                                                                                 tacgcgccgcCTGGGCGGCGTGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGAt          >  2:205599/5‑150 (MQ=255)
                                                                                                                                    gcgccgcCTGGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCa       >  1:545755/2‑150 (MQ=255)
                                                                                                                                      gccgccCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGa     >  1:518882/1‑150 (MQ=255)
                                                                                                                                         gccCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAgg  <  1:81722/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGCCGATGCAGAAGG  >  NC_002947/4165531‑4165817

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