Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,712,080 A→T 19.2% 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.2% 62.0 / 8.6 40intergenic (+71/‑187)ppc/adkphosphoenolpyruvate carboxylase/adenylate kinase
Reads supporting (aligned to +/- strand):  ref base A (20/10);  new base T (2/6);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 4.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

CTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGC  >  NC_002947/1711933‑1712225
                                                                                                                                                     |                                                                                                                                                 
ctctGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCNGCACAAGGACGGGGAGGGCGGTTGCGCCAg                                                                                                                                                  <  1:22436/149‑1 (MQ=255)
        aGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCTGGCGCAAcc                                                                                                                                          <  3:124928/149‑1 (MQ=255)
          cccGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGc                                                                                                                                        >  3:155810/1‑149 (MQ=255)
            cGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCTGGCGCAACCGccc                                                                                                                                      <  2:55413/149‑2 (MQ=255)
             gTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCt                                                                                                                                      >  7:13989/1‑148 (MQ=255)
             gTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCt                                                                                                                                      <  8:13989/148‑1 (MQ=255)
                        ccTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCACAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGATGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGcc                                                                                                                          >  1:203284/1‑149 (MQ=255)
                                               aTTGCCGCCGGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACt                                                                                                   >  7:53437/1‑149 (MQ=255)
                                                   ccgccgGGCTGCGCAACACCGGCTAAGGGGCAAGCCCGGGTCA‑TGGCTGGTGCCGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATACCTGCACCTTTg                                                                                                >  1:153394/1‑149 (MQ=255)
                                                   ccgccgGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTgg                                                                                               >  1:58508/1‑149 (MQ=255)
                                                    cgccgGGCTGCGCAACACCGGCTAA‑GGGCAAGCCGGGGTCA‑TGTCTGGTGACGGCAATACCGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTggg                                                                                              <  1:197471/149‑1 (MQ=255)
                                                    cgccgGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTgg                                                                                               >  3:208706/1‑148 (MQ=255)
                                                          gCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCtt                                                                                        >  3:223144/1‑149 (MQ=255)
                                                                   cacCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTTCGCCAGGCGCAACCGcc                                                                                                                                       >  2:141735/1‑93 (MQ=255)
                                                                   cacCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGcc                                                                                                                                       <  1:141735/93‑1 (MQ=255)
                                                                    acCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGc                                                                              >  6:223959/1‑149 (MQ=255)
                                                                     ccGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGcc                                                                             >  8:229312/1‑149 (MQ=255)
                                                                         cTAA‑GGGCAAGCCCGGGTCA‑TGTCTGATGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGgc                                                                         >  1:143527/1‑149 (MQ=255)
                                                                               tgcgctgccGGGTTATTGTGTGGTCACGGCCATACAGGCCGGGAAAATGTAGGGGTTGGGCGTTGGCGCCAGGCGAAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCt                                                                      <  2:68201/141‑1 (MQ=255)
                                                                                                   tGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCt                                                 >  6:108188/1‑149 (MQ=255)
                                                                                                    ggTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCt                                                                      >  1:182017/1‑127 (MQ=255)
                                                                                                    ggTGACGGCAATACAGGCCGGGCAAAAGGAGGGGGAGGGCGGTTGCGCCAGGCGAAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCt                                                                      <  2:182017/127‑1 (MQ=255)
                                                                                                        agggCCATACTGGCCGGTCACATGCAGGGGGCGGGGGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGAcc                                             <  8:44796/146‑1 (MQ=255)
                                                                                                            ccaacccggCCCGCCACAAGGCGGGGGCGGGGGGTGGCGCCCGGCGCACCCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGAttt                                        <  8:53437/142‑1 (MQ=255)
                                                                                                               tACTGGCCCGGCAAATGGAGGGGTTGGGGGTTTGCGCCTGGCGAAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATg                                     <  2:58508/149‑1 (MQ=255)
                                                                                                                aCAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATgg                                    >  7:13097/1‑149 (MQ=255)
                                                                                                                          cacaAGGACGGGGAGGGCGGTTTCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGgcg                                                                                                                >  1:105855/1‑63 (MQ=255)
                                                                                                                              aGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCACCTTTGGGACGCTTGTCTGGCTGCCGGGCGCGGTGTATCTTGAGCAGCCGTCTGACCGATATATGGTCATACCCGAgttt                      >  5:208530/1‑149 (MQ=255)
                                                                                                                                 tgggggtggggggtggCGCCCGGCCAAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCg                   <  2:210867/147‑1 (MQ=255)
                                                                                                                                 agggggCGGGCGGTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCg                   <  2:153394/147‑1 (MQ=255)
                                                                                                                                   ggggAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCg                   >  1:210867/1‑147 (MQ=255)
                                                                                                                                     ggAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGTTCAGCCGCAAGTGGCGCTTTCGTGCAACTTTGGGACGCTTGTCTGGGTGCCGGGCGCTGTGAATATTGTGCAGCCTTCTGACCGATTTATGGACATACCCGATTTTCCgtaat               >  6:191627/1‑145 (MQ=255)
                                                                                                                                            ggTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCtt                                                   <  1:62712/106‑1 (MQ=255)
                                                                                                                                            ggTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCtt                                                   <  1:58959/106‑1 (MQ=255)
                                                                                                                                            ggTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCtt                                                   >  2:62712/1‑106 (MQ=255)
                                                                                                                                            ggTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCtt                                                   >  2:58959/1‑106 (MQ=255)
                                                                                                                                             gTTGGGCCAGGCGAAACGGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTg       <  8:13097/149‑1 (MQ=255)
                                                                                                                                             gTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTg       <  4:208706/149‑1 (MQ=255)
                                                                                                                                                  gcCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGc  >  8:93979/1‑149 (MQ=255)
                                                                                                                                                  gcCAGGCCCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGc  <  5:223959/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAA‑GGGCAAGCCCGGGTCA‑TGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGC  >  NC_002947/1711933‑1712225

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