Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,292,711 C→A 17.4% D132Y (GAC→TAC)  xthA ← exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,292,7110CA17.4% 64.0 / 2.4 23D132Y (GAC→TAC) xthAexodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI
Reads supporting (aligned to +/- strand):  ref base C (8/11);  new base A (2/2);  total (10/13)
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.33e-01

TCAGCCAGCGTTTGGCGTTGTCCGGGCCAATGCCGATATCACAGTCTTGCGGCGAGATGTTCATGTCGCCCATCACCAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATCGGCAAAGGTACCCCAGATGAAGCGGCGCTGG  >  NC_002947/3292585‑3292844
                                                                                                                              |                                                                                                                                     
tCAGCCAGCGTTTGGCGTTGTCCGGGCCAATGCCGATATCACAGTCTTGCGGCGAGATGTTCATGTCGCCCATCACCAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGaa                                                                                                                              >  2:456290/1‑136 (MQ=255)
tCAGCCAGCGTTTGGCGTTGTCCGGGCCAATGCCGATATCACAGTCTTGCGGCGAGATGTTCATGTCGCCCATCACCAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGaa                                                                                                                              <  1:352674/136‑1 (MQ=255)
              gCGTTGTCCGGGCCAATGCCGATATCACAGTCTTGCGGCGAGATGTTGATGTCGCCCATCACCAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCg                                                                                                                <  1:508568/136‑1 (MQ=255)
                               gCCGATATCACAGTCTTGCGGCGAGATGTTCATGTCGCCCATCACCAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTggtggggtgg                                                                                               <  1:176315/136‑1 (MQ=255)
                                               tGCGGCGAGATGTTCATGTCGCCCATCACCAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCtt                                                                               >  2:405890/1‑136 (MQ=255)
                                               tGCGGCGAGATGTTCATGTCGCCCATCACCAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCtt                                                                               >  2:116372/1‑136 (MQ=255)
                                                                 tCGCCCATCACCAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCa                                                             <  2:448833/136‑1 (MQ=255)
                                                                            cAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGGGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATc                                                  <  1:334444/136‑1 (MQ=255)
                                                                            cAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATc                                                  <  1:223497/136‑1 (MQ=255)
                                                                            cAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATc                                                  <  1:491631/136‑1 (MQ=255)
                                                                                       ggTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGGGTTGGGCGTGGGGGGAATAGCCGTTGATGATGGTGGTCGGGCTGCCggt                                       >  1:301316/1‑134 (MQ=255)
                                                                                       ggTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGGGTTAGCCGTGGGGGGAATAGCCGGTGATGATGGTGATCGGGCTGGCggt                                       >  1:291175/1‑134 (MQ=255)
                                                                                            ttGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGGGTTAGCCGTGGGGGAAATAGCCGGTCATGTTGGTTATCGGGCTGGCCTTGGCaa                                  >  1:277475/1‑135 (MQ=255)
                                                                                             tGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATc                                 <  1:116372/136‑1 (MQ=255)
                                                                                             tGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATc                                 <  1:405890/136‑1 (MQ=255)
                                                                                                 gAACTGGCTTTCCAGCAGTGCCTGCAGGTAGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATCGGCa                             >  2:373905/1‑136 (MQ=255)
                                                                                                 gAACTGGCTTTCCAGCAGTGCCTGCAGGTAGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATCGGCa                             >  2:376714/1‑136 (MQ=255)
                                                                                                               gcagTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTGGCCGTGGGGGAAATAGCCGTTGATGATGGTGATAGGGCTGGCGGTGGCATAGGCAAAGGTAACCCAGAt               >  2:292366/1‑136 (MQ=255)
                                                                                                                     ccTGCAGGTAGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATCGGCAAAGGTACCCCAGATGAAgcg         <  1:376714/136‑1 (MQ=255)
                                                                                                                     ccTGCAGGTAGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATCGGCAAAGGTACCCCAGATGAAgcg         <  1:373905/136‑1 (MQ=255)
                                                                                                                        gCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGGTGATGATGGTGATCGGGCTGCCGTTGGGATCGGCAAAGGAACCCCAGATGTAgcggcg      >  1:104462/1‑136 (MQ=255)
                                                                                                                            gTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATCGGCAAAGGTACCCCAGATGAAGCGGCGCTgg  <  1:300720/136‑1 (MQ=255)
                                                                                                                            gTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGGCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATCGGCAAAGGTACCCCAGATGAAGCGGCGCTgg  <  1:446215/136‑1 (MQ=255)
                                                                                                                              |                                                                                                                                     
TCAGCCAGCGTTTGGCGTTGTCCGGGCCAATGCCGATATCACAGTCTTGCGGCGAGATGTTCATGTCGCCCATCACCAGCACTGGCTGGTCGTTGCGGAACTGGCTTTCCAGCAGTGCCTGCAGGTCGCTGTAGAAGCGCTGCTTGGCCGGGAACTTGGTGGGGTGGTCGCGGCTTTCGCCTTGGGGGAAATAGCCGTTCATGATGGTGATCGGGCTGCCGTTGGCATCGGCAAAGGTACCCCAGATGAAGCGGCGCTGG  >  NC_002947/3292585‑3292844

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