Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,893,888 A→C 20.9% intergenic (+26/‑472) lrpC → / → pucI_2 HTH‑type transcriptional regulator LrpC/putative allantoin permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,893,8880AC20.9% 36.5 / 2.1 34intergenic (+26/‑472)lrpC/pucI_2HTH‑type transcriptional regulator LrpC/putative allantoin permease
Reads supporting (aligned to +/- strand):  ref base A (7/19);  new base C (5/3);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 9.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCT  >  KT2440/4893754‑4894019
                                                                                                                                       |                                                                                                                                        
caTGGCGCAGCTGGACCCCCTCCTCGACCTCATGAATGGCCTGGCCGCCACCAATCCCGCCAGCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAg                                                                                                                                         <  3:228858/135‑1 (MQ=255)
    gAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGt                                                                                                                                     >  4:54150/1‑136 (MQ=255)
          cTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCCGGCGGGGCc‑gc                                                                                                                              >  4:25211/1‑136 (MQ=255)
           tGGACACCCTGCTCGACCTCATGACTGTGCTGGCCGAGACCAATACCGCGATCATCAAGACGCCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCc‑gc                                                                                                                              <  6:304489/135‑1 (MQ=255)
               cacCCTCCCCCACCCCATGAATGTGCTGGCCGAGCCCAATCCCGCGATCATCAAGACCACCCCGGGCAAGCGGCGGTTGCCGCCAATGGATTGAGTGGCTCCTGGGCTGGCC‑TCTTCCCCGGCGTGCCC‑GCGAag‑a                                                                                                                        <  4:104835/136‑1 (MQ=255)
                       tCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑GAGTCGCACGGGTGGCC‑CCGAAC‑AGGG‑CAG‑t                                                                                                               >  2:200564/1‑135 (MQ=255)
                                              gagaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCATCGCAGGCGTGCCC‑GCCAAC‑AGGC‑CAG‑TACACGAAAATCTGAGTAtttttt                                                                                        >  8:362998/1‑135 (MQ=255)
                                                 cccAATACCGCGATCATCAAGACGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTGCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTa                                                                                     <  6:133841/134‑1 (MQ=255)
                                                    aaTACCGCGATCATCACGACGACCCCGGCCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCTTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTacac                                                                                  <  1:207502/135‑1 (MQ=255)
                                                      tACCGCTATCAGCAAGACTCCCCCGGTCAAGCGGCGGTTGCCGCAAATGGATTTATTGGCTCCTGGGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACaaa                                                                               <  4:56816/136‑1 (MQ=255)
                                                              tcatcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGGCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGCGTATTTTTTTTACAAAAAAACG‑TAAc                                                                      >  4:44359/1‑136 (MQ=255)
                                                              tcatcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TAATCGCACGCGTGGCC‑GCGCAC‑ACGG‑CAGAGACA‑GAAAAGCAGATGATTTTTTTTTCACaaaaaactta                                                                        >  6:84547/1‑130 (MQ=255)
                                                              tcatcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TAATAGCAGGCGTGCCC‑GCGAAG‑AGGG‑CAC‑GACAGGCAAAGCTGAGTATTTTTTTTACAAAAAAAccgaa                                                                        >  7:257854/1‑131 (MQ=255)
                                                                atcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTAGCAGGCGTGGCC‑GCGGAG‑AGGG‑CAG‑GTCAGGTAAAGCTGTGTTTTTTTTTTTAACAAAAACG‑TAAc                                                                      >  5:153055/1‑134 (MQ=255)
                                                                 tcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCAt                                                                   <  4:274213/136‑1 (MQ=255)
                                                                   aagaagACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCGTGGCCGGGGTGCCC‑CCGCAG‑ACGCGCCC‑TA‑AGGAAAAACAGGGGTTTTTTTTTAAACACAAACC‑TAAcacaa                                                                  >  8:203256/1‑131 (MQ=255)
                                                                        gaCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTTGCTCCTGTGCTGGCCAAATTCCCCGGCGGCCCC‑GCCATG‑AGCC‑CAG‑AAAAGAAAAACATATTAATTTTTTTTAAAAAAAAACG‑AAACCAAAgttc                                                              >  2:293913/1‑132 (MQ=255)
                                                                         aCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGcc                                                            <  6:58697/135‑1 (MQ=255)
                                                                               gTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCGGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTtg                                                     <  4:128467/136‑1 (MQ=255)
                                                                                  aaGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑ACATGGGCGGCGTGCGC‑CCGGAGTCGGC‑CGG‑T‑CAAGAAAAACCGCGGATTTTTTTTTAACCAAAACGTTAACAAACGGGGCCTTTTTtgg                                                    >  6:302300/1‑134 (MQ=255)
                                                                                   aGCGGCTGTTGCCGCCAATGTATTTATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTtgtgt                                                  <  6:140800/135‑1 (MQ=255)
                                                                                    gcggcgGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGc                                                <  3:54150/136‑1 (MQ=255)
                                                                                                    aTGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGgcg                                <  4:22898/136‑1 (MQ=255)
                                                                                                    aTGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGgc                                 <  6:299805/135‑1 (MQ=255)
                                                                                                        attgattgGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcccc                            <  3:25211/136‑1 (MQ=255)
                                                                                                         ttgattgGCTCCTGTGCTGGCC‑TCTGGCCAGGGGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcccc                            <  7:322881/135‑1 (MQ=255)
                                                                                                          ggattgGCTCCTGTGCTGGCC‑TCTTCGCAGGCGGGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcccc                            <  2:229276/133‑1 (MQ=255)
                                                                                                                gccccTGTGGGGGCC‑GCTGCGCCGGCGGGCCC‑GCGAAG‑CGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGtt                    <  3:207117/133‑1 (MQ=255)
                                                                                                                gCTCCTCTGCCGGCC‑TCTTCGCCGGC‑TGCCCTGCTTTT‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGt                     <  5:33917/135‑1 (MQ=255)
                                                                                                                 cTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGtt                    <  6:168867/135‑1 (MQ=255)
                                                                                                                 cTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGtt                    <  7:362998/135‑1 (MQ=255)
                                                                                                                    cTGTGCTGGCC‑ACATCGGCGGCGGGGCC‑GCGGAG‑AGGC‑CAG‑GAACGGGAATGGGGTGGTTTTTTTGTACACAAAAACG‑TAAACAGAGGTGCCCTTTTTTTGTGCACGTTGTTGAGTGGAGGCGCAAATCggtta                  >  6:182492/1‑130 (MQ=255)
                                                                                                                        gCTGGCC‑GGTTCGCAGGCGGGCCC‑GCTAAT‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACt             <  5:148107/135‑1 (MQ=255)
                                                                                                                                    cGCAGGCGTGCCG‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCctct  <  3:356597/135‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
GATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCC‑TCTTCGCAGGCGTGCCC‑GCGAAG‑AGGC‑CAG‑TACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACG‑TAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCT  >  KT2440/4893754‑4894019

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