Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,719,745 T→G 8.7% Q17P (CAG→CCG)  PP_4176 ← 5‑oxo‑L‑prolinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,719,7450TG 8.7% 131.3 / 5.1 46Q17P (CAG→CCG) PP_41765‑oxo‑L‑prolinase
Reads supporting (aligned to +/- strand):  ref base T (23/19);  new base G (2/2);  total (25/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 = 6.61e-01

GCCACCGGTATGCCTGAAATCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCTATGCTTGAGGCACTTCCAC  >  NC_002947/4719598‑4719892
                                                                                                                                                     |                                                                                                                                                   
gCCACCGGTATGCCTGAAATCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCt                                                                                                                                                     <  2:65890/148‑1 (MQ=255)
gCCACCGGTATGCCTGAAATCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTg                                                                                                                                                    >  2:347900/1‑149 (MQ=255)
    ccGGTATGCCTGAAATCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGca                                                                                                                                                 >  3:121874/1‑148 (MQ=255)
     cGGTATGCCTGAAATCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGcaca                                                                                                                                               <  8:210703/149‑1 (MQ=255)
            ccTGAAATCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACagccagcc                                                                                                                                        <  2:311155/149‑1 (MQ=255)
              tGAAATCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCg                                                                                                                                       >  1:108129/1‑148 (MQ=255)
               gAAATCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGaa                                                                                                                                     >  5:314584/1‑149 (MQ=255)
                aaaTCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGa                                                                                                                                      <  2:266441/147‑1 (MQ=255)
                   tcatcGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGaacaac                                                                                                                                 <  4:121874/149‑1 (MQ=255)
                     atcGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAacaaca                                                                                                                                <  2:108129/148‑1 (MQ=255)
                      tcGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGaacaa                                                                                                                                  >  5:291059/1‑145 (MQ=255)
                      tcGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGaacaa                                                                                                                                  <  6:291059/145‑1 (MQ=255)
                        gACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGcac                                                                                                                            <  6:314584/149‑1 (MQ=255)
                             cTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGTGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAg                                                                                                                       >  2:116123/1‑149 (MQ=255)
                              tCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGc                                                                                                                      >  3:148172/1‑149 (MQ=255)
                                acccccTGCAGCATCCATGGCCTTGCTCAGCGCCGGCATGCCTTCAGTTTCCTGGAAGAAGTCCACGTAGTGCAGGTGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACCGCCCGCCGAACCACAGCACACACGc                                                                                                                      >  2:241517/3‑147 (MQ=255)
                                        cagcaTCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTg                                                                                                            >  7:221777/1‑149 (MQ=255)
                                          gcaTCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGc                                                                                                          >  3:357965/1‑149 (MQ=255)
                                                          tCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCAtcctc                                                                                          >  6:284013/1‑149 (MQ=255)
                                                             gcgcCGTCATGCCTTCACTTTCTTGGAAGAAGTCCACTTAGTGAAGATGGTTCAACACTGT‑TCCCTCCACTGCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATcctcct                                                                                        <  1:116123/149‑1 (MQ=255)
                                                                 cGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATg                                                                                    >  2:347558/1‑148 (MQ=255)
                                                                   gCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCaa                                                                                 <  1:347900/149‑1 (MQ=255)
                                                                        ccTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTAc                                                                             >  1:218605/1‑148 (MQ=255)
                                                                             aCTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAAcc                                                                        >  3:21147/1‑148 (MQ=255)
                                                                                    tGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGAc                                                                >  7:28964/1‑149 (MQ=255)
                                                                                      gaagaaGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGCACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGAc                                                                >  5:20842/1‑147 (MQ=255)
                                                                                         gaagTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGCACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAggc                                                           <  6:20842/149‑1 (MQ=255)
                                                                                             tCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGgcc                                                       <  4:148172/149‑1 (MQ=255)
                                                                                                aCGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGgccgcc                                                    <  4:21147/149‑1 (MQ=255)
                                                                                                  gTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCGGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGgccgcc                                                    >  2:300456/1‑147 (MQ=255)
                                                                                                  gTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCGGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCaa                                                  >  8:334942/1‑149 (MQ=255)
                                                                                                     gTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGa                                                <  1:347558/148‑1 (MQ=255)
                                                                                                        cAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGt                                            >  8:328900/1‑149 (MQ=255)
                                                                                                         aGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGtt                                           >  2:175499/1‑149 (MQ=255)
                                                                                                                      cTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCAcggcgg                               <  4:357965/149‑1 (MQ=255)
                                                                                                                           gCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCaa                           <  2:218605/148‑1 (MQ=255)
                                                                                                                              gTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCg                       >  4:111878/1‑149 (MQ=255)
                                                                                                                              gTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCAGGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCg                       >  7:192628/1‑149 (MQ=255)
                                                                                                                               tACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGc                      <  8:28964/149‑1 (MQ=255)
                                                                                                                                 cTC‑GCGCGCCTGGGCCGCCGGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCt                     <  1:300456/148‑1 (MQ=255)
                                                                                                                                 cTC‑GCGCGCCTGGGCCGCCGGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCt                     <  7:334942/148‑1 (MQ=255)
                                                                                                                                     gcgcgcCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCTATGc                 >  2:22169/1‑149 (MQ=255)
                                                                                                                                           cTGTGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCTATGCTTGAgg           <  8:221777/149‑1 (MQ=255)
                                                                                                                                           cTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCTATGCTTGAgg           <  3:111878/149‑1 (MQ=255)
                                                                                                                                                ccgccTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCTATGCTTGAGGCACtt      >  3:105272/1‑149 (MQ=255)
                                                                                                                                                 cgccTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCTATGCTTGAGGCACTTc     >  7:136637/1‑149 (MQ=255)
                                                                                                                                                    cTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAATTACGAACCCAGTAGACCAGGCGGCAGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCTATGCTTGAGGCACTTCCAc  <  4:224568/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCACCGGTATGCCTGAAATCATCGACTGCTCAACCCCTGCAGCATCCATGGCCTTGATCAGCGCCGGCATGCCTTCACTTTCCTGGAAGAAGTCCACGTAGTGCAGATGG‑GCATCACTGTAGCGGTACTC‑GCGCGCCTGGGCCGCCTGGCACAGCCAGCCGAACAACAGCACACAAGCCAGCGCTCTGGCGGTCATGGTCATCCTCCTCATGCAAGTACGAACCCAGTAGACCAGGCGGCCGCCAAGACGGTTCAGCCACGGCGGGCAAACCGCTATGCTTGAGGCACTTCCAC  >  NC_002947/4719598‑4719892

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