Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,353,340 C→G 15.8% P141R (CCG→CGG)  yrpB → 2‑nitropropane dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,353,3400CG15.8% 94.5 / 12.4 38P141R (CCG→CGG) yrpB2‑nitropropane dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (14/18);  new base G (3/3);  total (17/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 = 8.54e-01

AGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCA  >  NC_002947/4353196‑4353486
                                                                                                                                                |                                                                                                                                                  
aGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAgg                                                                                                                                               <  1:897106/150‑1 (MQ=255)
aGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAgg                                                                                                                                               >  1:81258/1‑150 (MQ=255)
  gCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGcc                                                                                                                                             >  1:970459/1‑150 (MQ=255)
  gCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGcc                                                                                                                                             <  1:173320/150‑1 (MQ=255)
                 taCAGCGAAGTTGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                              <  2:585898/150‑1 (MQ=255)
                  aCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGCGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCGGCAGGCCGaaa                                                                                                                                         <  2:773881/138‑2 (MQ=255)
                  aCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGCGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCGGCAGGCCGaaa                                                                                                                                         >  1:773881/1‑137 (MQ=255)
                    aGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTa                                                                                                                           >  1:1128951/1‑150 (MQ=255)
                      cGAAGTGGGCGCCGATTTCACGGCGCCCGCGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTaaa                                                                                                                         <  1:782426/150‑1 (MQ=255)
                       gAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAg                                                                                                                        <  1:909186/150‑1 (MQ=255)
                                gcCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGt                                                                                                               <  1:91508/150‑1 (MQ=255)
                                   gATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGcc                                                                                                            >  2:1039887/1‑150 (MQ=255)
                                         aCGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                      <  1:212840/150‑1 (MQ=255)
                                          cGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTAc                                                                                                     >  2:117675/1‑150 (MQ=255)
                                              gcCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAATGCCAGCGGTGCCAAGGTACTGTc                                                                                                 <  1:240665/150‑1 (MQ=255)
                                                              aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                 <  1:117675/150‑1 (MQ=255)
                                                                   tGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaaga                                                                            >  2:129760/1‑150 (MQ=255)
                                                                              acgaGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctggc                                                                 >  1:657575/1‑150 (MQ=255)
                                                                                          gcctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATg                                                     >  1:724211/1‑150 (MQ=255)
                                                                                           cctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGc                                                                      >  2:14132/1‑132 (MQ=255)
                                                                                           cctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGc                                                                      <  1:14132/132‑1 (MQ=255)
                                                                                             tgctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCt                                                  >  2:258019/1‑150 (MQ=255)
                                                                                              gctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtg                                                 <  2:99571/150‑1 (MQ=255)
                                                                                              gctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtg                                                 >  1:206487/1‑150 (MQ=255)
                                                                                               ctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGTCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgt                                                <  1:527893/150‑2 (MQ=255)
                                                                                               ctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgc                                                <  1:283808/150‑1 (MQ=255)
                                                                                                  gTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgat                                             >  2:763397/1‑150 (MQ=255)
                                                                                                           ttGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGatcatc                                    <  2:206487/150‑1 (MQ=255)
                                                                                                               gcCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACGGCTGAAGCGGTTATATGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCa                                <  2:426670/150‑1 (MQ=255)
                                                                                                                 ccGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCCTg                              >  1:163744/1‑150 (MQ=255)
                                                                                                                  cGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATgg                             <  2:928383/150‑1 (MQ=255)
                                                                                                                     aGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGt                          <  2:81258/150‑1 (MQ=255)
                                                                                                                        tggtgAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATg                       >  1:1130100/1‑150 (MQ=255)
                                                                                                                                       tCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg              <  2:1052322/144‑1 (MQ=255)
                                                                                                                                       tCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg              >  1:1052322/1‑144 (MQ=255)
                                                                                                                                           ccTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGa                                                             >  1:919848/1‑93 (MQ=255)
                                                                                                                                           ccTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGa                                                             <  2:919848/93‑1 (MQ=255)
                                                                                                                                             tGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCTGTGCCAAGGTAGTGTCGAGCGCCACGGCGGGGGAGGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCa  <  2:417441/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCA  >  NC_002947/4353196‑4353486

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