Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 315,782 T→A 12.6% Q503L (CAG→CTG)  cntO_1 ← Metal‑pseudopaline receptor CntO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440315,7820TA12.6% 74.6 / 3.9 32Q503L (CAG→CTG) cntO_1Metal‑pseudopaline receptor CntO
Reads supporting (aligned to +/- strand):  ref base T (11/17);  new base A (2/2);  total (13/19)
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.46e-01

GCCGCTGAGCAAGATGATGTTCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTA  >  KT2440/315652‑315915
                                                                                                                                  |                                                                                                                                     
gCCGCTGAGCAAGATGATGTTCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCg                                                                                                                                   >  8:150867/1‑135 (MQ=255)
     tGAGCAAGATGATGTTCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGtct                                                                                                                              >  8:228964/1‑135 (MQ=255)
          aaGATGATGTTCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTcc                                                                                                                         <  2:49346/135‑1 (MQ=255)
            gatgatGTTCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGTCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGACTATTTCTTTCTGCCGGGTCTCTTCCgg                                                                                                                       <  6:276037/135‑1 (MQ=255)
               gatgTTCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCGGTTTTCGCGGATTTATTTGTGCCCGGTGTCTTGCGGGtg                                                                                                                    >  5:330916/1‑134 (MQ=255)
                    tCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGttt                                                                                                               <  5:239221/135‑1 (MQ=255)
                    tCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTATCTGCCGGGTCTCTTCCGGGTCGGttt                                                                                                               <  2:186135/135‑1 (MQ=255)
                    tCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTACAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCGATATCTTTCAGCCGGGTCTCTTCCGGGTCGGtg                                                                                                                >  3:320398/1‑133 (MQ=255)
                    tCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAACTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCGATTTCTTTCACCCGGGTCTCTTCCGGGTccggtt                                                                                                               >  3:325747/1‑130 (MQ=255)
                    tCGGGTCGTTTGGGTCGCGGGTGCGGAGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGttt                                                                                                               <  1:125762/135‑1 (MQ=255)
                         tCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTg                                                                                                          <  5:14352/135‑1 (MQ=255)
                              tGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTAATTCTTTCTGCCGGGTCTCTTCCGGGTCGGGTTTTTTGGTGtt                                                                                                     >  8:155526/1‑135 (MQ=255)
                                 gTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTgtc                                                                                                  <  2:119726/135‑3 (MQ=255)
                                 gTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGcc                                                                                                  <  4:140461/135‑1 (MQ=255)
                                        gTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGcc                                                                                           <  4:130022/135‑1 (MQ=255)
                                                      aaGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGttt                                                                                                               >  6:72787/1‑101 (MQ=255)
                                                      aaGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGttt                                                                                                               <  5:72787/101‑1 (MQ=255)
                                                              gCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTTTTGGTGTTGCCCTGTGCCCCTGGGGTGATGGCGCTCAAgc                                                                     >  7:9836/1‑135 (MQ=255)
                                                                        ggCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGccacccac                                                           <  6:330916/135‑1 (MQ=255)
                                                                             gggTGGTGGTCAGGCGCTGATCCACCCAGTCGCGCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGGTTTTTTGGTGGTGCCCTGGGCCCCAGGGGGGATGCCGATCAAGCCGCCACCCCCCgg                                                        >  6:142111/1‑133 (MQ=255)
                                                                             gggTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTg                                                      <  8:171724/135‑1 (MQ=255)
                                                                                 ggtggtCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTATTCCGGGTCGGTTTTATTTGTGTTTGCCTTGGCACCTGGGGTGATGCCGCTAAACCCGCCACCCAACGGtttga                                                   >  6:290488/1‑130 (MQ=255)
                                                                                  gtggtCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGaaa                                                 <  2:241006/135‑1 (MQ=255)
                                                                                             cTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGGTTTATTTGTGTTGCCCTGTGCCCCTGGGGTGATGGCGATCAAGCCGCCACCCACCGGTGAGCAAGTGTTTCTCCa                                      >  5:97518/1‑135 (MQ=255)
                                                                                              tGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCC‑GGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACg                                    <  1:139736/135‑1 (MQ=255)
                                                                                                  ccAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGGTTTATTGGTGTTTCCCAGTGGCCCAGGGGGGATGCCGCTCAAGCCGCCCCCCACCCGTGAGTAAGTCTTGCGGCACGCgg                                 >  3:115306/1‑135 (MQ=255)
                                                                                                             ctctTTAAGCCTATTTCTTTCTGCCGGGTAACTTACCCGTCGGTTTTATTTGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAACCCACCACCCACCGGTGAGAAACTCTTGCACCACGATGCCTACAACGaa                      >  6:296509/1‑135 (MQ=255)
                                                                                                                   aCGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTc                 <  5:296509/134‑1 (MQ=255)
                                                                                                                       cTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGcc            <  6:76014/135‑1 (MQ=255)
                                                                                                                       cTATTTCTTTCAGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGcc            <  4:325747/135‑1 (MQ=255)
                                                                                                                       cTATTTCTTTCAGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGcc            <  4:320398/135‑1 (MQ=255)
                                                                                                                       cTAGTTCTAGCTGCCGGGGCGCTTCCGGGTCGGTGTTATGGGTGGTGCCCTGGGCCCCTGGGGTGATGCCGAGCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGcc            <  2:169953/135‑1 (MQ=255)
                                                                                                                           ttctttctGCCGGGTCTCTTCCGGGTCGGGTTTATTTGTGTTGCCCTGTGCCCCAGGGGTGATGCCGATCAAGCCGCCCCCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCCTgg        >  4:106208/1‑135 (MQ=255)
                                                                                                                             ctttctGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGCTCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGCACGAATGGTCGCGCCATGGCg      >  7:139059/1‑135 (MQ=255)
                                                                                                                             ctttctGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGCTCAAGCCGCCACCCACCGGTGAGAAAGTCTTGATCCACGAGGCGTAGAACGAATGGTCGCGCCAAGGCg      >  8:211801/1‑135 (MQ=255)
                                                                                                                                  tGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATTGTCGCGCCATGGCGTGTa  <  5:290488/134‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
GCCGCTGAGCAAGATGATGTTCGGGTCGTTTGGGTCGCGGGTGCGGCGGTTGTAAAGTTCCAGCTCGTAGATGGCCAGGGTGGTGGTCAGGCGCTGATCCAACCAGTCGCTCTTTACGCCTATTTCTTTCTGCCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTA  >  KT2440/315652‑315915

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