Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,890,959 T→G 25.9% E103D (GAA→GAC)  rdgC ← Recombination‑associated protein RdgC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,890,9590TG25.9% 72.9 / 3.4 42E103D (GAA→GAC) rdgCRecombination‑associated protein RdgC
Reads supporting (aligned to +/- strand):  ref base T (5/25);  new base G (9/3);  total (14/28)
Fisher's exact test for biased strand distribution p-value = 6.28e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCGCGCTTGGCGCTGGCCGAGTTGACCAGGATCACACCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGCGCGCGGCAATCAGCAGGTACTCGCCGCTG  >  KT2440/1890825‑1891091
                                                                                                                                      |                                                                                                                                    
tCGGCGCGCTTGGCGCTGGCCGAGTTGACCAGGATCACACCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGt                                                                                                                                      <  6:132088/135‑1 (MQ=255)
  ggCGCGCTTGGCGCTGGCCGAGTTGACCAGGATCACACCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGt                                                                                                                                      <  7:235727/133‑1 (MQ=255)
        cTTGTCGCGGGCCGAGTTGACCAGGATCACACCCTGGCGCGGGGCATTGGCAGCGAAAATCATCGTGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCttttt                                                                                                                              <  5:200398/135‑1 (MQ=255)
         ttGGCGCTGGCCGAGTTGACCAGGATCACACCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTGCgtttt                                                                                                                              >  8:8906/1‑134 (MQ=255)
           ggCGCTGGCCGAGTTGACCAGGATCACCCCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTACttttt                                                                                                                              >  6:166027/1‑132 (MQ=255)
                 ggCCGAGTTGACCAGGATCACACCCAGCCGCGGGGCAATGGCACCCAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCtt                                                                                                                     <  1:320362/135‑1 (MQ=255)
                     gAGTTGACCAGGATCACACCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTGCTTTTTTTATACCTTgcgc                                                                                                                 >  2:239782/1‑135 (MQ=255)
                         tGACCCGGACCCCACCCCGGGGCGGGCAAATGGCAGCTAAAATCATCGCGCGGCGGAGGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCAttt                                                                                                             <  5:293464/135‑1 (MQ=255)
                         tGACCAGGATCACACCCAGGCGCGGGGCAATGGCATCGAAAATCATCGAGCGGCGGAGGAACGCAGGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTCCTTTTTATAGACCTTGCGCAttt                                                                                                             <  7:87743/135‑1 (MQ=255)
                         tGACCAGGATCACACCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCAttt                                                                                                             <  6:93060/135‑1 (MQ=255)
                         tGACCAGGATCACACCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCAttt                                                                                                             <  8:209244/135‑1 (MQ=255)
                                 atcaaccccAGGCGCGGGGCAATGGCATCTAAAATCATCGCGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGtc                                                                                                     <  2:157010/129‑1 (MQ=255)
                                 aTCACACCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGtc                                                                                                     <  2:307946/135‑1 (MQ=255)
                                   cacaCCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTCCGTTTTTTAGACCGTGCGCATTTGCTGGGtgtc                                                                                                   >  8:352977/1‑135 (MQ=255)
                                        ccAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTTCTTATTATAAAACTTTCGCATTTGCTCGGGATCGa                                                                                                 >  5:129625/1‑132 (MQ=255)
                                             cgcgGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTGCTTTTTTTAGACAGTGCGCCATTGCTCGGTCTCGAtctgatcg                                                                                         >  8:178951/1‑130 (MQ=255)
                                              gcgGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTCCTTTTTATATACCTTTCGCAATTGCTCGGTCTCGCTCTCttggg                                                                                        >  4:157680/1‑132 (MQ=255)
                                                  ggCAATGGCAGCGAAAATCATCGAGCGGCGGATGCACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCtt                                                                                    <  4:5582/135‑1 (MQ=255)
                                                     aaTGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTGCGTTTTTTATTACTTTCGCATTTTGTCCGTGTGGCTCTCGTAGACCattt                                                                                  >  1:4024/1‑134 (MQ=255)
                                                         gCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTGCGTTTTATAGACCTTGCGCATTTGCTCGGTATCGATCTCTTAGACCGTTTACtt                                                                              >  8:6515/1‑134 (MQ=255)
                                                          cAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTg                                                                             <  1:239782/134‑1 (MQ=255)
                                                          cAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGa                                                                            <  8:282077/135‑1 (MQ=255)
                                                              gAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACcgc                                                                        <  6:187977/135‑1 (MQ=255)
                                                                   tcatcGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTGGTTTTTATATACCTTTTGGATTTTCTAGTTCTGAATCTGTTCGACCTTTTACTTTACCgcggcgg                                                                   >  8:342019/1‑134 (MQ=255)
                                                                        gAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTAATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGAcc                                                              <  1:129434/135‑1 (MQ=255)
                                                                        gAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTCCTTTTTATAGACCTTGCGCATGTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGAcc                                                              <  7:389356/135‑1 (MQ=255)
                                                                         aGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGAcc                                                              <  4:208578/134‑1 (MQ=255)
                                                                          gcggcgGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTTCGTTTTAAATACCGTGCGCCTATGCTCGGTCTGGATCTCTTCGACCTTTTACTTGACCCCGGCCTTGACCAc                                                            >  7:58297/1‑135 (MQ=255)
                                                                                gATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTGCTTTTTTTAGTCCTGTGGGATTTTCTAGGGGTCGGTCTCTGAGTCCTTTTGCTTGCGCGCGGCGGTGACCCAgctggg                                                      >  4:95748/1‑133 (MQ=255)
                                                                                 aTGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCTGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTc                                                     <  1:246699/135‑1 (MQ=255)
                                                                                      cgcacgcGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAg                                                <  4:271485/135‑1 (MQ=255)
                                                                                        cccgcgGCAGGATCGCCTGGATGATTTCGTCCTTGATCTGGTCGCGGTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGa                                              <  1:386087/133‑1 (MQ=255)
                                                                                                       ccTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTGGATCTCTTCGACCTTTTCCTTGGCCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTgc                               <  2:179985/135‑1 (MQ=255)
                                                                                                       ccTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTgc                               <  8:31621/135‑1 (MQ=255)
                                                                                                            atgatTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCCCGCTGCTCGGCAGGATGCGTTCTTGCTTGcgcgcg                          >  8:244243/1‑135 (MQ=255)
                                                                                                            atgatTTCGTCCTTGAGCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGcgcgcg                          <  2:240741/135‑1 (MQ=255)
                                                                                                                                 tCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGCGCGCGGCAATCAGCAGGTACTcg        <  6:326856/132‑1 (MQ=255)
                                                                                                                                  cgcgTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGCGCGCGGCAATCAGCAGGTACTcgccgc    <  4:310026/135‑1 (MQ=255)
                                                                                                                                  cgcgTTCCTTTTTATAGACCTTGCGCAGTTGCTCGGTCGCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGCGCGCGGCAATCAGCAGGTACTcgccgc    <  6:245585/135‑1 (MQ=255)
                                                                                                                                   gcgTTACGTTTTATATACCTTGCGCATTTTCTAGGTCTAGATCTATTCGCCCTTTTCCTTGAACCCGTCGTTTACCACCCGGCTCGGCAGGATGCGTTCTTACTTTCGCGCGGCAATCAGCAGGTAATagccgc    >  8:51978/1‑134 (MQ=255)
                                                                                                                                    cgTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGCGCGCGGCTATCAGCAGGTACTCGCCGCTg  <  3:384053/135‑1 (MQ=255)
                                                                                                                                    cgTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGCGCGCGGCAATCAGCAGGTACTCGCCGCTg  <  7:178951/135‑1 (MQ=255)
                                                                                                                                    cgTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGCGCGCGGCAATCAGCAGGTACTCGCCGCTg  <  2:142391/135‑1 (MQ=255)
                                                                                                                                    cgTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGCGCGCGGCAATCAGCAGGTACTCGCCGCTg  <  2:146286/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
TCGGCGCGCTTGGCGCTGGCCGAGTTGACCAGGATCACACCCAGGCGCGGGGCAATGGCAGCGAAAATCATCGAGCGGCGGATGAACGCACGCGGCAGGAACGCCTGGATGATTTCGTCCTTGATCTGGTCGCGTTCCTTTTTATAGACCTTGCGCATTTGCTCGGTCTCGATCTCTTCGACCTTTTCCTTGACCGCGTCGTTGACCACGCTGCTCGGCAGGATGCGTTCTTCCTTGCGCGCGGCAATCAGCAGGTACTCGCCGCTG  >  KT2440/1890825‑1891091

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