Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,901,881 G→A 8.0% A449T (GCC→ACC)  PP_4309 → NCS1 family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,901,8810GA 8.0% 143.9 / 2.1 50A449T (GCC→ACC) PP_4309NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base G (26/20);  new base A (2/2);  total (28/22)
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.06e-01

GGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATCGCGCCACGTAACCGCCATTACCACGACG  >  NC_002947/4901733‑4902028
                                                                                                                                                     |                                                                                                                                                   
gggCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGgccg                                                                                                                                                     >  3:48352/1‑149 (MQ=255)
      ccTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCGCTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCg                                                                                                                                               >  8:147879/1‑149 (MQ=255)
        tGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTg                                                                                                                                             >  6:315190/1‑149 (MQ=255)
        tGCTAGGCCCGCTGTACGGCGTGATCATGACCGACTACTGCTTGCTGCACAAAGGCTGCATCAACGTGCCGCAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTg                                                                                                                                             >  8:78343/1‑149 (MQ=255)
         gCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTgc                                                                                                                                            >  3:323975/1‑149 (MQ=255)
          cTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCGCTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTgc                                                                                                                                            <  7:147879/148‑1 (MQ=255)
            aGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACACGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTgccgg                                                                                                                                         >  7:122625/1‑149 (MQ=255)
                        cGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATGAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTgg                                                                                                                             >  1:246797/1‑149 (MQ=255)
                         ggCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGc                                                                                                                            <  8:122625/149‑1 (MQ=255)
                                        gACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCt                                                                                                             >  1:240773/1‑149 (MQ=255)
                                        gACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCt                                                                                                             >  2:224931/1‑149 (MQ=255)
                                                       cgcAAAGGCTGCATCTACGTGCCGGAGCTGTATACCGAGCACCCGGCCTGCGCCTATCACTACAGCAAGGGGATCATGTTGCGCGCCGTGGCCGCCTTCGTGCCGGCAGCACTGCTGGGCATCGTTCTGGGCCTGGTGCCCAACTTcc                                                                                               <  8:110907/148‑1 (MQ=255)
                                                       cgcAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTc                                                                                                <  1:224931/147‑1 (MQ=255)
                                                              gCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATc                                                                                       <  4:48352/149‑1 (MQ=255)
                                                                tGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATcgc                                                                                     >  5:287934/1‑149 (MQ=255)
                                                                  cATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATcgc                                                                                     <  2:240773/147‑1 (MQ=255)
                                                                     cAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGt                                                                                                          <  1:110101/123‑1 (MQ=255)
                                                                     cAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGt                                                                                                          >  2:110101/1‑123 (MQ=255)
                                                                      aaCGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGtt                                                                               >  2:192103/1‑149 (MQ=255)
                                                                       aCGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTtc                                                                              <  7:78343/149‑1 (MQ=255)
                                                                         gTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTtctc                                                                            <  8:74520/149‑1 (MQ=255)
                                                                            ccGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAg                                                                                             >  1:186215/1‑129 (MQ=255)
                                                                            ccGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAg                                                                                             <  2:186215/129‑1 (MQ=255)
                                                                             cGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTgg                                                                        <  2:246797/149‑1 (MQ=255)
                                                                              ggAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGc                                                                       <  1:192103/149‑1 (MQ=255)
                                                                                   tGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCACCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGAtt                                                                  >  2:14530/1‑149 (MQ=255)
                                                                                   tGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCACCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGAtt                                                                  >  2:20069/1‑149 (MQ=255)
                                                                                     tataCCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGAtt                                                                  <  3:166485/147‑1 (MQ=255)
                                                                                     tataCCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGAtt                                                                  <  6:287934/147‑1 (MQ=255)
                                                                                          cGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCg                                                           >  7:149837/1‑149 (MQ=255)
                                                                                               aCCCGGCCGGCGCCTATCACTACAGCAAGGGACATCAACCTCCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAgg                                                                                            >  1:212665/1‑112 (MQ=255)
                                                                                               aCCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAgg                                                                                            <  2:212665/111‑1 (MQ=255)
                                                                                                cccGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATcg                                                     >  7:175681/1‑149 (MQ=255)
                                                                                                cacggccggcGCCTATCACTACAGCAAGGGAATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGCCCATCGTGCTGGCCCTGCTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGACGGCATcg                                                     >  6:168333/3‑149 (MQ=255)
                                                                                                      cggcGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATcgccgccg                                               >  5:57269/1‑149 (MQ=255)
                                                                                                             tATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTtc                                                                              <  1:55375/111‑1 (MQ=255)
                                                                                                             tATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTtc                                                                              >  2:55375/1‑111 (MQ=255)
                                                                                                               tCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCGGTCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCtata                                        <  5:315190/147‑1 (MQ=255)
                                                                                                               tCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCACCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATAc                                       <  1:14530/148‑1 (MQ=255)
                                                                                                               tCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCACCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATAc                                       <  1:20069/148‑1 (MQ=255)
                                                                                                                     cAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGGCGGCCGCACTGCTGGCCATCGTGCTTGCCCTGGTGCCAAACTTCCAGGGCATGGCGCCGTTCTGCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGGTTa                                <  5:168333/149‑1 (MQ=255)
                                                                                                                      aGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTAt                               >  8:69200/1‑149 (MQ=255)
                                                                                                                        cAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCg                                                           >  4:74370/1‑119 (MQ=255)
                                                                                                                        cAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCg                                                           <  3:74370/119‑1 (MQ=255)
                                                                                                                        cAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCtat                                         >  6:258886/1‑137 (MQ=255)
                                                                                                                        cAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCtat                                         <  5:258886/137‑1 (MQ=255)
                                                                                                                         aaGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATcgc                            >  2:280532/1‑149 (MQ=255)
                                                                                                                             gCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATCGCGCCa                        >  6:265807/1‑149 (MQ=255)
                                                                                                                                    aCCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATCGCGCCACGTAAc                  >  6:35435/1‑149 (MQ=255)
                                                                                                                                          gcgcCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCTTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATCGCGCCACGTAACCGCCAt            >  6:288102/1‑149 (MQ=255)
                                                                                                                                            gcCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATCGCGCCACGTAACCGCCATTa          >  8:278706/1‑149 (MQ=255)
                                                                                                                                                    cgccTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATCGCGCCACGTAACCGCCATTACCacgacg  >  6:89222/1‑149 (MQ=255)
                                                                                                                                                    cgccTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATCGCGCCACGTAACCGCCATTACCacgacg  >  5:224894/1‑149 (MQ=255)
                                                                                                                                                    cgccTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATCGCGCCACGTAACCGCCATTACCacgacg  >  3:113074/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGCCGGAGCTGTATACCGAGCACCCGGCCGGCGCCTATCACTACAGCAAGGGCATCAACCTGCGCGCCGTGGCCGCCTTCGTGCCGGCCGCACTGCTGGCCATCGTGCTGGCCCTGGTGCCCAACTTCCAGGGCATCGCGCCGTTCTCCTGGCTGATTGGTGCCGGCATCGCCGCCGCCCTATACCTGCTTATCGCGCCACGTAACCGCCATTACCACGACG  >  NC_002947/4901733‑4902028

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: