Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,141,770 A→C 8.9% M245L (ATG→CTG) ‡ PP_3645 → L‑proline amide hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,141,7700AC 8.9% 131.5 / 5.2 45M245L (ATG→CTG) ‡PP_3645L‑proline amide hydrolase
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base C (2/2);  total (24/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 = 4.81e-01

ACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGC  >  NC_002947/4141624‑4141912
                                                                                                                                                  |                                                                                                                                              
aCCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCatat                                                                                                                                               >  1:299701/1‑148 (MQ=255)
            acCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGaa                                                                                                                                  <  5:58572/149‑1 (MQ=255)
                      tATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGCCCCGCCAGCTCGCCGCCAGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGt                                                                                                                         >  6:317651/1‑148 (MQ=255)
                        tCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGaa                                                                                                                      <  7:159707/149‑1 (MQ=255)
                          cAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGt                                                                                                                    <  5:265010/149‑1 (MQ=255)
                           aGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGt                                                                                                                    >  3:12467/1‑148 (MQ=255)
                             gTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATCATATGCCGTTgct                                                                                                                                   <  1:209568/131‑3 (MQ=255)
                             gTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATCATATGCCGTTgct                                                                                                                                   >  2:209568/1‑129 (MQ=255)
                                ccGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTcaac                                                                                                                >  6:116079/1‑147 (MQ=255)
                                  gACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTc                                                                                                                   >  7:294784/1‑142 (MQ=255)
                                  gACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTc                                                                                                                   <  8:294784/142‑1 (MQ=255)
                                   aCGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGAt                                                                                                           >  5:42763/1‑149 (MQ=255)
                                       tGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCACCGGCATCATATGCCGTTgctc                                                                                                                                  <  6:252454/122‑4 (MQ=255)
                                       tGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCACCGGCATCATATGCCGTTgctc                                                                                                                                  >  5:252454/1‑119 (MQ=255)
                                          tCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGctgctg                                                                                                    <  4:101501/149‑1 (MQ=255)
                                           cGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGctgctg                                                                                                    >  4:130813/1‑148 (MQ=255)
                                                  ggCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCCGt                                                                                                                                      >  8:298567/1‑107 (MQ=255)
                                                   gCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCt                                                                                            >  6:16117/1‑148 (MQ=255)
                                                   gCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTc                                                                                           <  4:12467/149‑1 (MQ=255)
                                                         tcgactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGctgctg                                                                                                    <  3:187585/134‑1 (MQ=255)
                                                         tcgactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGctgctg                                                                                                    >  4:187585/1‑134 (MQ=255)
                                                         tcgactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCAccc                                                                                       >  6:63810/1‑147 (MQ=255)
                                                           gactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGT‑GAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTc                                                                                              >  7:148033/1‑137 (MQ=255)
                                                           gactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGT‑GAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTc                                                                                              <  8:148033/137‑1 (MQ=255)
                                                             ctcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGccc                                                                                 <  3:130813/149‑1 (MQ=255)
                                                              tcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCg                                                                                >  3:223116/1‑149 (MQ=255)
                                                              tcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCg                                                                                >  4:181387/1‑149 (MQ=255)
                                                                 ggCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCAc                                                                             >  6:37650/1‑149 (MQ=255)
                                                                  gCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCAcc                                                                            <  5:116079/149‑1 (MQ=255)
                                                                      cACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCa                                                                        <  5:16117/149‑1 (MQ=255)
                                                                          ccGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCg                                                                     >  7:186669/1‑148 (MQ=255)
                                                                                    cccGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAg                                                          <  5:317651/149‑1 (MQ=255)
                                                                                    cccGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAg                                                          <  3:181387/149‑1 (MQ=255)
                                                                                     ccGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAgg                                                         >  7:172221/1‑149 (MQ=255)
                                                                                                ccgccTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCa                                              <  4:172109/149‑1 (MQ=255)
                                                                                                            cTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCtcgttc                                  >  6:270566/1‑149 (MQ=255)
                                                                                                                                cTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGAc              <  5:63810/149‑1 (MQ=255)
                                                                                                                                cTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGAc              >  4:7781/1‑149 (MQ=255)
                                                                                                                                cTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGAc              >  1:107927/1‑149 (MQ=255)
                                                                                                                                 tCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACg             >  5:299055/1‑149 (MQ=255)
                                                                                                                                        gCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATg      <  4:223116/149‑1 (MQ=255)
                                                                                                                                         cAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATgg     >  1:103396/1‑149 (MQ=255)
                                                                                                                                          aaCGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGgc    <  5:37650/149‑1 (MQ=255)
                                                                                                                                            cGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGgcgc  <  6:42763/149‑1 (MQ=255)
                                                                                                                                             ggCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGc                                               >  7:202818/1‑103 (MQ=255)
                                                                                                                                             ggCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGc                                               <  8:202818/103‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
ACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGC  >  NC_002947/4141624‑4141912

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