Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,142,007 G→C 8.9% M114I (ATG→ATC)  PP_4526 → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,142,0070GC 8.9% 133.1 / 4.6 45M114I (ATG→ATC) PP_4526MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (24/17);  new base C (2/2);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

GCCAGGTACGGCGACGCCTGCATGGGGGGCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTGGCTCGGCAGCCCTGT  >  NC_002947/5141868‑5142155
                                                                                                                                           |                                                                                                                                                     
gccaggTACGGCGACGCCTGCATGGGGGGCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGc                                                                                                                                              <  5:195463/149‑1 (MQ=255)
gccaggTACGGCGACGCCTGCATGGGGGGCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGc                                                                                                                                              <  2:30145/149‑1 (MQ=255)
 ccaggTACGGCGACGCCTGCATGGGGGGCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGGCAACGGCTTTTGTGGCCACCTGCTGCTGCGCTTTGGCACATGCTTAGCCAGCGCCTGGGTGCTGGTGATCATCACCAGca                                                                                                                                             <  4:13939/149‑2 (MQ=255)
  caggTACGGCGACGCCTGCATGGGGGGCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATCATCACCAGcac                                                                                                                                            <  5:137483/149‑3 (MQ=255)
               gCCTGCATGGGGGGCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTa                                                                                                                               <  8:74310/149‑1 (MQ=255)
                ccTGCATGGGTGGATGTGGCTGTGTGTTCTGCTGACCGTTTACTCCTGGGCTGCCAACGGCTTTTGGGGCCTCCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCTCCTGTGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAg                                                                                                                              <  5:91460/149‑1 (MQ=255)
                      tGGGGGGCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTg                                                                                                                        >  8:299670/1‑149 (MQ=255)
                        gggggCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTgc                                                                                                                       <  7:299670/148‑1 (MQ=255)
                          gggCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGca                                                                                                                    >  7:158393/1‑149 (MQ=255)
                             ctgtggctgtgTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACaa                                                                                                                 >  3:128532/1‑149 (MQ=255)
                                tggctgtgTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCg                                                                                                              >  2:57685/1‑149 (MQ=255)
                                  gctgtgTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGc                                                                                                             >  2:130146/1‑148 (MQ=255)
                                                tgACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGcac                                                                                                                   <  6:310815/128‑1 (MQ=255)
                                                tgACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGcac                                                                                                                   >  5:310815/1‑128 (MQ=255)
                                                   cccTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTg                                                                                          >  1:224001/1‑149 (MQ=255)
                                                    ccTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTgg                                                                                         >  4:122970/1‑149 (MQ=255)
                                                        gCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTAtt                                                                                     >  6:243400/1‑149 (MQ=255)
                                                                   tGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCg                                                                          >  6:261812/1‑149 (MQ=255)
                                                                      cAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGcc                                                                       >  8:43640/1‑149 (MQ=255)
                                                                      cAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGcc                                                                       >  5:284356/1‑149 (MQ=255)
                                                                       aaCGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCt                                                                      >  3:36440/1‑149 (MQ=255)
                                                                             ttttGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTAt                                                                >  4:69304/1‑149 (MQ=255)
                                                                                tGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCTAGGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATc                                                               <  1:57685/148‑1 (MQ=255)
                                                                                  gggCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGcc                                                            <  2:224001/148‑1 (MQ=255)
                                                                                   ggCCACCTGCTGCTGCGCTTTTGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGt                                                          <  8:198787/149‑1 (MQ=255)
                                                                                    gCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGt                                                          >  8:200011/1‑148 (MQ=255)
                                                                                        ccTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCg                                                     >  8:222748/1‑149 (MQ=255)
                                                                                            ctgctgcGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCggg                                                   <  8:158393/147‑1 (MQ=255)
                                                                                            ctgctgcGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGGCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGc                                                  >  8:111946/1‑148 (MQ=255)
                                                                                                    cTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCAc                                         <  3:122970/149‑1 (MQ=255)
                                                                                                     tttGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACt                                        >  2:278708/1‑149 (MQ=255)
                                                                                                      ttGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGa                                                        >  6:265966/1‑132 (MQ=255)
                                                                                                      ttGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGa                                                        <  5:265966/132‑1 (MQ=255)
                                                                                                       tGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTg                                       <  7:222748/148‑1 (MQ=255)
                                                                                                               gCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCAt                              >  2:219101/1‑149 (MQ=255)
                                                                                                                 gTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCAtgt                            <  5:243400/149‑1 (MQ=255)
                                                                                                                  tcaccAGCACCCAGGCGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCAtgtg                           >  3:154311/4‑149 (MQ=255)
                                                                                                                          gcCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCt                   >  6:244031/1‑149 (MQ=255)
                                                                                                                            gctgaggCTGGTGATCATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTg                 >  6:126666/7‑149 (MQ=255)
                                                                                                                            cTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTg                 >  8:197769/1‑149 (MQ=255)
                                                                                                                                gTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTGGCTc             <  5:300284/149‑1 (MQ=255)
                                                                                                                                aggCTGGTGATCATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTGGCTc             >  4:38436/3‑149 (MQ=255)
                                                                                                                                   cTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTGGCTCgg           <  1:278708/148‑1 (MQ=255)
                                                                                                                                    tGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCg                                                                    <  7:126733/90‑1 (MQ=255)
                                                                                                                                    tGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCg                                                                    >  8:126733/1‑90 (MQ=255)
                                                                                                                                       tgatgaTCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTGGCTCGGCAGcc      <  4:36440/149‑1 (MQ=255)
                                                                                                                                          tgatCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTGGCTCGGCAGCCCtg   <  6:284356/149‑1 (MQ=255)
                                                                                                                                           gatCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGc                                                                             <  4:179217/74‑1 (MQ=255)
                                                                                                                                           gatCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGc                                                                             >  3:179217/1‑74 (MQ=255)
                                                                                                                                           gatCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCgg                                                    >  6:272674/1‑99 (MQ=255)
                                                                                                                                           gatCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCgg                                                    <  5:272674/99‑1 (MQ=255)
                                                                                                                                           gatCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTGGCTCGGCAGCCCtgt  >  7:269730/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
GCCAGGTACGGCGACGCCTGCATGGGGGGCTGTGGCTGTGTGTGCTGCTGACCCTGGCCTCCTGGGCTGCCAACGGCTTTTGGGGCCACCTGCTGCTGCGCTTTGGCACAGGCGTAGCCAGCGCCTGGGTGCTGGTGATGATCACCAGCCTCAGCCAGCAAGTAGCCACTGCGCACAACCGCCAACGTCT‑GGGAGCGCTGGTATTCGCCGGGCCCGGCCTCGGTATCGCCGTGACCGGGCTACTGGCACTGGTTGCGCATGTGCTGGGGCTTGGCTCGGCAGCCCTGT  >  NC_002947/5141868‑5142155

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: ‑