Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,503,197 T→G 22.4% P33P (CCT→CCG) ‡ oprC → copper receptor OprC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,503,1970TG22.4% 18.5 / 5.2 18P33P (CCT→CCG) ‡oprCcopper receptor OprC
Reads supporting (aligned to +/- strand):  ref base T (4/10);  new base G (2/2);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

ACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCG  >  NC_002947/5503048‑5503330
                                                                                                                                                     |                                                                                                                                     
aCGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCt                                                                                                                                       >  2:19108/1‑150 (MQ=255)
  gAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTgg                                                                                                                                     <  2:192637/150‑1 (MQ=255)
     gCCAGGAGGCTTCGCGCTAACTCATGACTGATCGCCTGGAATCATTATGTCCGGCTGCACCCCTTCTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCGCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCa                                                                                                                                  <  1:186786/150‑1 (MQ=255)
      ccAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTCTCCCTGCGTGGCACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCACGGCCAGGGCCGGCGGCCAGGGCCAt                                                                                                                                 <  2:46443/150‑1 (MQ=255)
           aggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGa                                                                                                                            <  1:370191/150‑1 (MQ=255)
                                                           cTGCACCCCGGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGt                                                                            <  1:19108/150‑1 (MQ=255)
                                                                         ttGCCTTGTCCCTGCGTGGGACCACCGCCGCGCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGTGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgc                                                              <  2:59872/150‑1 (MQ=255)
                                                                              ttGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCATGGCCCCGAGTTCGCCGCTg                                                         >  1:171846/1‑150 (MQ=255)
                                                                                    cTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCTGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACtgtg                                                   <  2:45390/150‑1 (MQ=255)
                                                                                        gTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCa                                               <  1:145084/150‑1 (MQ=255)
                                                                                                       cccTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAcc                                >  1:9349/1‑150 (MQ=255)
                                                                                                          tgtgCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCcgc                             <  1:229755/150‑1 (MQ=255)
                                                                                                              cGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCgcca                         <  2:9349/150‑1 (MQ=255)
                                                                                                                         tCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCt                                                          <  2:140947/106‑1 (MQ=255)
                                                                                                                         tCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCt                                                          >  1:140947/1‑106 (MQ=255)
                                                                                                                               tcATGGCCAGGGCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAg          <  1:84064/147‑1 (MQ=255)
                                                                                                                               tcATGGCCAGGGCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAg          >  2:84064/2‑148 (MQ=255)
                                                                                                                                cATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGa       >  1:35561/1‑150 (MQ=255)
                                                                                                                                     ccagggCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgcg  >  1:262629/5‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                     
ACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCG  >  NC_002947/5503048‑5503330

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