Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,890,184 T→C 9.5% Y595H (TAT→CAT)  PP_5164 → penicillin amidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,890,1840TC 9.5% 117.1 / 3.1 42Y595H (TAT→CAT) PP_5164penicillin amidase family protein
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base C (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 = 7.82e-01

TGGCACCGCTGCTAGCGACGTTGCGTCAGATCGCCCAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAACACCCGCGGCGCCAACCAGCGCACTGACCGCAG  >  NC_002947/5890038‑5890332
                                                                                                                                                  |                                                                                                                                                     
tGGCACCGCTGCTAGCGACGTTGCGTCAGATCGCCCAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATg                                                                                                                                                    <  2:146617/150‑1 (MQ=255)
    aCCGCTGCTAGCGACGTTGCGTCAGATCGCCCAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCcgc                                                                                                                                               >  2:838544/1‑150 (MQ=255)
                        gTCAGATCGCCCAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCg                                                                                                                           <  1:441943/150‑1 (MQ=255)
                             aTCGCCCAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCCATGCCGCGTTGCATGAc                                                                                                                                     >  2:657582/1‑135 (MQ=255)
                             aTCGCCCAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCCATGCCGCGTTGCATGAc                                                                                                                                     <  1:657582/135‑1 (MQ=255)
                               cGCCCAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGAGCGATATCta                                                                                                                    >  1:441473/1‑149 (MQ=255)
                                 cccAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGt                                                                                                                  <  1:784268/150‑1 (MQ=255)
                                 cccAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACAAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGt                                                                                                                  >  2:422039/1‑150 (MQ=255)
                                  ccAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGtt                                                                                                                 <  2:18051/150‑1 (MQ=255)
                                   cAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTc                                                                                                                >  2:196175/1‑150 (MQ=255)
                                        tGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGt                                                                                                           >  2:337036/1‑150 (MQ=255)
                                               gaagaggaaCTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGAt                                                                                                     >  2:374258/1‑149 (MQ=255)
                                               gaagaggaaCTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGAt                                                                                                     <  1:374258/149‑1 (MQ=255)
                                                     gaaCTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCAcc                                                                                              >  2:837246/1‑150 (MQ=255)
                                                            tCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGAt                                                                                                     <  1:254431/136‑1 (MQ=255)
                                                            tCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGAt                                                                                                     >  2:254431/1‑136 (MQ=255)
                                                                  agcTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCAccgggaa                                                                                         <  1:565028/142‑6 (MQ=255)
                                                                  agcTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCAccgggaa                                                                                         >  2:565028/1‑137 (MQ=255)
                                                                         tgCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGAtgctgc                                                                             >  2:346416/1‑147 (MQ=255)
                                                                         tgCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGAtgctgc                                                                             <  1:346416/147‑1 (MQ=255)
                                                                                 gcggcgACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGcc                                                                  <  2:272493/150‑1 (MQ=255)
                                                                                  cggcgACTATCCACTGGACTCAACCGGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCAGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGGGCTGCCGCGCCt                                                                 <  1:116820/150‑1 (MQ=255)
                                                                                   ggcgAATATCCTCTGAACTCAACCAGCGCCACGCTGTTCAACCAGTTCCACTGCCGGCGCCTTGCGCAGCCGCTGCGGAACGAATTGTGCGATAACAGGTTCCCGGTGATGTTCAGCGGCCGTGCCATCGATGCTGCGCTGCCGCGCCTg                                                                <  1:196175/150‑1 (MQ=255)
                                                                                    gcgACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTgg                                                               >  2:540117/1‑150 (MQ=255)
                                                                                        cTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCgcgc                                                                   >  2:429110/1‑142 (MQ=255)
                                                                                        cTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCgcgc                                                                   <  1:429110/142‑1 (MQ=255)
                                                                                             cACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATg                                                      >  2:801839/1‑150 (MQ=255)
                                                                                              aCTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGc                                                     <  1:797147/150‑1 (MQ=255)
                                                                                                   aCTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGt                                                <  1:540117/150‑1 (MQ=255)
                                                                                                                cACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCATGGACGCCGGTGCCGAGTCGCCCTGGTGGaa                                   >  2:623438/1‑150 (MQ=255)
                                                                                                                   gCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAacac                                >  1:159582/1‑150 (MQ=255)
                                                                                                                          aaCCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAACACCCgcggc                         >  2:443661/1‑150 (MQ=255)
                                                                                                                            ccAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAACACCCGCGgcgc                       >  1:591494/1‑150 (MQ=255)
                                                                                                                               gTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCtggt                                       >  2:423693/1‑131 (MQ=255)
                                                                                                                               gTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCtggt                                       <  1:423693/131‑1 (MQ=255)
                                                                                                                                   ctctACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTCGCGGAACAGCAGCAGCGCCaccaag                >  2:819964/1‑145 (MQ=255)
                                                                                                                                    tctACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAACACCCGCGGCGCCAACCAgc               <  1:838544/150‑1 (MQ=255)
                                                                                                                                      tACGAGCTGGCCCATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCtggtg                                      <  1:759816/125‑1 (MQ=255)
                                                                                                                                      tACGAGCTGGCCCATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCtggtg                                      >  2:759816/1‑125 (MQ=255)
                                                                                                                                           gCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCCCCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCAGACCATGCCGAGTCGCCCTGGTGGAACACCCGCGGCGCCAACCAGCGCACTGa        <  1:662473/150‑1 (MQ=255)
                                                                                                                                            cTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAACACCCGCGGCGCCAACCAGCGCACTGAc       <  2:343888/150‑1 (MQ=255)
                                                                                                                                                 cTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAACACCCGCGGCGCCAACCAGCGCACTGACCGCAg  <  1:837246/150‑1 (MQ=255)
                                                                                                                                                  tATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAACACCcgc                            <  2:452383/123‑1 (MQ=255)
                                                                                                                                                  tATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAACACCcgc                            >  1:452383/1‑123 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGGCACCGCTGCTAGCGACGTTGCGTCAGATCGCCCAGGGTGACGAGGAAGAGGAACTGGTCGAGCAGCTGGCTGCCTGGCGCGGCGACTATCCACTGGACTCAACCAGCGCCACGCTGTTCAACCAGTTCCTCTACGAGCTGGCCTATGCCGCGTTGCATGACGAATTGGGCGATATCTGGTTCCCGGTGCTGATCAGCACCCGTGCCATCGATGCTGCGCTGCCGCGCCTGGCCGCCGATGCCGAGTCGCCCTGGTGGAACACCCGCGGCGCCAACCAGCGCACTGACCGCAG  >  NC_002947/5890038‑5890332

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