Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,663,222 G→T 21.0% L210I (CTC→ATC)  ynfL ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,663,2220GT21.0% 45.3 / 4.5 19L210I (CTC→ATC) ynfLLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (7/8);  new base T (2/2);  total (9/10)
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.77e-01

TTACGCGCAGCCGGACTTTGTTCATGATGTTTCGGCCAGACCAACCACATTTCAGAAACCGCATCCTCTTCAGCAATCGGCACCCAGCGCATTTCGTTGAGCTGAACACGTTTAAATGACGCAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACATTCGGGTTATTTGCCAGGGGATGATCGTGCGGGATCATCGCCATAAGCGGTTCATGGACGATGACTGCGTGTTCAAGCGACTCCGGTAACGCTGTGTTACGCAGCA  >  CP009273/1662990‑1663407
                                                                                                                                                                                                                                        |                                                                                                                                                                                         
ttACGCGCAGCCGGACTTTGTTCATGATGTTTCGGCCAGACCAACCACATTTCAGAAACCGCATCCTCTTCAGCAATCGGCACCCAGCGCATTTCGTTGAGCTGAACACGTTTAAATGACGCAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCaaaaaa                                                                                                                                                 >  1:139211/1‑275 (MQ=255)
                                                                            tCGGCACCCAGCGCATTTCGTTGAGCTGAACACGTTTAAATGACGCAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGATAATATCGTCATACAGCCCTGTCCCGACGTGCGGATc                                                                                                                                                       >  1:260834/1‑193 (MQ=255)
                                                                                cACCCAGCGCATTTCGTTGAGCTGAACACGTTTAAATGACGCAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGATAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCaaaa                                                                                                                                                   <  2:260829/193‑1 (MQ=255)
                                                                                       cgcATTTCGTTGAGCTGAACACGTTTAAATGACGCAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAggg                                                                                                                                        >  2:236068/1‑197 (MQ=255)
                                                                                           tttCGTTGAGCTGAACACGTTTAAATGACGCAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCt                                                                                                                                    <  1:236073/197‑1 (MQ=255)
                                                                                               cgtaaGCAGAACAACTTTAAATGACGCAGGCAAAATTGACACACCCAGACCGGCGGAAACCAGGCCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACATTCGGGTTATTTGCCAGGGGATGATCGTGCGGGATCATCGCCATAAGCGGTTCATGGa                                                  <  2:139206/271‑1 (MQ=255)
                                                                                                                           ggCAAAATTGAAACACCCTGACCGGCTGAAACCAGTCTGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGAGTGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGAGGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACATTCGGGTTATTTGCCAGGGGATGATCGTGCGGGATCATCGCCATAAGCGGTTCATGGACGATGACTGCGTGTTCAAGCGACTCCgg                      <  1:135274/275‑1 (MQ=255)
                                                                                                                                      aaCACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAata                                                                                                                                                                                       >  1:274067/1‑103 (MQ=255)
                                                                                                                                      aaCACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCCCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAACGGGTTCTTTCGCCAGTTCAGCCAGCGTTACAATCGGGTTATTTGCCAGGGGATGATCGTGCGGGATCATCGCCCTAAGCTGTTCATGGACGATGACTGCGTGTTCAAGCGCCTCCGGTAACGCTgtgt           >  1:225112/1‑275 (MQ=255)
                                                                                                                                          cccAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGt                                                                                                                                                                                   <  2:274062/103‑1 (MQ=255)
                                                                                                                                               aCCGGCGGAAACCAGACCGTTGATGGTCATTGCCTCGCCGACCTCCTGAGTTATGACGGGCGTCAAATGGTAACGTGGCATCAGACCCAGAATATCGACATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACATTCGGGTTATTTGCCAGGGGATGATCGTGCGGGATCATCGCCATAAGCGGTTCATGGACGATGACTGCGTGTTCAAGCGACTCCGGTAACGCTGTGTTACgcagca  <  1:205011/275‑1 (MQ=255)
                                                                                                                                                         aCCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACATTCGGGTTATTTGCCAGGGGATGATCGTGCGGGATCATCGCCATAAGCgg                                                          >  1:77060/1‑209 (MQ=255)
                                                                                                                                                             gACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACATTCGGGTTATTTGCCAGGGGATGATCGTGCGGGATCATCGCCATAAGCGGTTCa                                                      <  2:77057/209‑1 (MQ=255)
                                                                                                                                                                                             tGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGATAATATCGTCATAc                                                                                                                                                                              >  1:340425/1‑57 (MQ=255)
                                                                                                                                                                                               aGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGt                                                                                                                                                                                   >  2:95391/1‑50 (MQ=255)
                                                                                                                                                                                                 tgatgaCGGGCGTCAAATGGTAACGTCGCATCAGCCCGATAATATCGTCATACAGcc                                                                                                                                                                          <  2:340419/57‑1 (MQ=255)
                                                                                                                                                                                                   atgaCGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATa                                                                                                                                                                               <  1:95395/50‑1 (MQ=255)
                                                                                                                                                                                                        gggCGTCAAATGGTAACGTCGCATCAGCCCGAGaa                                                                                                                                                                                         <  1:63771/35‑1 (MQ=255)
                                                                                                                                                                                                                                     cGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACa                                                                                                                                             >  1:111802/1‑50 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                         
TTACGCGCAGCCGGACTTTGTTCATGATGTTTCGGCCAGACCAACCACATTTCAGAAACCGCATCCTCTTCAGCAATCGGCACCCAGCGCATTTCGTTGAGCTGAACACGTTTAAATGACGCAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACATTCGGGTTATTTGCCAGGGGATGATCGTGCGGGATCATCGCCATAAGCGGTTCATGGACGATGACTGCGTGTTCAAGCGACTCCGGTAACGCTGTGTTACGCAGCA  >  CP009273/1662990‑1663407

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: