Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,278,282 T→C 9.3% intergenic (+22/‑208) PP_4652 → / → purP membrane protein/adenine:H+ symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,278,2820TC 9.3% 112.6 / 2.4 43intergenic (+22/‑208)PP_4652/purPmembrane protein/adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base T (24/15);  new base C (2/2);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACTGGTTTGGTACAGTGCCGGCGGTGGTGCTGGGCGGGGTGGGCACTTTGGTGGTGACCGGGGTGTGGATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTC  >  NC_002947/5278140‑5278430
                                                                                                                                              |                                                                                                                                                    
ggaCTGGTTTGGTAGTGTGCGGGGGGTGGTGGTGGGCGGGTTGGGCACTTTGGTGGTGACCGGGGTGTGGATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCgg                                                                                                                                                <  4:261480/147‑1 (MQ=255)
       tttGGTACAGTGCCGGCGGTGGTGCTGGGCGGGGTGGGCACTTTGGTGGTGACCGGGGTGTGGATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCAGCTTCGCGGGCAAGcc                                                                                                                                         >  1:156588/1‑149 (MQ=255)
                                     ggggTGGGCACTTTGGTGGTGACCGGGGTGTGGATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAgg                                                                                                           >  5:271966/1‑149 (MQ=255)
                                                tttGGTGGTGACCGGGGTGTGGATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAAc                                                                                                >  4:216763/1‑149 (MQ=255)
                                                  tggtggTGACCGGGGTGTGGATATAAGTGTTCGCTACCCTGGCGAAGCGGGACCGGTTGCATAACGTGTGATGGGTCACCTGTGCGCGCCGCTTCGGGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGaaaatgg                                                                                               <  4:158502/148‑4 (MQ=255)
                                                  tggtggTGACCGGGGTGTGGATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAAc                                                                                                <  3:216763/147‑1 (MQ=255)
                                                                gtgtGGATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGc                                                                                 >  7:6649/1‑148 (MQ=255)
                                                                  gtgGATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCa                                                                               >  2:335934/1‑148 (MQ=255)
                                                                    ggATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTc                                                                                          >  4:256756/1‑135 (MQ=255)
                                                                    ggATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTc                                                                                          <  3:256756/135‑1 (MQ=255)
                                                                     taTAAAAGTGTTCCCTACAGTGGCAACGCGGGACCGGTTGCATAAGGGGTTATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGTCCAGAAAAACCCATTCCCCAACGGCCATAg                                                                            <  2:156588/147‑1 (MQ=255)
                                                                     gATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAgg                                                                           >  5:122282/1‑149 (MQ=255)
                                                                      aTAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGc                                                                          >  6:80612/1‑149 (MQ=255)
                                                                      aTAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGc                                                                          >  3:182305/1‑149 (MQ=255)
                                                                               ttCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCcat                                                                 <  2:128366/149‑1 (MQ=255)
                                                                                  ccTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCcatgca                                                              >  8:200325/1‑149 (MQ=255)
                                                                                  ccTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCCTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCcatgca                                                              >  8:195772/1‑149 (MQ=255)
                                                                                  ccTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCCTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCcatgca                                                              >  8:190410/1‑149 (MQ=255)
                                                                                       ccTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCaa                                                                                                                                            >  1:192215/1‑66 (MQ=255)
                                                                                       ccTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCaa                                                                                                                                            <  2:192215/66‑1 (MQ=255)
                                                                                       ccTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTg                                                         <  6:271966/149‑1 (MQ=255)
                                                                                        cTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTgg                                                        >  8:181785/1‑149 (MQ=255)
                                                                                         tGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAACCCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTgg                                                        >  4:311860/1‑148 (MQ=255)
                                                                                            ccAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTat                                                     >  4:221888/1‑148 (MQ=255)
                                                                                              aaGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTatgat                                                  >  2:2112/1‑149 (MQ=255)
                                                                                                      ccGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCt                                          <  1:335934/149‑1 (MQ=255)
                                                                                                         gTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCtttt                                       >  7:159569/1‑149 (MQ=255)
                                                                                                               tAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTtcctc                                 >  7:25303/1‑149 (MQ=255)
                                                                                                               tAACGGGTGATGGGTCACCTGTGCCCGCCGCCTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTtcctc                                 <  7:195772/149‑1 (MQ=255)
                                                                                                               tAACGGGTGATGGGTCACCTGTGCCCGCCGCCTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTtcctc                                 <  7:190410/149‑1 (MQ=255)
                                                                                                                   gggTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCca                             >  7:138023/1‑149 (MQ=255)
                                                                                                                   gggTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCca                             <  1:2112/149‑1 (MQ=255)
                                                                                                                   gggTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCca                             >  1:325724/1‑149 (MQ=255)
                                                                                                                          gggTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACaaaa                      >  2:225529/1‑149 (MQ=255)
                                                                                                                               aCCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACg                 >  5:249901/1‑149 (MQ=255)
                                                                                                                               aCCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTACCTCCCCACACAAAACCAc                  <  7:181785/148‑1 (MQ=255)
                                                                                                                                ccTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACgg                <  1:196270/149‑1 (MQ=255)
                                                                                                                                    tgCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCa              <  8:159569/147‑1 (MQ=255)
                                                                                                                                      cccGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCg          >  1:337845/1‑149 (MQ=255)
                                                                                                                                      cccGCCGCTTCGCGGGCAAGACCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCg          >  6:53751/1‑149 (MQ=255)
                                                                                                                                       ccgccgCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTACCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCgg         >  8:204860/1‑149 (MQ=255)
                                                                                                                                        cgccgcTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGt        <  2:325724/149‑1 (MQ=255)
                                                                                                                                             cTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGt   <  6:122282/149‑1 (MQ=255)
                                                                                                                                              ttCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTc  >  4:153431/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCACTGGTTTGGTACAGTGCCGGCGGTGGTGCTGGGCGGGGTGGGCACTTTGGTGGTGACCGGGGTGTGGATAAAACTGTTCCCTACCCTGGCCAAGCGGGACCGGTTGCATAACGGGTGATGGGTCACCTGTGCCCGCCGCTTCGCGGGCAAGCCCGCGCCTACGGGGGAACTTCCCCGCACAGGGCCAGAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTC  >  NC_002947/5278140‑5278430

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