Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,712,080 A→T 19.8% intergenic (+71/‑187) ppc → / → adk phosphoenolpyruvate carboxylase/adenylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,712,0800AT19.8% 69.7 / 6.8 38intergenic (+71/‑187)ppc/adkphosphoenolpyruvate carboxylase/adenylate kinase
Reads supporting (aligned to +/- strand):  ref base A (24/5);  new base T (4/4);  total (28/10)
Fisher's exact test for biased strand distribution p-value = 7.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCC  >  NC_002947/1711933‑1712226
                                                                                                                                                     |                                                                                                                                                   
ctctGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAg                                                                                                                                                    >  5:41136/1‑149 (MQ=255)
          cccGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑c                                                                                                                                         >  8:278590/1‑149 (MQ=255)
            cGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGCAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCt                                                                                                                                       <  1:175924/149‑1 (MQ=255)
                    cAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTc                                                                                                                               >  7:242565/1‑149 (MQ=255)
                                             gTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCaa                                                                                                      >  7:3535/1‑149 (MQ=255)
                                                   ccgccgGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTc                                                                                                                               >  5:233203/1‑118 (MQ=255)
                                                    cgccgGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑TCGCCACGCGCAACCC‑CCTGCCCGGTCAGCCGCAAGGGGCGCAATCCGGCCACTTTggg                                                                                               >  8:253254/1‑149 (MQ=255)
                                                    cgccgGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCCGTCAGCCGCAAGTGGCGCATTCCTGCAACttttg                                                                                                >  6:71008/1‑146 (MQ=255)
                                                     gccgGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGa                                                                                              >  2:207826/1‑149 (MQ=255)
                                                        ggggTGCGCAAAACCGGTTATGGGCACGCCGGGGTAATGTCTGGTGACGGCCATACCGGCCCGGC‑AAAAGGAGGGGGTGGGGGGTT‑GGGCCTGGCGAAACGG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGc                                                                                           <  2:24954/149‑1 (MQ=255)
                                                                      cGGTTATGGGCATGCCGGGGTCTTGTGTGGTGAGGGCAATACAGGCCCGGC‑AAATGGAGGGGGTGGGGGGTT‑GCGCCTGGCGAAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCtt                                                                                         <  2:235467/137‑1 (MQ=255)
                                                                      cGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACGAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCtt                                                                                         >  1:235467/1‑137 (MQ=255)
                                                                        gCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTg                                                                                                 >  1:262815/1‑127 (MQ=255)
                                                                            aGGGCACGCCGGGGTCTTGTCTGGTCCCGGCCAAACCGCCCCGCC‑ACATGGAGGGGGAGGGCGGTG‑GCGCCTGGCCCACCCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCt                                                                       <  7:278590/149‑1 (MQ=255)
                                                                            aGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑AAAACGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTAAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCATGTATGCCTGCCGGGCGCt                                                                       >  8:159486/1‑149 (MQ=255)
                                                                              ggCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCtgt                                                                     >  7:273385/1‑149 (MQ=255)
                                                                                   gCCGGGGCCCTGGCTGGTGCCGCCCAT‑CCCGCCCGGCAACATGGA‑GGGGCGGGCGGTTGGCGCCCGCCCCCACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATc                                                                <  3:180023/149‑1 (MQ=255)
                                                                                                cTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTc                                                   >  8:31270/1‑149 (MQ=255)
                                                                                                              aCAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCGGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATgg                                     >  5:262156/1‑149 (MQ=255)
                                                                                                                    ccGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATAc                               >  7:274948/1‑149 (MQ=255)
                                                                                                                            aaGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGAttt                        >  8:36180/1‑149 (MQ=255)
                                                                                                                              ggACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTc                      >  2:198054/1‑149 (MQ=255)
                                                                                                                               tccGGGGCGGTCGGTG‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTcc                     <  4:182580/147‑1 (MQ=255)
                                                                                                                                       gggCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTgg             >  3:226538/1‑149 (MQ=255)
                                                                                                                                        gggggTT‑GCGCCAGGCGAAACGG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTgg             <  5:71008/145‑1 (MQ=255)
                                                                                                                                         gCGGTT‑GCGCCTGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGcc           >  7:177744/1‑149 (MQ=255)
                                                                                                                                          cGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGccc          >  4:173945/1‑149 (MQ=255)
                                                                                                                                           ggTT‑GCGCCTGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCt         >  3:128493/1‑149 (MQ=255)
                                                                                                                                           ggTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTg              >  1:162435/1‑144 (MQ=255)
                                                                                                                                           ggTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTg              <  2:162435/144‑1 (MQ=255)
                                                                                                                                            gTT‑GCGCCTGGCGCAACCGCCCTCCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCAGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCt         >  5:265652/1‑149 (MQ=255)
                                                                                                                                            gTT‑GCGCCTGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTg        >  5:86300/1‑149 (MQ=255)
                                                                                                                                            gTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTg        >  1:142926/1‑149 (MQ=255)
                                                                                                                                                gcgcCTGGCCCACCCC‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTc     <  6:41136/149‑1 (MQ=255)
                                                                                                                                                 cgcCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTc     <  2:142926/148‑1 (MQ=255)
                                                                                                                                                  gcCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCtt                                                    >  2:130868/1‑100 (MQ=255)
                                                                                                                                                  gcCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGc   >  2:154138/1‑149 (MQ=255)
                                                                                                                                                    cAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGcc  >  5:304954/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGC‑ACAAGGACGGGGAGGGCGGTT‑GCGCCAGGCGCAACCG‑CCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCC  >  NC_002947/1711933‑1712226

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