Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,274,792 C→A 23.5% E182* (GAG→TAG)  radA ← DNA repair enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,274,7920CA23.5% 36.4 / 3.6 17E182* (GAG→TAG) radADNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base C (5/8);  new base A (2/2);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

AGCCAGCGAGCCCTCCTTGGTCACGTGGCCGACCAGGAAGATTGCCGTACCGCTTTGCTTGGCATAGCGCACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACCCGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTGAGTTGGTCCTGGGGTAGGCCCAGCCGCCGCGAGC  >  NC_002947/5274661‑5274926
                                                                                                                                   |                                                                                                                                      
agccagcGAGCCCTCCTTGGTCACGTGGCCGACCAGGAAGATTGCCGTACCGCTTTGCTTGGCATAGCGCACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTg                                                                                                                                    <  1:65995/136‑1 (MQ=255)
    agcGAGCCCTCCTTGGTCACGTGGCCGACCAGGAAGATTGCCGTACCGCTTTGCTTGGCATAGCGCACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTGaaaa                                                                                                                                >  1:413859/1‑136 (MQ=255)
              cctTGGTCACGTGGCCGACCAGGAAGATTGCCGTACCGCTTTGCTTGGCATAGCGCACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCTGTGAAAATGGTCTGGAt                                                                                                                      <  2:392919/136‑1 (MQ=255)
                 tGGTCACGTGGCCGACCAGGAAGATTGCCGTACCGCTTTGCTTGGCATAGCGCACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTAGGTGAAAATGGTCTGGATCGa                                                                                                                   >  1:128366/1‑136 (MQ=255)
                 tGGTCACGTGGCCGACCAGGAAGATTGCCGTACCGCTTTGCTTGGCATAGCGCACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTAGGTGAAAATGGTCTGGATCGa                                                                                                                   >  1:342981/1‑136 (MQ=255)
                    tCACGTGGCCGACCAGGAAGATTGCCGTACCGCTTTGCTTGGCATAGCGCACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTc                                                                                                                <  1:280254/136‑1 (MQ=255)
                            ccGACCAGGAAGATTGCCGTACCGCTTTGCTTGGCATAGCGCACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTCGATCACCa                                                                                                        >  1:26707/1‑136 (MQ=255)
                                                                    gcACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACCCGCGCCGTGGCAATGATCGt                                                                <  2:413859/136‑1 (MQ=255)
                                                                                     cTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACCCGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCgg                                               <  2:254697/136‑1 (MQ=255)
                                                                                                cTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTCGATCCCCATTACCCGCGGCTTCACGACCCGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTg                                    >  1:164066/1‑136 (MQ=255)
                                                                                                cTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACCCGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTg                                    <  1:366552/136‑1 (MQ=255)
                                                                                                cTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTAGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACACGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTg                                    <  2:128366/136‑1 (MQ=255)
                                                                                                cTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTAGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACACGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTg                                    <  2:342981/136‑1 (MQ=255)
                                                                                                             cccGGGGCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGCCCCGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTGAGTTGGTCCTggg                       >  2:415881/1‑136 (MQ=255)
                                                                                                                 gggCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACCCGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTGAGTTGGTCCTGGGGTAg                   >  2:45621/1‑136 (MQ=255)
                                                                                                                                gCTCGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACCCGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTGAGTTGGTCCTGGGGTAGGCCCAGCCGCCGCGa    <  2:96803/136‑1 (MQ=255)
                                                                                                                                  tCGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACCCGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTGAGTTGGTCCTGGGGTAGGCCCAGCCGCCGCGAgc  <  1:45621/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
AGCCAGCGAGCCCTCCTTGGTCACGTGGCCGACCAGGAAGATTGCCGTACCGCTTTGCTTGGCATAGCGCACCAGCAACGCTGTGCTTTCGCGCACCTGGGCCACGCCGCCCGGGGCCGATTGCAACTGCTCGGTGAAAATGGTCTGGATCGAGTCGATCACCATTACCCGCGGCTTCTCGACCCGCGCCGTGGCAATGATCGTCTCGATGCAGGTTTCGGTCATCACCTTGAGTTGGTCCTGGGGTAGGCCCAGCCGCCGCGAGC  >  NC_002947/5274661‑5274926

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