Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,588,062 T→C 8.9% G158G (GGT→GGC)  PP_1392 → NAD‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,588,0620TC 8.9% 119.0 / 3.0 45G158G (GGT→GGC) PP_1392NAD‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base C (2/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

GACAGCCCGGATTTCCCAGAGGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGTTTACAGTCGGCTATTGA  >  NC_002947/1587922‑1588203
                                                                                                                                            |                                                                                                                                             
gACAGCCCGGATTTCCCAGAGGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCGTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATgc                                                                                                                                       >  8:4605/1‑149 (MQ=255)
    gCCCGGTTTTCCCCGAGGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTg                                                                                                                                   <  7:321614/149‑1 (MQ=255)
     cccGGATTTCCCAGAGGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTgg                                                                                                                                  <  6:352189/149‑1 (MQ=255)
        ggATTTCCCAGAGGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGc                                                                                                                                >  1:173006/1‑148 (MQ=255)
        ggATTTCCCAGAGGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGc                                                                                                                                >  4:203685/1‑148 (MQ=255)
        ggATTTCCCAGAGGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCa                                                                                                                               <  6:118208/149‑1 (MQ=255)
                cAGAGGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGatcat                                                                                                                       <  7:224363/149‑1 (MQ=255)
                 agagGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGatcatc                                                                                                                      >  8:37518/1‑149 (MQ=255)
                   agGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCCTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTg                                                                                                                    <  7:4605/149‑1 (MQ=255)
                         aaaaaaCCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTTCGGCGCGAGCAGAACCTGGACTGTACCTTCCTCTCGACTTCGGCTGAGTTCGTTGAGTGTGAGAGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATc                                                                                                              >  5:369722/3‑149 (MQ=255)
                              aCCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGc                                                                                                         >  5:354326/1‑149 (MQ=255)
                                      cAGTGCCGGCGTGCGTATTCTGGAGGCGCTGCGGCGCGAGCCGAACCTCGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGAAGAGGCTGAGCGCAGAGGTCAGTATGCGCTGGACAAGGATCATCTGATGATCGGTGCGGACGGc                                                                                                  >  6:258324/1‑148 (MQ=255)
                                               cgttcgttTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTATCGCCTTCGGCTGAGTTCGTTGAGGGTTTGCGCAGAGGTCACTATGCGCTGTGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGa                                                                                        <  5:258324/149‑1 (MQ=255)
                                                gtgcgtTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGAt                                                                                       >  8:202711/1‑149 (MQ=255)
                                                   cgtTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCAc                                                                                    <  5:80054/149‑1 (MQ=255)
                                                     ttttCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGt                                                                                  >  8:156679/1‑149 (MQ=255)
                                                     ttttCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGCCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGt                                                                                  >  3:20434/1‑149 (MQ=255)
                                                     ttttCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTAGTTGAGGGTGAGCGCAGAGGCCACTATTCGCTGGGCAAGGAGCATATGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGt                                                                                  >  3:25310/1‑149 (MQ=255)
                                                       ttCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGttt                                                                                >  1:248894/1‑149 (MQ=255)
                                                          tGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGcc                                                                             >  5:287663/1‑149 (MQ=255)
                                                            gAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTTGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCATGTTTGcc                                                                             <  1:273623/147‑1 (MQ=255)
                                                             aGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGAc                                                                          <  6:354326/149‑1 (MQ=255)
                                                              ggCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACt                                                                         <  6:287663/149‑1 (MQ=255)
                                                                           gAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGc                                                                                                  <  5:156053/111‑1 (MQ=255)
                                                                           gAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGc                                                                                                  >  6:156053/1‑111 (MQ=255)
                                                                                       gactggacCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGc                                                 <  3:203685/148‑1 (MQ=255)
                                                                                         ctggacCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTg                                               >  8:128374/1‑148 (MQ=255)
                                                                                            gacCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGgag                                            >  8:262961/1‑148 (MQ=255)
                                                                                             acCTTCCTCTCGCATTAGGCAGAGTTCGTTGAGGGTGAGCGAAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGgag                                            >  2:287365/1‑147 (MQ=255)
                                                                                               cTTCCTCTCGCCTTCTGCTGAGTTCGTTGCGGGTGATCGCAGAGGCCACTATGCGCTGGGCAAGGATCATCTTTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGc                                        <  4:25310/149‑1 (MQ=255)
                                                                                               cTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGCCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGc                                        <  4:20434/149‑1 (MQ=255)
                                                                                                    tctcGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGttt                                                                                <  4:284162/104‑1 (MQ=255)
                                                                                                    tctcGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGttt                                                                                >  3:284162/1‑104 (MQ=255)
                                                                                                     ctcGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGaaa                                   >  1:240438/1‑148 (MQ=255)
                                                                                                             cGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTcgcg                           <  2:240438/148‑1 (MQ=255)
                                                                                                              ggCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCgg                          >  8:9409/1‑148 (MQ=255)
                                                                                                                     ttCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGtt                  <  7:156679/149‑1 (MQ=255)
                                                                                                                     ttCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGtt                  >  3:25190/1‑149 (MQ=255)
                                                                                                                        gTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGtt                  <  2:6181/146‑1 (MQ=255)
                                                                                                                        gTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGtt                  >  1:6181/1‑146 (MQ=255)
                                                                                                                        gTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGTTTa                >  6:310406/1‑148 (MQ=255)
                                                                                                                         ttGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGTTTACa              <  7:128374/149‑1 (MQ=255)
                                                                                                                                   gcgcAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGTTTACAGTCGGCTAtt    <  7:202711/149‑1 (MQ=255)
                                                                                                                                   gcgcAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGAGCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGTTTACAGTCGGCTAtt    <  6:369722/149‑1 (MQ=255)
                                                                                                                                     gcAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGTTTACAGTCGGCTATTGa  <  2:173006/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
GACAGCCCGGATTTCCCAGAGGCGTATAAAACCGAGGCCAGTGCCGGCGTGCGTTTTCTGGAGGCGCTGCGGCGCGAGCCGAACCTGGACTGGACCTTCCTCTCGCCTTCGGCTGAGTTCGTTGAGGGTGAGCGCAGAGGTCACTATGCGCTGGGCAAGGATCATCTGTTGATCGGTGCGGACGGCAAGAGCTGGATCACGTTTGCCGACTATGCGATTGCCATGCTGGATGAGCTGGAGAAGCCGAAACACTCGCGGCAACGGTTTACAGTCGGCTATTGA  >  NC_002947/1587922‑1588203

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