Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,044,678 C→G 9.3% G344R (GGG→CGG) ‡ PP_3566 ← MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,044,6780CG 9.3% 129.2 / 4.6 43G344R (GGG→CGG) ‡PP_3566MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (16/23);  new base G (2/2);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGGCCGCCGGCGC  >  NC_002947/4044541‑4044825
                                                                                                                                          |                                                                                                                                                   
cAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCGGGCGGCTgaatg                                                                                                                                          <  7:268408/149‑5 (MQ=255)
 aGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCGGGCGGCTgaatgc                                                                                                                                         <  8:121930/149‑6 (MQ=255)
  gcgcCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCggg                                                                                                                                        <  2:56479/149‑1 (MQ=255)
        cacCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGTGTCCAGt                                                                                                                                  <  2:141418/149‑1 (MQ=255)
                aaCCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTacac                                                                                                                          <  6:384952/149‑1 (MQ=255)
                aaCCCCGCCCCCACCACTGCCCCGCCTTGCCAGCCAAACTGGCTGAACCGCGGACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGCAGATGCTGGTCATGTACACCGCATTTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTaca                                                                                                                           <  1:67069/149‑1 (MQ=255)
                  ccccGCCCCCCCCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCCCCCGCGGCTGGCCGGGTCCAGTGCGTACACt                                                                                                                         >  2:15699/1‑148 (MQ=255)
                          ccaccaGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGAGCGTACACTTCACGCTgg                                                                                                                <  7:28100/149‑1 (MQ=255)
                           caccaGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGc                                                                                                               <  3:50167/149‑1 (MQ=255)
                            accaGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGcc                                                                                                              <  3:357310/149‑1 (MQ=255)
                                     aCGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCAt                                                                                                     >  7:223746/1‑149 (MQ=255)
                                     aCGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCAt                                                                                                     >  3:311591/1‑149 (MQ=255)
                                     aCGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCAt                                                                                                     <  4:311591/149‑1 (MQ=255)
                                       gCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGt                                                                                                   <  2:134090/149‑1 (MQ=255)
                                                  caAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCAcc                                                                                        <  1:248148/149‑1 (MQ=255)
                                                             aaCAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAg                                                                             >  7:352917/1‑149 (MQ=255)
                                                             aaCAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAg                                                                             >  6:194055/1‑149 (MQ=255)
                                                                       tGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTaca                                                                                                                           <  3:339383/93‑1 (MQ=255)
                                                                       tGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTaca                                                                                                                           >  4:339383/1‑93 (MQ=255)
                                                                         gCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCg                                                                   <  7:62345/147‑1 (MQ=255)
                                                                         gCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCg                                                                   <  8:223746/147‑1 (MQ=255)
                                                                            aCCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTc                                                                                                 >  4:18474/1‑114 (MQ=255)
                                                                            aCCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTc                                                                                                 <  3:18474/114‑1 (MQ=255)
                                                                                 cgACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCa                                                         <  1:346832/149‑1 (MQ=255)
                                                                                            ccccGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCATCACGCCGGTAAGGCCCAGGCCAATCACg                                              <  1:211517/149‑1 (MQ=255)
                                                                                            ccccGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACg                                              <  5:319783/149‑1 (MQ=255)
                                                                                             cccGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCCGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGc                                             >  6:323068/1‑149 (MQ=255)
                                                                                              ccgccgAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCt                                            <  8:352917/149‑1 (MQ=255)
                                                                                               cgccgAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAgg                                                             <  5:61391/131‑1 (MQ=255)
                                                                                               cgccgAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAgg                                                             >  6:61391/1‑131 (MQ=255)
                                                                                                      gAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCg                                    >  5:156212/1‑149 (MQ=255)
                                                                                                             cTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCTGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGt                             >  8:106524/1‑149 (MQ=255)
                                                                                                               ggTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCAc                                               <  8:93699/129‑1 (MQ=255)
                                                                                                               ggTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCAc                                               >  7:93699/1‑129 (MQ=255)
                                                                                                                    tGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGcccc                      <  7:259245/149‑1 (MQ=255)
                                                                                                                    tGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGccc                       >  1:13703/1‑148 (MQ=255)
                                                                                                                      tACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCa                     <  5:194055/148‑1 (MQ=255)
                                                                                                                       acacCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAg                    >  5:356186/1‑148 (MQ=255)
                                                                                                                           cGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGcaac                >  2:195651/1‑148 (MQ=255)
                                                                                                                             cA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAg             <  3:338961/149‑1 (MQ=255)
                                                                                                                               ggccAGCCGCGGGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAgg            >  8:117131/4‑149 (MQ=255)
                                                                                                                                 ccAGCCGCGGGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGgcc          >  7:243586/2‑149 (MQ=255)
                                                                                                                                    gccgccCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGgccgcc       >  1:291991/1‑149 (MQ=255)
                                                                                                                                      cgcccgcGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGGCCGCCgg     >  5:8368/1‑149 (MQ=255)
                                                                                                                                         ccgcGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGGCCGCCGgcgc  <  6:8368/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
CAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCA‑TTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGGCCGCCGGCGC  >  NC_002947/4044541‑4044825

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