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

GAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCATGGCAGCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCAGCGTCAACGGCAGTTTGTCTGCCG  >  KT2440/868210‑868469
                                                                                                                                    |                                                                                                                                    
gAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGGTTGCTTGGCGGcccc                                                                                                                                >  7:117852/1‑135 (MQ=255)
   ttCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCGGCATGGCAGcccccccc                                                                                                                            >  7:342981/1‑134 (MQ=255)
        gAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGGTGCGTGGGGGcccccccccgttt                                                                                                                       >  1:300017/1‑129 (MQ=255)
         tgaaCTCCCCCGGCGGCCGCGAGCCGCCTTGCTCCGCGGCCGGCGGCGATTCCGAAGGTGATGCCCCTGTATCAGGCTGCCGCCAAGGCGTAATTGCCTGATTCCGGGGGCTGCATGGCAGCCCCCCCGGTTTTc                                                                                                                      <  8:282528/132‑1 (MQ=255)
          gTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCGGCGGGGGCGCCCCCCCGGgtttggc                                                                                                                    >  3:92297/1‑133 (MQ=255)
              tCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCATGGCAGCCCCCCGGTTTTCCCTTTa                                                                                                                <  8:310214/135‑1 (MQ=255)
                aCCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGGTGCATGGCCGCCCCCCCCGGGTTcccctgtg                                                                                                               >  8:212146/1‑130 (MQ=255)
                aCCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCATGGCCGCCCCCCCGGGTTTGCCGTTAc                                                                                                               >  7:297518/1‑135 (MQ=255)
                aCCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCGGGCTGGCCGCCCCCCCCGGTTTGCCGTTaa                                                                                                               >  7:78549/1‑134 (MQ=255)
                        tCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGGGGCCTGGGCCCCCCCCCCGGTTTTGCATTTACCCAAAcc                                                                                                        >  7:233737/1‑135 (MQ=255)
                        tCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGGGGCATGGCCGCCCCCCCCGGTTTGCCGGTTCCCaacaca                                                                                                       >  2:301402/1‑131 (MQ=255)
                                      tgctgcTCTGCCGGCGGCGTTTCGTACGGTGTTGCCGCTGTATCCGTCTCCCTCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCAGGGAACCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCa                                                                                         <  6:304429/135‑1 (MQ=255)
                                                             ccagggTGAGGCGGCTTTACCCGGCGGCCCCAACCCCTAAAAGCCTTCAATCGGGGGGGTGCAGGGCCCCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGAt                                                                   <  5:159761/130‑1 (MQ=255)
                                                                   tGAGGCCGCTGTATCCGGCGCCCCCAAACCCTAAAAGGCCTAAAGCGTGGGGCTGCAGGGAACCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGt                                                            <  5:364298/135‑1 (MQ=255)
                                                                        ccGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCCTTGGAGCCCCCCCCGGGTTGCCGGTTACCAATACCGCCCCTGGCCGTAATTTGCGAAAAGACGCGCCCGAATTCCCGTACg                                                         >  4:102759/1‑133 (MQ=255)
                                                                                                 cGTAAATGCCTGACTCCTGGGGCGGCATGCCAGCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGaa                             <  4:147175/136‑1 (MQ=255)
                                                                                                                       gggCTGCCTGCCCCCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCAGCGTCAACGGCAGtt           <  8:154474/135‑1 (MQ=255)
                                                                                                                          cTGCATGGCAGCCCCCCCCGGTTTGCCGTTTACCCAACCACCCCCAGGCCGTCATTTGCGAATAGACCGCCCCCCGTGGCACGACCAAGAACGCCCTTGGCCCCAGGCGGCGCAACACCCGCAACCGCAGTTTGt        >  2:254312/1‑135 (MQ=255)
                                                                                                                                 ggCCGCCCCCCCGGTTTTCCCTTTACCCAATCCCGCCACTGGCCTTCATTCCCGAATCGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCAGCGTCAACGGCAGTTTGTCTGCcg  <  2:77983/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCATGGCAGCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCAGCGTCAACGGCAGTTTGTCTGCCG  >  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: