Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,978,089 A→C 25.5% M288R (ATG→AGG)  yqeF ← short chain acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,978,0890AC25.5% 86.0 / 8.2 49M288R (ATG→AGG) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (9/27);  new base C (10/3);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 2.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACC  >  CP009273/2977945‑2978237
                                                                                                                                                |                                                                                                                                                     
cGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAAt                                                                                                                                                  <  1:725715/150‑1 (MQ=255)
 gACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACAACCGCGCCATTTCCCTCCAtt                                                                                                                                                 >  1:739290/1‑149 (MQ=255)
         tCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCt                                                                                                                                        <  1:908283/150‑1 (MQ=255)
         tCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCt                                                                                                                                        <  2:643357/150‑1 (MQ=255)
            cACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTaca                                                                                                                                     <  2:440844/150‑1 (MQ=255)
             aCTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCCGCGCAATTCCCATCAATG‑CCGGGTCTacac                                                                                                                                    >  1:679649/1‑150 (MQ=255)
                 aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCaa                                                                                                                                 <  2:93881/149‑1 (MQ=255)
                 aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCACACAACCGCGCACTTCCCATCAAAA‑CCCGATATACACCaa                                                                                                                                 >  1:636708/1‑149 (MQ=255)
                  aGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCAATTCCCCTCAATT‑CCCGATCTACACCAAcc                                                                                                                               >  2:604098/1‑149 (MQ=255)
                   gCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGTTCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCGGCGCCATTCCCCTCAATT‑GCCGATCTACACCAAcac                                                                                                                              >  1:241526/1‑148 (MQ=255)
                          gCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCCTCCATG‑CCCGGTCTACACCAACCCCGGCaa                                                                                                                        >  2:466988/1‑149 (MQ=255)
                             aaCCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAAATCCCATCAATG‑CCGGATCTACACCAACGCCGGCCAAAgc                                                                                                                    >  2:611396/1‑150 (MQ=255)
                               ccGAAAGCGCCTGTGGAGCAAACGCTTCATTAGCCACGATAAAATCGACTTCAGCCAACTGCCAGCCTACACGCACCAGGCAACGGCGGGGCGCATACACTGGCGCAATTCCCATCCATG‑CCGGATCTACACCAACGCTGGCAAATgcgc                                                                                                                  <  2:560886/150‑1 (MQ=255)
                                 gAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCTAGGCAACGGCGGGTCGCACACACCCCCGCCAATCCCCTCAATGACCGCACCT‑CAACCACACTTGCAAATGcgcg                                                                                                                 >  2:305182/1‑149 (MQ=255)
                                 gAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAAATCCCCTCAATG‑CCCGAACTACACCAACACTGGCAAAAGcgcg                                                                                                                 >  2:885232/1‑149 (MQ=255)
                                      cgcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCTTCGGAACGATACGCAGCACTC‑CGAGATCTACACCAACGCTGGCAAATGCGCGAATGCg                                                                                                           <  2:593527/150‑1 (MQ=255)
                                      cgcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCg                                                                                                           <  1:575276/150‑1 (MQ=255)
                                       gcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCgg                                                                                                          <  1:178585/150‑1 (MQ=255)
                                        ccTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGAGCGATTCGCAGCAAAG‑CGAGATCTACACCAACGCTGGCAAATGCGCGAATGCggg                                                                                                         <  2:919694/150‑1 (MQ=255)
                                            tgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTCTTGAGCGATTCGCCGCACTC‑CGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAg                                                                                                     <  2:789572/150‑1 (MQ=255)
                                            tgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAg                                                                                                     <  1:866786/150‑1 (MQ=255)
                                                gcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACg                                                                                                 <  1:764003/150‑1 (MQ=255)
                                                gcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCGCGCCACTCCCATCAATA‑CCGGATCTACACCCACGCTGGCAAAAGCGCGAATGCGGGCCAGCACg                                                                                                 >  2:296399/1‑150 (MQ=255)
                                                 caAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCTTTTTAATGATACGGCGCCCTC‑CGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACgg                                                                                                <  2:761557/150‑1 (MQ=255)
                                                  aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTCCGGAGCAATACCCAGCAATG‑CGGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACggg                                                                                               <  2:13312/150‑1 (MQ=255)
                                                       cTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAAt                                                                                          <  1:739949/150‑1 (MQ=255)
                                                          cATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCCCGCAATTTCCCTCAAAG‑CCGGATCTACACCCAAGCTGGCAAATGCGCCAAAACGGGGCAGCACGGGGAAAttt                                                                                        >  1:992612/1‑148 (MQ=255)
                                                           aTTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCa                                                                                       <  1:434892/149‑1 (MQ=255)
                                                            ttAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAAc                                                                                     <  1:604098/150‑1 (MQ=255)
                                                            ttAGCCACGATAAGATCGACTTCAGCCAGCTGCCAGCCTGCACGCTCCAGGCAACGGCGGGTCGCATACTCCGTAATGATTCGGAGCAAAG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAAc                                                                                     <  2:807622/150‑1 (MQ=255)
                                                               gCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCGCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGc                                                                                   <  2:39466/149‑1 (MQ=255)
                                                                           tCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCCATTCCCCTCAATG‑CCGGGTCTACACCAACGCTGGCAAAAGCGCGAATGCGGGCCAGCACGGGGAAATTCAACACTCGTGCTTTGGCt                                                                      >  2:30767/1‑150 (MQ=255)
                                                                              aCTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCg                                                                   <  1:376033/150‑1 (MQ=255)
                                                                                     ccAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCCATTCCCCCCAATG‑CCGGAGCTACACCAACGCTGGCAAAAGCGCGGATGCGGGGCAGCACGGGGAAAATCCACGCCCCTGCTTTTGCTTCGCtcctca                                                            >  2:932710/1‑150 (MQ=255)
                                                                                     ccAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACATCGTAACAATACCCCGCAATG‑CCAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatca                                                            <  2:662992/150‑1 (MQ=255)
                                                                                     ccAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatca                                                            <  1:676403/150‑1 (MQ=255)
                                                                                            ccAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTa                                                      <  1:466988/149‑1 (MQ=255)
                                                                                             cAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTCCGGCGCAATACCCAGCAATG‑CCAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTAcc                                                    <  2:470379/150‑1 (MQ=255)
                                                                                                 cTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAg                                                <  1:611396/150‑1 (MQ=255)
                                                                                                    cacGCTCCAGGCAACGGCGGGTCGCATACACCGTCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTg                                             <  2:60779/150‑1 (MQ=255)
                                                                                                    cacGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCAAATCCCCTCAAAG‑CCCCGTCTACACCAACCCTTGCAAATGTGCCAATGCGGGCCAGCACGGGTCAAATCAAAGCCCGTGCTTTGGGTTCGCCCCTCACCATTACcgccgcgc                                             >  1:884890/1‑148 (MQ=255)
                                                                                                            aGGCAACGGCGGGTCGCATTTTCCGTAGCAATTCCCATCAATG‑CGGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCg                                     <  2:739290/150‑1 (MQ=255)
                                                                                                                       ggTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatga                          <  1:296399/150‑1 (MQ=255)
                                                                                                                              tACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTa                   >  1:917892/1‑150 (MQ=255)
                                                                                                                              tACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTa                   <  1:836588/150‑1 (MQ=255)
                                                                                                                                  ccGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCg               >  1:635071/1‑150 (MQ=255)
                                                                                                                                   cGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGc              >  1:210581/1‑150 (MQ=255)
                                                                                                                                    ggCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCt             >  2:763996/1‑150 (MQ=255)
                                                                                                                                    ggCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCt             >  1:326156/1‑150 (MQ=255)
                                                                                                                                                aTCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAcc  >  2:34986/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATG‑CCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACC  >  CP009273/2977945‑2978237

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑