Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,836,062 T→A 11.0% Q114L (CAG→CTG)  hycI ← protease involved in processing C‑terminal end of HycE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,836,0620TA11.0% 86.6 / 4.1 35Q114L (CAG→CTG) hycIprotease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base T (15/16);  new base A (2/2);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

AATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCC  >  CP009273/2835923‑2836205
                                                                                                                                           |                                                                                                                                               
aaTGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCACCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGa                                                                                                                                        <  1:308002/149‑1 (MQ=255)
   gAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCAGAATGcc                                                                                                                                          >  2:332295/1‑144 (MQ=255)
   gAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCAGAATGcc                                                                                                                                          >  2:332199/1‑144 (MQ=255)
                    ccACCGCTAACGGCGCGAAGCCGCCATTTCCTGCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGc                                                                                                                         <  2:235885/144‑1 (MQ=255)
                           tAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGGGGGAAATACCTTCGCCAATCTcttctt                                                                                                             >  1:49664/1‑149 (MQ=255)
                           caccTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATcttctt                                                                                                             <  1:35614/146‑1 (MQ=255)
                                    gAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAGTCACTTCGCCAATATCTTCTTTCaa                                                                                                         <  2:235830/144‑1 (MQ=255)
                                     caGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTc                                                                                                   <  1:164968/148‑1 (MQ=255)
                                                    tccCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTaa                                                                                         <  2:321587/144‑1 (MQ=255)
                                                          ccTTCCAGTCGTTGATAAACGGGTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGGAAATCACTTCGCCAATATGTTATTTCAACTGGTCGATAGGATACATCAACGACAt                                                                              >  1:110722/1‑149 (MQ=255)
                                                            ttCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACgg                                                                                 >  2:189736/1‑144 (MQ=255)
                                                              ccAGTCGTTGATAAACGGTTTCGACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCa                                                                               <  2:319147/144‑1 (MQ=255)
                                                               cAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATAACTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTAt                                                                         >  1:28126/1‑149 (MQ=255)
                                                                    gTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCCACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAAGATCTTCTTTCAACTGGTCGATCAGGTAATTCAACGGCagtatag                                                                         >  2:136586/1‑138 (MQ=255)
                                                                        aTAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCATCTGGTCGATAAGGTAATTCAACGGCATGTTATggga                                                                     >  2:207161/1‑143 (MQ=255)
                                                                         tAAACGTTTTCTACCTCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTa                                                                    <  2:156135/144‑1 (MQ=255)
                                                                           aaCGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGtcatca                                                             <  1:237333/149‑1 (MQ=255)
                                                                            aCGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGGTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCCACTGGGCGATAAGGTACTTAAACGGCATGTTATGGTTGGTcatcat                                                            >  1:66006/1‑149 (MQ=255)
                                                                            aCGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACGGGTCGATAAGGTCATTTAACGGCATGTTATGGGTAGtc                                                                 >  2:324424/1‑144 (MQ=255)
                                                                               gTTTGTACCGCATCTTTAATCGGCTGGGTCGTCGGGTAGTAAAGGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGCTCGATAAGGTAATTCATCGGCATGTTATGGGTAGtcatc                                                              <  2:32216/144‑1 (MQ=255)
                                                                               gTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGtcatc                                                              <  2:317395/144‑1 (MQ=255)
                                                                                      ccGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACa                                                       >  2:245445/1‑144 (MQ=255)
                                                                                          ctCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTTTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGa                                              <  1:19547/148‑1 (MQ=255)
                                                                                           tCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCAGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGat                                             <  1:368659/149‑1 (MQ=255)
                                                                                           tCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCAGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGat                                             <  1:368561/149‑1 (MQ=255)
                                                                                               tAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGACTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTAGATAAGGTAATTCAACGGCCTGTTATGGGTAGTCATCATACACATCTGGGCgg                                              >  2:154817/1‑143 (MQ=255)
                                                                                                     gCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATatcat                                        >  2:35649/1‑144 (MQ=255)
                                                                                                           tCATCGGGTAGTAAAAGCCAACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAGGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAATATCTCGCCGATATCATCGGTGTCgatg                             >  1:203363/1‑149 (MQ=255)
                                                                                                               cGGGTAGTAAGAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCgat                              <  2:279073/144‑1 (MQ=255)
                                                                                                                 ggTAGTAAAAGTCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGGCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCgatga                            <  2:4139/144‑1 (MQ=255)
                                                                                                                   tagtaAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGAt                     <  1:102344/149‑1 (MQ=255)
                                                                                                                      taAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGAtctc                  >  1:180137/1‑149 (MQ=255)
                                                                                                                           gCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGGAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCgg             <  1:252792/149‑1 (MQ=255)
                                                                                                                             cgacgaTATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGc                >  2:302972/1‑144 (MQ=255)
                                                                                                                                      ccGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATATCGGCTCTATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTa       >  2:98018/1‑144 (MQ=255)
                                                                                                                                           tGATTGCCGAGGACAATCACTTCGCCAATATCTTCTCTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATCGGTAGTCTTCATAACCATCTCGGCGACATCATCCTGGTCGATGATGCGGATCTCGCCGGGGTTTAGcccc  <  2:146581/144‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
AATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCGCCAATATCTTCTTTCAACTGGTCGATAAGGTAATTCAACGGCATGTTATGGGTAGTCATCATAAACATCTCGGCGATATCATCCGGGTCGATGATGCGGATCTCGCCGGGGTTTAGCCCC  >  CP009273/2835923‑2836205

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