Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,283,902 A→G 8.1% L108P (CTG→CCG)  PP_2885 ← transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,283,9020AG 8.1% 135.7 / 2.8 49L108P (CTG→CCG) PP_2885transporter
Reads supporting (aligned to +/- strand):  ref base A (24/21);  new base G (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

CGGCTCCACATGGTCGGCGATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCACGCCAGCCACGGTCACGGCGGCGCCGGCTT  >  NC_002947/3283760‑3284046
                                                                                                                                              |                                                                                                                                                
cGGCTCCACATGGTCGGCGATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGcaca                                                                                                                                            <  1:186875/149‑1 (MQ=255)
  gCTCCACATGGTCGGCGATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCAcag                                                                                                                                           >  2:119724/1‑148 (MQ=255)
    tCCACATGGTCGGCGATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACa                                                                                                                                                  <  2:46952/139‑1 (MQ=255)
    tCCACATGGTCGGCGATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACa                                                                                                                                                  >  1:46952/1‑139 (MQ=255)
      cacaTGGTCGGCGATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACagcagc                                                                                                                                       <  1:188519/148‑1 (MQ=255)
         aTGGTCGGCGATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAg                                                                                                                                    >  6:170326/1‑148 (MQ=255)
                 cGATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAgg                                                                                                                           <  4:343078/149‑1 (MQ=255)
                  gATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAggg                                                                                                                          >  6:223925/1‑149 (MQ=255)
                           aTTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAACGCCgg                                                                                                                 >  5:326303/1‑149 (MQ=255)
                           aTTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCgg                                                                                                                 >  5:332097/1‑149 (MQ=255)
                           aTTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCgg                                                                                                                 <  5:39422/149‑1 (MQ=255)
                            ttGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGgc                                                                                                                <  6:326303/149‑1 (MQ=255)
                            ttGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGgc                                                                                                                <  5:223925/149‑1 (MQ=255)
                            ttGTAGACTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGgc                                                                                                                <  6:332097/149‑1 (MQ=255)
                              gTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAgg                                                                                                                           <  7:383405/136‑1 (MQ=255)
                              gTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAgg                                                                                                                           >  8:383405/1‑136 (MQ=255)
                                      ccGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCcagc                                                                                                      <  7:80812/149‑1 (MQ=255)
                                             ggcggcATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGa                                                                                               >  4:149920/1‑149 (MQ=255)
                                              gcggcATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGa                                                                                               <  1:119724/148‑1 (MQ=255)
                                                   aTAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGcc                                                                                         >  2:122960/1‑149 (MQ=255)
                                                   aTAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGACGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGTCGCCTTCGACGAAAAGGCAAAGCAGCCAGTAGGCAAGGGTaaa                                                                                                                      >  2:333677/1‑120 (MQ=255)
                                                            gcCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGcc                                                                                  >  5:359923/1‑147 (MQ=255)
                                                              cGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACGGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAgg                                                                              >  1:15516/1‑149 (MQ=255)
                                                               gATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGc                                                                             >  7:72580/1‑149 (MQ=255)
                                                                    aaCCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCt                                                                        >  2:197365/1‑149 (MQ=255)
                                                                         cgcgcgATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCcgc                                                                   <  3:149920/149‑1 (MQ=255)
                                                                         cgcgcgATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCcgc                                                                   <  5:170326/149‑1 (MQ=255)
                                                                          gcgcgATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCcgcg                                                                  >  3:276338/1‑149 (MQ=255)
                                                                          gcgcgATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCcgcg                                                                  >  5:231160/1‑149 (MQ=255)
                                                                            gcgATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCcgcgcg                                                                <  1:122960/149‑1 (MQ=255)
                                                                             cgATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCgg                                                               <  6:231160/149‑1 (MQ=255)
                                                                             cgATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGGCTGCCGCGCgg                                                               <  4:276338/149‑1 (MQ=255)
                                                                              gATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGc                                                              >  1:188270/1‑149 (MQ=255)
                                                                                  atcaGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGccc                                                            >  7:226977/1‑147 (MQ=255)
                                                                                   tcaGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAg                                                         <  8:72580/149‑1 (MQ=255)
                                                                                        ccGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAgggccg                                                    <  2:82884/149‑1 (MQ=255)
                                                                                                cGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCcagca                                            <  1:197365/149‑1 (MQ=255)
                                                                                                   cAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAgggccgg                                                   <  3:6571/139‑1 (MQ=255)
                                                                                                   cAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAgggccgg                                                   >  4:6571/1‑139 (MQ=255)
                                                                                                     aGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACGGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCAcc                                       <  2:15516/149‑1 (MQ=255)
                                                                                                            ccccGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCa                                >  6:390694/1‑149 (MQ=255)
                                                                                                                       aGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCcagca                                            >  2:352053/1‑126 (MQ=255)
                                                                                                                              ggCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCACgcca                           <  1:79827/136‑1 (MQ=255)
                                                                                                                              ggCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCACgcca                           >  2:79827/1‑136 (MQ=255)
                                                                                                                               gcgcCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCACGCCAGCCACGGTCAcggc             >  7:250856/1‑149 (MQ=255)
                                                                                                                                   cTTCGACGAACGGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCAc                               >  2:242774/1‑127 (MQ=255)
                                                                                                                                   cTTCGACGAACGGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCAc                               <  1:242774/127‑1 (MQ=255)
                                                                                                                                          gaACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCACGCCAGCCACGGTCACGGCGGCGCCGGCtt  >  2:191494/1‑149 (MQ=255)
                                                                                                                                          gaACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCACGCCAGCCACGGTCACGGCGGCGCCGGCtt  >  6:111556/1‑149 (MQ=255)
                                                                                                                                            aCAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAgcagca                                         >  3:210251/1‑108 (MQ=255)
                                                                                                                                            aCAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAgcagca                                         <  4:210251/108‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CGGCTCCACATGGTCGGCGATCAGCGCATTGTAGCCTACCGGGATGGCGGCATAAAAGGCGCCGATGCAACCGCGCGCGATCATCAGGCCGATGAACGCCAAGGTCGCCCCCGGCATCCAGCGCAGGGCGCCTTCGACGAACAGGCACAGCAGCCAGTAGGCCAGGGTAAAGCCGGCGCTGCCCAGCAGCAGGATGCGCCGGCGGCCCAGGCGGTCTGCCGCGCGGCCCCAGGGCCGGGCCAGCAGCACCCAGACCACGCCAGCCACGGTCACGGCGGCGCCGGCTT  >  NC_002947/3283760‑3284046

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