Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,772,673 T→A 12.1% L87Q (CTA→CAA)  lldR → dual role activator/repressor for lldPRD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,772,6730TA12.1% 74.8 / 2.4 33L87Q (CTA→CAA) lldRdual role activator/repressor for lldPRD operon
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base A (2/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CAACTCGGCGTATCACGTAATTCACTGCGCGAGGCGCTGGCAAAACTGGTGAGTGAAGGCGTGCTGCTCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACAAAGAAAAGATTCAGCTTTGCTTTGAAGCAAC  >  CP009273/3772531‑3772817
                                                                                                                                              |                                                                                                                                                
cAACTCGGCGTATCACGTAATTCACTGCGCGAGGCGCTGGCAAAACTGGTGAGTGAGGGCGTGCTGCTCAGTCGACGCGGCGGCGGGACGATTATTCGCTGGCGTCATGACAAGTGGTCGGAGCAAAACATCGTCCAGCCGCTaa                                                                                                                                                >  2:235651/1‑145 (MQ=255)
          tagcACGTAATTCACTGCGCGAGGCGCTGGCAAAACTGGTGAGTGAAGGCGTGCTGCTCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGAt                                                                                                                                      <  2:199427/142‑1 (MQ=255)
          tATCACGTAATTCACTGCGCGAGGCGCTGGCAAAACTGGTGAGTGAAGGCGTGCTGCTCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGAt                                                                                                                                      >  2:228787/1‑145 (MQ=255)
                       aCTGCGCGAGGCGCTGGCAAAACTGGTGAGTGAAGGCGTGCTGCTCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCg                                                                                                                         <  1:227812/145‑1 (MQ=255)
                       aCTGCGCGAGGCGCTGGCAAAACTGGTGAGTGAAGGCGTGCTGCTCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCAAAAAACACTGATGGCCGATGATCCg                                                                                                                         >  2:123991/1‑145 (MQ=255)
                       aCTGCGCGAGGCGCTGGCAAAACTGGTGAGGGAAGGCGTGCTGCTCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCAAAAAACACTGATGGCCGATGATCCg                                                                                                                         >  2:332612/1‑145 (MQ=255)
                        cTGCGCGGGGCTCTGGCAAACCTGGTGAGTGCAGGCGTGCTTCTCAGTCGCCGCGGCGGCGGGACGTTTATTCGCTGGCGTCATTACACATGGTCGGAGCAAAATATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCgg                                                                                                                        <  1:79795/145‑1 (MQ=255)
                                  cgcTGGCAAAACTGGTGAGTGAAGGCGTGCTGCTCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTAATGACACCTGGTCGCACCAAACCACCGAACAGACGCTAAAATTATTGATGGCCGATGATCCGGTTTATAGttt                                                                                                              <  1:101831/145‑1 (MQ=255)
                                               gttgggtgaAGGCGTGCTGCTCCGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGaa                                                                                                 <  2:386791/140‑1 (MQ=255)
                                                         ggCGTGCTGCTCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACg                                                                                       <  2:401559/145‑1 (MQ=255)
                                                               ctgctCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATTATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTg                                                                                 <  2:376728/145‑1 (MQ=255)
                                                               ctgcctaGTGGACGCGGCGGTGGGACGTTTATTCGCTGGTGTCTTGACACATGGACGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTg                                                                                 <  2:17396/139‑1 (MQ=255)
                                                                       tCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGc                                                                         >  2:383930/1‑145 (MQ=255)
                                                                       tCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGGCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGAGCCGGATTACCGTTTCGATATCCTGGAAGGCCGCGACGCCATTGGAGCCAGc                                                                         >  2:178749/1‑145 (MQ=255)
                                                                            gcggcggcgGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCgc                                                                    >  2:104157/1‑145 (MQ=255)
                                                                                       aCGTTTATTCGCTGTGGTGATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCgg                                                         <  2:120344/145‑1 (MQ=255)
                                                                                       aCGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCgg                                                         <  2:384465/145‑1 (MQ=255)
                                                                                         gTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCa                                                       >  1:103039/1‑145 (MQ=255)
                                                                                           ttATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAAt                                                     >  2:107029/1‑145 (MQ=255)
                                                                                                          aTGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACAGCTgg                                      >  1:214575/1‑145 (MQ=255)
                                                                                                             acacGTGGTGGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGCAGCCCGCTACGCCATTGAAGCCAGCACCGCATGTCATGCGGCAATGCGCGCCACCCCTGGCGa                                   >  1:216781/1‑145 (MQ=255)
                                                                                                              cacaTGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGAc                                  <  1:228787/145‑1 (MQ=255)
                                                                                                              cacaTGGTCGGAGCAAAACATCGTCCAGCCGCAAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGAc                                  <  1:332612/145‑1 (MQ=255)
                                                                                                              cacaTGGTCGGAGCAAAACATCGCCCAGCCGCAAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGAc                                  <  1:123991/145‑1 (MQ=255)
                                                                                                                     tCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACaaagaaa                           >  1:232244/1‑145 (MQ=255)
                                                                                                                     tCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGAGTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACaaagaca                           >  1:105110/1‑143 (MQ=255)
                                                                                                                     tCGGAGCAAAACATCCTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACaaagcaa                           >  1:366393/1‑145 (MQ=255)
                                                                                                                                cATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCAGTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACAAAGAAAAGATTCAgctt                <  2:256412/145‑1 (MQ=255)
                                                                                                                                    gTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACAAAGAAAAGATTCAgctttgct            >  1:52879/1‑145 (MQ=255)
                                                                                                                                         gCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCCCATGGCATGCGGCAATGCGCGCCACACCTGGCGACAAAGAAAAGATTCAGCTTTGCTTTGaa       <  2:304642/145‑1 (MQ=255)
                                                                                                                                          ccGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACAAAGAAAAGATTCAGCTTTGCTTTGAAg      >  2:202004/1‑145 (MQ=255)
                                                                                                                                            gCTAAAACCACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACAAAGAAAAGATTCAGCTTTGCTTTGAAGCa    <  1:182792/145‑1 (MQ=255)
                                                                                                                                             cTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACAAAGAAAAGATGCAGCTTTGCTTTGAAGCaa   <  2:307276/145‑1 (MQ=255)
                                                                                                                                              tAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACAAAGAAAAGATTCAGCTTTGCTTTGAAGCAAc  >  2:197695/1‑145 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CAACTCGGCGTATCACGTAATTCACTGCGCGAGGCGCTGGCAAAACTGGTGAGTGAAGGCGTGCTGCTCAGTCGACGCGGCGGCGGGACGTTTATTCGCTGGCGTCATGACACATGGTCGGAGCAAAACATCGTCCAGCCGCTAAAAACACTGATGGCCGATGATCCGGATTACAGTTTCGATATTCTGGAAGCCCGCTACGCCATTGAAGCCAGCACCGCATGGCATGCGGCAATGCGCGCCACACCTGGCGACAAAGAAAAGATTCAGCTTTGCTTTGAAGCAAC  >  CP009273/3772531‑3772817

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: