Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,209,536 G→T 13.8% A216E (GCA→GAA)  rnd ← ribonuclease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,209,5360GT13.8% 72.5 / 2.7 29A216E (GCA→GAA) rndribonuclease D
Reads supporting (aligned to +/- strand):  ref base G (15/10);  new base T (2/2);  total (17/12)
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.38e-01

TCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGGCGCCGTCTTCCAGCACCCAGGCATACTTGTCGTCACTCAGCCGCGGGCGCAGCAC  >  NC_002947/5209390‑5209684
                                                                                                                                                  |                                                                                                                                                    
tctcGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATCAGATCACGCGACAAAACAAGGTGGGTCTGTTGGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCa                                                                                                                                                   <  1:471999/150‑1 (MQ=255)
    gtcgATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAgttg                                                                                                                                               >  1:529168/1‑150 (MQ=255)
        aTCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGa                                                                                                                                           <  2:75985/150‑1 (MQ=255)
             ggCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGNTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGgcgc                                                                                                                                      <  1:49333/150‑1 (MQ=255)
                  aGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTCCGCGCTTCGTGTTGCATCCATTCGCACAGCGCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCt                                                                                                                                 >  2:16591/1‑150 (MQ=255)
                      cGACAGGTTGTTCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCa                                                                                                                             >  1:368797/1‑150 (MQ=255)
                       gACAGGTTGTTCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGGTCGCGCCAGGCGCACAGCTCAAGTAGCACCGCCAGTTGTTGAGGCGCCAGCTTCCAg                                                                                                                            <  1:16591/150‑1 (MQ=255)
                                ttCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTg                                                                                                                   >  1:274752/1‑150 (MQ=255)
                                   ggctggctTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACg                                                                                                                >  1:267664/1‑150 (MQ=255)
                                     ctggctTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTc                                                                                                              <  2:368797/150‑1 (MQ=255)
                                          ttttGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGt                                                                                                         >  2:146935/1‑150 (MQ=255)
                                                   tGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTTCCAGTTGTTGAGGCGCCa                                                                                                                                    >  2:531776/1‑114 (MQ=255)
                                                   tGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTTCCAGTTGTTGAGGCGCCa                                                                                                                                    >  2:531928/1‑114 (MQ=255)
                                                   tGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTTCCAGTTGTTGAGGCGCCa                                                                                                                                    <  1:531776/114‑1 (MQ=255)
                                                   tGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTTCCAGTTGTTGAGGCGCCa                                                                                                                                    <  1:531928/114‑1 (MQ=255)
                                                       ccACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCggg                                                                                            >  2:367087/1‑150 (MQ=255)
                                                         acaGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCt                                                                                          >  1:223453/1‑150 (MQ=255)
                                                               aGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACtt                                                                                    <  2:274752/150‑1 (MQ=255)
                                                                   cTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACg                                                                                >  1:550783/1‑150 (MQ=255)
                                                                    tCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACgg                                                                               >  2:165263/1‑150 (MQ=255)
                                                                             aTACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCa                                                                      <  1:367087/150‑1 (MQ=255)
                                                                                           cacaTCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCACGCCAGCTTGACGCCGCGGTACAGGCTTTTTTGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGgc                                                        <  2:529168/150‑1 (MQ=255)
                                                                                                                                   aGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGGCGCCGTCTTCCAGCACCCAGGCATACTTGTCGTCACTCAGc                >  2:154834/1‑150 (MQ=255)
                                                                                                                                     cTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGGCGCCGTCTTCCAGCACCCAGGCATACTTGTCGTCACTCAGCcg              >  2:561392/1‑150 (MQ=255)
                                                                                                                                        aCGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGGCGCCGTCTTCCAGCACCCAGGCATACTTGTCGTCACTCAGCCGCgg           >  2:472754/1‑150 (MQ=255)
                                                                                                                                           tAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGGCGCCGTCTTCCAGCACCCAGGCATACTTGTCGTCACTCAGCCGCGGgcg        <  1:165263/150‑1 (MQ=255)
                                                                                                                                                 tGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGGCGCCGTCTTCCAGCACCCAGGCATACTTgtcgt                        <  2:86732/128‑1 (MQ=255)
                                                                                                                                                 tGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGGCGCCGTCTTCCAGCACCCAGGCATACTTgtcgt                        >  1:86732/1‑128 (MQ=255)
                                                                                                                                                 tGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGGCGCCGTCTTCCAGCACCCAGGCATACTTATCGTCACTCAGCCGCGGGCGCAGCAc  >  1:218189/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCTCGTCGATCTTGGCCAAGGCCGACAGGTTGTTCGGCTGGCTTTTGGCCATGGGCCACAGCGAGTGCTCCTTGAGGATACGATTACGCGGCACATCGCGGTTGCGCGCTTCGCGTTCGCGCCAGGCGCACAGCTCACGTAGCACTGCCAGTTGTTGAGGCGCCAGCTTCCAGGCCAGCTTGACGTCGCGGTACAGGCTTTCGGGCTCGACTTCACGGCGCAACGCAGCCACCAGCTCGGCGCCGTCTTCCAGCACCCAGGCATACTTGTCGTCACTCAGCCGCGGGCGCAGCAC  >  NC_002947/5209390‑5209684

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