Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,982,924 A→G 15.4% Q278R (CAG→CGG)  PP_1775 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,982,9240AG15.4% 51.6 / 4.2 26Q278R (CAG→CGG) PP_1775metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base A (11/11);  new base G (2/2);  total (13/13)
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.41e-02

CGACAAGGGCACCCTCCTGGTTCCTGCATTCAGTCTTGGTCGCACGCAGGAGCTGTTGTACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGG  >  NC_002947/1982781‑1983044
                                                                                                                                               |                                                                                                                        
cGACAAGGGCACCCTCCTGGTTCCTGCATTAAGTCTTGGTCGCACGCAGGAGCTGTTGTACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGcc                                                                                                                    <  1:138450/150‑1 (MQ=255)
                      ccTGCATTCAGTCTTGGTCGCACGCAGGAGCTGTTGTACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTg                                                                                              <  1:17348/150‑1 (MQ=255)
                       cTGCATTCAGTCTTGGTCGCACGCAGGAGCTGTTGTACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTGCCCCTTgg                                                                                             >  1:218193/1‑150 (MQ=255)
                            ttCAGTCTTGGTCGCACGCAGGAGCTGTTGTACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAg                                                                                        >  1:509270/1‑150 (MQ=255)
                             tCAGTCTTGGTCGCACGCAGGAGCTGTTGTACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAgc                                                                                       >  2:249708/1‑150 (MQ=255)
                                              caGGAGCTGTTGTACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTAc                                                                      <  2:212576/150‑1 (MQ=255)
                                                      gttgtACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCAGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCg                                                                                                    >  1:186709/1‑112 (MQ=255)
                                                      gttgtACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCAGACGGTGAGCCAACCGACTGGTCACAGTTGCCCATCATCCTCGATTCg                                                                                                    <  2:186709/112‑1 (MQ=255)
                                                          tACGAGCTGGAAGACACACTGCATCCCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCAc                                                          <  2:144121/150‑1 (MQ=255)
                                                            cGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGa                                                        <  2:255743/150‑1 (MQ=255)
                                                                cTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTTCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAatat                                                    >  1:599910/1‑150 (MQ=255)
                                                                        cATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGc                                            <  1:249708/150‑1 (MQ=255)
                                                                         aTTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAAGCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCa                                           >  2:209353/1‑150 (MQ=255)
                                                                                    cAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCCGCACGAATCTAAGAATGAAGACACCCGGgc                                >  2:35780/1‑150 (MQ=255)
                                                                                              ctgAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTgca                      <  1:328387/150‑1 (MQ=255)
                                                                                               tgAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTgcag                     >  1:458119/1‑150 (MQ=255)
                                                                                                 aaCACTACCGGCCCGGCTCAGGAGGGTGATCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGTATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTgcagca                   <  2:509270/150‑1 (MQ=255)
                                                                                                         cggcccggcTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGa           >  1:187155/1‑150 (MQ=255)
                                                                                                                  tCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCGTCATCCTCGATTCGCCCTCGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGGGGGTTGCAGCAGGGCGCGATCCACTCgg  >  1:276145/1‑150 (MQ=255)
                                                                                                                  tCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCgg  <  2:187155/150‑1 (MQ=255)
                                                                                                                    cGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGAt          >  2:77227/1‑140 (MQ=255)
                                                                                                                    cGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGAt          <  1:77227/140‑1 (MQ=255)
                                                                                                                                      aCTGGTCACGGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTgcagc                    >  2:428534/1‑112 (MQ=255)
                                                                                                                                      aCTGGTCACGGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTgcagc                    >  2:428587/1‑112 (MQ=255)
                                                                                                                                      aCTGGTCACGGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTgcagc                    <  1:428587/112‑1 (MQ=255)
                                                                                                                                      aCTGGTCACGGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTgcagc                    <  1:428534/112‑1 (MQ=255)
                                                                                                                                               |                                                                                                                        
CGACAAGGGCACCCTCCTGGTTCCTGCATTCAGTCTTGGTCGCACGCAGGAGCTGTTGTACGAGCTGGAAGACATTCTGCACCGCAAAGCCCTGCTGAACACTACCGGCCCGGCTCCGGACGGTGAGCCAATCGACTGGTCACAGTTGCCCATCATCCTCGATTCGCCCTTGGCTCAGCGCATCACCCGCGTGTACCGGGAGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGG  >  NC_002947/1982781‑1983044

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: