Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 868,338 A→C 56.0% intergenic (+27/+19) loiP → / ← mneA Metalloprotease LoiP/Putative manganese exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440868,3380AC56.0% ‑4.7 / 6.3 19intergenic (+27/+19)loiP/mneAMetalloprotease LoiP/Putative manganese exporter
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base C (6/3);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGCTGCAT‑GGCAG‑CCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCAGCGTCAACGGCAGTTTGTCTGCCG  >  KT2440/868210‑868469
                                                                                                                                    |                                                                                                                                    
gAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCTGGGGGTTGCTT‑GGCGG‑cccc                                                                                                                                >  7:117852/1‑135 (MQ=255)
   ttCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGCGGCAT‑GGCAG‑cccccccc                                                                                                                            >  7:342981/1‑134 (MQ=255)
        gAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGGTGCGT‑GGGGG‑cccccccccgttt                                                                                                                       >  1:300017/1‑129 (MQ=255)
         tgaaCTCCCCCGGCGGCCGCGAGCCGCCTTGCTCCGCGGCCGGCGGCGATTCCGAAGGTGATGCCCCTGTATCAGGCTGCCGCCAAGGCGT‑AATTGCC‑TGATTCCG‑GGGGCTGCAT‑GGCAG‑CCCCCCCGGTTTTc                                                                                                                      <  8:282528/132‑1 (MQ=255)
          gTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGCGGCGG‑GGGCG‑CCCCCCCGGgtttggc                                                                                                                    >  3:92297/1‑133 (MQ=255)
              tCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGCTGCAT‑GGCAG‑‑CCCCCCGGTTTTCCCTTTa                                                                                                                <  8:310214/135‑1 (MQ=255)
                aCCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGGTGCAT‑GGCCG‑CCCCCCCCGGGTTcccctgtg                                                                                                               >  8:212146/1‑130 (MQ=255)
                aCCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGCTGCAT‑GGCCG‑CCCCCCCGGGTTTGCCGTTAc                                                                                                               >  7:297518/1‑135 (MQ=255)
                aCCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGCGGGCT‑GGCCG‑CCCCCCCCGGTTTGCCGTTaa                                                                                                               >  7:78549/1‑134 (MQ=255)
                        tCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGGGGCCTGGGCCC‑CCCCCCCGGTTTTGCATTTACCCAAAcc                                                                                                        >  7:233737/1‑135 (MQ=255)
                        tCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGGGGCAT‑GGCCG‑CCCCCCCCGGTTTGCCGGTTCCCaacaca                                                                                                       >  2:301402/1‑131 (MQ=255)
                                      tgctgcTCTGCCGGCGGCGTTTCGTACGGTGTTGCCGCTGTATCCGTCTCCCTCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGCTGCAG‑GGAAC‑CCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCa                                                                                         <  6:304429/135‑1 (MQ=255)
                                                             ccagggTGAGGCGGCTTTACCCGGCGGCC‑CCAACCCCT‑AAAAGCCTTCAATCGG‑GGGGGTGCAG‑GGCCC‑CCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGAt                                                                   <  5:159761/130‑1 (MQ=255)
                                                                   tGAGGCCGCTGTATCCGGC‑GCCCCCAAACCCTAAAAGGCC‑TAAAGCGT‑GGGGCTGCAG‑GGAAC‑CCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGt                                                            <  5:364298/135‑1 (MQ=255)
                                                                        ccGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGCTGCCT‑TGGAGCCCCCCCCGGGTTGCCGGTTA‑CCAATACCGCCCCTGGCCGTAATTTGCGAAAAGACGCGCCCGAATTCCCGTACg                                                         >  4:102759/1‑133 (MQ=255)
                                                                                                 cGT‑AAATGCC‑TGACTCCT‑GGGGCGGCAT‑GCCAG‑CCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGaa                             <  4:147175/136‑1 (MQ=255)
                                                                                                                       gggCTGCCT‑GCCCC‑CCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCAGCGTCAACGGCAGtt           <  8:154474/135‑1 (MQ=255)
                                                                                                                          cTGCAT‑GGCAGCCCCCCCCGGTTTGCCGTTTACCCAA‑CCACCCCCAGGCCGTCATTTGCGAATAGACCGCCCCCCGTGGCACGACCAAGAACGCCCTTGGCCCCAGGCGGCGCAACACCCGCAACCGCAGTTTGt        >  2:254312/1‑135 (MQ=255)
                                                                                                                                 ggCCG‑CCCCCCCGGTTTTCCCTTTACCCAATCCCGCCACTGGCCTTCATTCCCGAATCGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCAGCGTCAACGGCAGTTTGTCTGCcg  <  2:77983/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGT‑AAATGCC‑TGAATCCT‑GGGGCTGCAT‑GGCAG‑CCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCAGCGTCAACGGCAGTTTGTCTGCCG  >  KT2440/868210‑868469

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