Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,932,047 T→G 9.1% E477A (GAG→GCG)  PP_3466 ← ABC transporter permease/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,932,0470TG 9.1% 119.6 / 2.9 42E477A (GAG→GCG) PP_3466ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base G (2/2);  total (21/21)
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.69e-01

GCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACGGAACAACAC  >  NC_002947/3931909‑3932192
                                                                                                                                                   |                                                                                                                                                 
gCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGGCGAGCGCGGCCAGCag                                                                                                                                         <  8:8711/149‑2 (MQ=255)
   gcgATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGc                                                                                                                                       >  6:164904/1‑148 (MQ=255)
       tCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATcc                                                                                                                                  <  3:182438/149‑1 (MQ=255)
         atttCACGCTCGGGGGGCTGATCCACATCAGAAGTTGGCTCGTCCAGCACCAGGACGATCCCCCCGCCGGGGATATGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCa                                                                                                                                          <  1:98553/146‑1 (MQ=255)
          cTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAgg                                                                                                                               >  3:268219/1‑149 (MQ=255)
                     gTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAg                                                                                                                    >  6:321152/1‑149 (MQ=255)
                      tGGCCTCATCCACTGCCGAAGTGGCCACGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCacg                                                                                                                   >  6:283973/1‑147 (MQ=255)
                                 aCTGCCGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACATCTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAg                                                                                                        <  1:246132/149‑1 (MQ=255)
                                   tGCCGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGc                                                                                                      >  2:23293/1‑149 (MQ=255)
                                     ccGAAGTGGCCTCTTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCATCAAGCCATCCAGGCCCAGCGGCAGTCTGTCGACCAGGCt                                                                                                     <  5:283973/148‑1 (MQ=255)
                                      cGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGAAAAACTGAACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGt                                                                                                   >  5:9431/1‑149 (MQ=255)
                                                       gCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGgcg                                                                                   >  6:115869/1‑148 (MQ=255)
                                                         aCCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTc                                                                                                              >  5:237804/1‑119 (MQ=255)
                                                         aCCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTc                                                                                                              <  6:237804/119‑1 (MQ=255)
                                                         aCCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGcac                                                                                <  7:86413/149‑1 (MQ=255)
                                                          ccAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCGGCGCTGCTCGGGCAGTTGTCATGCGGAGGTTCTGCGACCCCGACAACTGGAGGAAGCGAGCGCCGCGCAAGACCTGCCTGCACAGAGGCAGGCTGTCGTCAAGGCTGTCTCGCCGCGCCAGGCGcag                                                                                >  4:192588/1‑147 (MQ=255)
                                                           cAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGcacca                                                                              >  8:67295/1‑149 (MQ=255)
                                                           cAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGcacca                                                                              >  2:39224/1‑149 (MQ=255)
                                                            aG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGcaccac                                                                             <  1:23293/149‑1 (MQ=255)
                                                                            cGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCa                                                                     <  8:273688/149‑1 (MQ=255)
                                                                              ggg‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGt                                                                   <  5:321152/149‑1 (MQ=255)
                                                                                g‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTtc                                                                 <  5:115869/149‑1 (MQ=255)
                                                                                                         gCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGc                                          >  7:259759/1‑148 (MQ=255)
                                                                                                           gATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGTAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCAt                                       <  5:164904/149‑1 (MQ=255)
                                                                                                              gCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATaa                                     <  8:259759/148‑1 (MQ=255)
                                                                                                                 ttgttTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGc                                   >  6:158166/1‑147 (MQ=255)
                                                                                                                 ttgttTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCAACCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGc                                   >  1:231054/1‑147 (MQ=255)
                                                                                                                   gttTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCaag                                                                                     >  1:131446/1‑93 (MQ=255)
                                                                                                                   gttTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCaag                                                                                     <  2:131446/95‑3 (MQ=255)
                                                                                                                    tttCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAgg                               <  6:64882/148‑1 (MQ=255)
                                                                                                                            ccccGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTgcgc                       >  5:142035/1‑148 (MQ=255)
                                                                                                                              ccGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAg                    <  6:142035/149‑1 (MQ=255)
                                                                                                                                  aCAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTg                >  4:169137/1‑149 (MQ=255)
                                                                                                                                       tgctgGGCGAGCGCGGCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGc                                              <  3:153577/114‑1 (MQ=255)
                                                                                                                                       tgctgGGCGAGCGCGGCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGc                                              >  4:153577/1‑114 (MQ=255)
                                                                                                                                         ctgGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACgg          <  6:9431/148‑1 (MQ=255)
                                                                                                                                          tgGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATaa                                     <  4:10879/120‑1 (MQ=255)
                                                                                                                                          tgGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATaa                                     >  3:10879/1‑120 (MQ=255)
                                                                                                                                             ccGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACGGaacaa     <  5:158166/149‑1 (MQ=255)
                                                                                                                                              cgcgCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTc                            >  6:83156/1‑125 (MQ=255)
                                                                                                                                              cgcgCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTc                            <  5:83156/125‑1 (MQ=255)
                                                                                                                                                 gcTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCGCGGAACAacac  >  6:77906/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAG‑‑‑GAT‑‑‑‑‑CGCCGGGG‑CCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACGGAACAACAC  >  NC_002947/3931909‑3932192

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