Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,131,369 C→G 8.5% P83R (CCG→CGG)  PP_2748 → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,131,3690CG 8.5% 135.4 / 4.6 47P83R (CCG→CGG) PP_2748ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base G (2/2);  total (23/24)
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.88e-02

TTCAGCGTCAGGCGCGGCGAGATCTGCGCCCTGATCGGGCCGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAACAGGCCCTGGGCCTGCCACGTGCC  >  NC_002947/3131224‑3131517
                                                                                                                                                 |                                                                                                                                                    
ttCAGCGTCAGGCGCGGCGAGATCTGCGCCCTGATCGGGCCGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCt                                                                                                                                                   >  7:183176/1‑149 (MQ=255)
           gcgcGGCGAGATCTGCGCCCTGATCGGGCCGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGgcggcgc                                                                                                                                        >  2:7962/1‑149 (MQ=255)
                gcgAGATCTGCGCCCTGATCGGGCCGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTg                                                                                                                                   <  1:7962/149‑1 (MQ=255)
                       cTGCGCCCTGATCGGGCCGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATc                                                                                                                             >  2:63282/1‑148 (MQ=255)
                       cTGCGCCCTGATCGGGCCGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATc                                                                                                                             <  1:63282/148‑1 (MQ=255)
                            cccTGATCGGGCCGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCg                                                                                                                        >  4:74038/1‑148 (MQ=255)
                                     ggCCGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTc                                                                                                                 >  5:262177/1‑146 (MQ=255)
                                       ccGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGaa                                                                                                            >  1:201953/1‑149 (MQ=255)
                                         gAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAAc                                                                                                           >  4:18390/1‑148 (MQ=255)
                                          aaTGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCAGCCGCGCCGCCTGCAGCGGGTGCGGCCggc                                                                                                                                 <  1:188531/125‑2 (MQ=25)
                                          aaTGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCAGCCGCGCCGCCTGCAGCGGGTGCGGCCggc                                                                                                                                 >  2:188531/1‑124 (MQ=25)
                                           aTGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAAt                                                                                                        <  8:151493/149‑1 (MQ=255)
                                                      aaGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCaa                                                                                             <  2:201953/149‑1 (MQ=255)
                                                           cTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCaag                                                                                            >  7:105099/1‑145 (MQ=255)
                                                           cTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCaag                                                                                            >  4:320113/1‑145 (MQ=255)
                                                           cTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAagaaga                                                                                        <  3:320113/149‑1 (MQ=255)
                                                              gctgctGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGa                                                                                     <  3:18390/149‑1 (MQ=255)
                                                                tgctgAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAgt                                                                                   <  8:305352/149‑1 (MQ=255)
                                                                  ctgAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAgtgt                                                                                 >  3:309496/1‑149 (MQ=255)
                                                                        aTCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGa                                                                           >  1:326840/1‑149 (MQ=255)
                                                                           cTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACaa                                                                        >  2:116905/1‑149 (MQ=255)
                                                                                         ccGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCg                                                             <  8:105099/146‑1 (MQ=255)
                                                                                         ccGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGt                                                           <  4:309496/148‑1 (MQ=255)
                                                                                         ccGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTc                                                          >  7:331809/1‑149 (MQ=255)
                                                                                         ccGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGCCTGCAGCGGGTGcggccggc                                                                                                                                          <  3:197024/69‑1 (MQ=21)
                                                                                         ccGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGCCTGCAGCGGGTGcggccggc                                                                                                                                          >  4:197024/1‑69 (MQ=21)
                                                                                          cGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCt                                                         >  3:243328/1‑149 (MQ=255)
                                                                                           gggCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGa                                                                                     <  7:53523/120‑1 (MQ=255)
                                                                                           gggCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGa                                                                                     >  8:53523/1‑120 (MQ=255)
                                                                                             gCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTc                                                      <  3:74038/149‑1 (MQ=255)
                                                                                                        ccAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAgc                                           >  7:204839/1‑149 (MQ=255)
                                                                                                                 cttcGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGTGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCa                                             <  1:274120/138‑1 (MQ=255)
                                                                                                                 cttcGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCa                                             >  2:274120/1‑138 (MQ=255)
                                                                                                                  ttcGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACtt                                 >  7:262970/1‑149 (MQ=255)
                                                                                                                     gCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCt                                                                              >  2:321905/1‑101 (MQ=255)
                                                                                                                     gCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCt                                                                              <  1:321905/101‑1 (MQ=255)
                                                                                                                     gCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTtc                               <  8:183176/148‑1 (MQ=255)
                                                                                                                         ccGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGaa                          <  8:204839/149‑1 (MQ=255)
                                                                                                                                   gcgcCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAACAGGCCCTg                 <  8:331809/148‑1 (MQ=255)
                                                                                                                                     gccggccgCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCt                                                                  >  1:238964/1‑97 (MQ=255)
                                                                                                                                     gccggccgCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCt                                                                  <  2:238964/97‑1 (MQ=255)
                                                                                                                                      tgtaaagCACCCGCTGCAGGCGGCGCTGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCGGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAACAGGCCCTggg               <  8:5913/141‑1 (MQ=255)
                                                                                                                                       cggccgCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAACAGGCCCTGGGCCt            <  6:262177/149‑1 (MQ=255)
                                                                                                                                       cggccgCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAACAGGCCCTGGGCCt            <  1:116905/149‑1 (MQ=255)
                                                                                                                                             cACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAAc                         <  3:208137/130‑1 (MQ=255)
                                                                                                                                             cACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAAc                         >  4:208137/1‑130 (MQ=255)
                                                                                                                                              aCCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAACAGGCCCTGGGCCTGCCACGt     <  2:326840/149‑1 (MQ=255)
                                                                                                                                                 cGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAACAGGCCCTGGGCCTGCCACGTGcc  >  4:225664/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTCAGCGTCAGGCGCGGCGAGATCTGCGCCCTGATCGGGCCGAATGGCGCCGGCAAGAGCTCGCTGCTGAACATCCTCAATGGCGTCTACCGGGCCGATGCCGGCCAGCTGTTCTTCGCCACCGAGCCGTTGCGCCGGCCGCACCCGCTGCAGGCGGCGCGGCTGGGGATCGGCCGTACCTTCCAGAACAATGCGCTGTTCAAGAAGATGAGTGTGCTGGACAACCTGCTCACCGGTCTGTCGCGTTTCCAGCGCAGTTACTTTCTCGAACAGGCCCTGGGCCTGCCACGTGCC  >  NC_002947/3131224‑3131517

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: