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

GATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCT  >  KT2440/4893754‑4894019
                                                                                                                                       |                                                                                                                                        
caTGGCGCAGCTGGACCCCCTCCTCGACCTCATGAATGGCCTGGCCGCCACCAATCCCGCCAGCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAg                                                                                                                                         <  3:228858/135‑1 (MQ=255)
    gAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGt                                                                                                                                     >  4:54150/1‑136 (MQ=255)
          cTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCCGGCGGGGCcgc                                                                                                                              >  4:25211/1‑136 (MQ=255)
           tGGACACCCTGCTCGACCTCATGACTGTGCTGGCCGAGACCAATACCGCGATCATCAAGACGCCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCcgc                                                                                                                              <  6:304489/135‑1 (MQ=255)
               cacCCTCCCCCACCCCATGAATGTGCTGGCCGAGCCCAATCCCGCGATCATCAAGACCACCCCGGGCAAGCGGCGGTTGCCGCCAATGGATTGAGTGGCTCCTGGGCTGGCCTCTTCCCCGGCGTGCCCGCGAaga                                                                                                                        <  4:104835/136‑1 (MQ=255)
                       tCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCGAGTCGCACGGGTGGCCCCGAACAGGGCAGt                                                                                                               >  2:200564/1‑135 (MQ=255)
                                              gagaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCATCGCAGGCGTGCCCGCCAACAGGCCAGTACACGAAAATCTGAGTAtttttt                                                                                        >  8:362998/1‑135 (MQ=255)
                                                 cccAATACCGCGATCATCAAGACGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTGCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTa                                                                                     <  6:133841/134‑1 (MQ=255)
                                                    aaTACCGCGATCATCACGACGACCCCGGCCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCTTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTacac                                                                                  <  1:207502/135‑1 (MQ=255)
                                                      tACCGCTATCAGCAAGACTCCCCCGGTCAAGCGGCGGTTGCCGCAAATGGATTTATTGGCTCCTGGGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaaa                                                                               <  4:56816/136‑1 (MQ=255)
                                                              tcatcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGGCCGCGAAGAGGCCAGTACAGGAAAAGCTGCGTATTTTTTTTACAAAAAAACGTAAc                                                                      >  4:44359/1‑136 (MQ=255)
                                                              tcatcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTAATCGCACGCGTGGCCGCGCACACGGCAGAGACAGAAAAGCAGATGATTTTTTTTTCACaaaaaactta                                                                        >  6:84547/1‑130 (MQ=255)
                                                              tcatcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTAATAGCAGGCGTGCCCGCGAAGAGGGCACGACAGGCAAAGCTGAGTATTTTTTTTACAAAAAAAccgaa                                                                        >  7:257854/1‑131 (MQ=255)
                                                                atcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTAGCAGGCGTGGCCGCGGAGAGGGCAGGTCAGGTAAAGCTGTGTTTTTTTTTTTAACAAAAACGTAAc                                                                      >  5:153055/1‑134 (MQ=255)
                                                                 tcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCAt                                                                   <  4:274213/136‑1 (MQ=255)
                                                                   aagaagACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCGTGGCCGGGGTGCCCCCGCAGACGCGCCCTAAGGAAAAACAGGGGTTTTTTTTTAAACACAAACCTAAcacaa                                                                  >  8:203256/1‑131 (MQ=255)
                                                                        gaCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTTGCTCCTGTGCTGGCCAAATTCCCCGGCGGCCCCGCCATGAGCCCAGAAAAGAAAAACATATTAATTTTTTTTAAAAAAAAACGAAACCAAAgttc                                                              >  2:293913/1‑132 (MQ=255)
                                                                         aCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGcc                                                            <  6:58697/135‑1 (MQ=255)
                                                                               gTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCGGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtg                                                     <  4:128467/136‑1 (MQ=255)
                                                                                  aaGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCACATGGGCGGCGTGCGCCCGGAGTCGGCCGGTCAAGAAAAACCGCGGATTTTTTTTTAACCAAAACGTTAACAAACGGGGCCTTTTTtgg                                                    >  6:302300/1‑134 (MQ=255)
                                                                                   aGCGGCTGTTGCCGCCAATGTATTTATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgtgt                                                  <  6:140800/135‑1 (MQ=255)
                                                                                    gcggcgGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGc                                                <  3:54150/136‑1 (MQ=255)
                                                                                                    aTGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGgcg                                <  4:22898/136‑1 (MQ=255)
                                                                                                    aTGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGgc                                 <  6:299805/135‑1 (MQ=255)
                                                                                                        attgattgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcccc                            <  3:25211/136‑1 (MQ=255)
                                                                                                         ttgattgGCTCCTGTGCTGGCCTCTGGCCAGGGGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcccc                            <  7:322881/135‑1 (MQ=255)
                                                                                                          ggattgGCTCCTGTGCTGGCCTCTTCGCAGGCGGGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcccc                            <  2:229276/133‑1 (MQ=255)
                                                                                                                gccccTGTGGGGGCCGCTGCGCCGGCGGGCCCGCGAAGCGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGtt                    <  3:207117/133‑1 (MQ=255)
                                                                                                                gCTCCTCTGCCGGCCTCTTCGCCGGCTGCCCTGCTTTTAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGt                     <  5:33917/135‑1 (MQ=255)
                                                                                                                 cTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGtt                    <  6:168867/135‑1 (MQ=255)
                                                                                                                 cTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGtt                    <  7:362998/135‑1 (MQ=255)
                                                                                                                    cTGTGCTGGCCACATCGGCGGCGGGGCCGCGGAGAGGCCAGGAACGGGAATGGGGTGGTTTTTTTGTACACAAAAACGTAAACAGAGGTGCCCTTTTTTTGTGCACGTTGTTGAGTGGAGGCGCAAATCggtta                  >  6:182492/1‑130 (MQ=255)
                                                                                                                        gCTGGCCGGTTCGCAGGCGGGCCCGCTAATAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACt             <  5:148107/135‑1 (MQ=255)
                                                                                                                                    cGCAGGCGTGCCGGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCctct  <  3:356597/135‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
GATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCT  >  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: