Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,894,035 T→G 32.8% intergenic (+31/‑467) 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
*ACB1224,894,0350TG32.8% 24.4 / 4.0 37intergenic (+31/‑467)lrpC/pucI_2HTH‑type transcriptional regulator LrpC/putative allantoin permease
Reads supporting (aligned to +/- strand):  ref base T (9/15);  new base G (6/7);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGA  >  ACB122/4893907‑4894168
                                                                                                                                 |                                                                                                                                     
tGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCcgc                                                                                                                                 <  4:307248/135‑1 (MQ=255)
tGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAACAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCcgc                                                                                                                                 <  3:218945/135‑1 (MQ=255)
gggACACCCTGCTCGACCTGTTTTTTGTCCTGCCCTACACCAAAAGGGCGATCTACACGACGCCCCGGTGCCAGGGGCGGTTGCCGCCAATGGATTGTTTGGCTCGTGTGCTGCCCTCTTCGCTGGGGGGCCcgc                                                                                                                                 <  2:400718/134‑1 (MQ=255)
  gACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCACTTCGCAGGCGGGCCcgc                                                                                                                                 >  3:99766/1‑133 (MQ=255)
     acCCTGCTCGACCTCATGAAGGTGCGGCCCGAGACAAATCCCCCGAGCATCAAGACCACCCCGGTCACGCGGCGGTTGCCGCCAATGGATTGATTGGCGCCTGTGCGGGCCTCTTCGCAGGCGGGCCCGCGAaga                                                                                                                            <  4:38358/135‑1 (MQ=255)
      cccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCCGGCGGGCCCGCGAagag                                                                                                                           <  4:99766/135‑1 (MQ=255)
      cccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCGCTTTGCCGGCGGGGCCCCGGagag                                                                                                                           >  1:19733/1‑135 (MQ=255)
         ggCTCGACCTCATGAAGGTGCGGCCCCCTCCAAAAACCGCGATCATCAAGACCCCCCCGGTAAAGCGGCGGTTGCCGCCAATGGATTGATTGGCGCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCTAAGAGGcc                                                                                                                        <  7:50767/134‑1 (MQ=255)
            tCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAg                                                                                                                      >  8:392753/1‑134 (MQ=255)
            tCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGt                                                                                                                     >  4:400907/1‑135 (MQ=255)
            tCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGt                                                                                                                     >  4:130075/1‑135 (MQ=255)
                      gAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTACCACGCGTGGCCCCGCAGCCGCCACTACACGaaaaa                                                                                                           >  4:32475/1‑134 (MQ=255)
                      caaTGGGCTGGCCGAGACAAATACCGCGATCAGCAAGAAGACCCCGGTCAAGGGGCGTTTGCCGCCATTGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAg                                                                                                           <  2:44204/134‑1 (MQ=38)
                       aaTGTGCTGGCCGATACCAATACCGCCATCAACAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCAAATGTATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaaa                                                                                                             <  4:132779/132‑1 (MQ=39)
                       aaTGTGCGGCCCGAGACCAAAACCGCCAGCAGCAAGACGCCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCGCCTGGGCTGGCCTCTTCGCAGGCGGGCCCGCGAAGAGGCCAGTACAGGAAAAg                                                                                                           <  1:265739/134‑1 (MQ=255)
                              tGGCCGAGCCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTAGTCCCCGGCCTGGCCGCGAACAGGCCCGTACCGGAAAAGCTGAGTAt                                                                                                   >  3:132601/1‑135 (MQ=255)
                                gCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCACGCGTGGCCGCGGAGAGGCCAGTACAGGAAAAGCTGAGTAttt                                                                                                 >  4:182621/1‑135 (MQ=255)
                                gCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCACATCGCCGGGGGGGCCGCGCAGAGGCCACTAAACGGAAAGCAGAGTtttt                                                                                                 >  8:22764/1‑135 (MQ=255)
                                gCCGACACCAAACCCGCCATCAACAACAACCCCCCGGTCAAGCGGGGGTTGCCGCAATTGTATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGGGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAttt                                                                                                 <  4:301206/135‑1 (MQ=37)
                                 ccGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCCCGCGTGCCCGCGGAGCGGCCAGTACAGGAAAACCTGCGTAtttt                                                                                                >  6:360183/1‑135 (MQ=255)
                                             ccGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaaaa                                                                                    <  7:392753/135‑1 (MQ=255)
                                                   tcatcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCGCTTAGCCGGCGGGCCCCCCCAGAGGCCACTACAGGAAAAGCGGAGTATTTTTTTTTAACAAAAAAGTaa                                                                              >  4:274931/1‑135 (MQ=255)
                                                   tcatcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCACATAGCCCGCGGGGCCCGGCACAGGGCAGGACAAGCAAAGCCGGGGTTTTTTTTTTCACAAAAAccata                                                                              >  3:174548/1‑131 (MQ=255)
                                                                   ggTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTAGCAGGCGTGCCCGCGGAGAGGGCCGGAACGGCAAAGCTGTGTATTTTTTTTTCAAAAAAACGTAACCATAGGTGGCCtttt                                                              >  4:96597/1‑135 (MQ=255)
                                                                      ccacgcggcgGTTGCCGCCAATGGATTGATGGGCGCCTGTGCGGGCCTCTCCGCCGGCGGGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtg                                                            <  1:409642/131‑1 (MQ=255)
                                                                        aaGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCGTGGGAGGCGGGGCCGCGACGCGGGCCGGACAGGTAAAGCTGTGTTTTTTTTTTTAACAAAAACGTAACCAAATGTGGCATTTTtttg                                                          >  7:128047/1‑135 (MQ=255)
                                                                                      gccAATGGATGGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAg                                            <  3:282023/135‑1 (MQ=255)
                                                                                            ggATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTg                                                                     <  7:155118/104‑1 (MQ=255)
                                                                                            ggATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTg                                                                     >  8:155118/1‑104 (MQ=255)
                                                                                              attgattgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcc                                     <  3:400907/134‑1 (MQ=255)
                                                                                               ttgattgGCTCCTGTGCTGGCCGCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcccc                                   <  7:365372/135‑1 (MQ=255)
                                                                                                              gcgggCCTCTTCGCCGGCGTGCCCGCGATGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACt                    <  4:199439/132‑1 (MQ=255)
                                                                                                               cgggCCTCTTCGCAGGCGGGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACtt                   <  4:16672/133‑1 (MQ=255)
                                                                                                               cTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACtt                   <  6:317187/135‑1 (MQ=255)
                                                                                                               cTGGCCTCTTCGCAGGCGGGCCCGCTAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACtt                   <  3:32475/135‑1 (MQ=255)
                                                                                                                   ccTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCa               <  3:281558/135‑1 (MQ=255)
                                                                                                                    ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAc              <  3:182621/135‑1 (MQ=255)
                                                                                                                                 tGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGa  >  8:279068/1‑134 (MQ=255)
                                                                                                                                 |                                                                                                                                     
TGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGA  >  ACB122/4893907‑4894168

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: