Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,709,297 T→A 28.0% K544* (AAA→TAA)  USA300HOU_2562 ← possible acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,709,2970TA28.0% 36.6 / 16.1 25K544* (AAA→TAA) USA300HOU_2562possible acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (6/12);  new base A (2/5);  total (8/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 = 8.83e-01

GCTTTACTACCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGTTCATTTAATTGATCTTTCGTAAATGCCCCATTTGTACCAA  >  CP000730/2709170‑2709416
                                                                                                                               |                                                                                                                       
gcTTTACTACCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTAcc                                                                                                             <  2:328857/140‑1 (MQ=255)
  tttACTACCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCtt                                                                                                           <  2:71468/140‑1 (MQ=255)
      ctacCTGCATATTATAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGATTTTCAGCTGACCCACCAGTTAATGTATTTATTCTACCTTCATa                                                                                                       >  2:73819/1‑140 (MQ=255)
       tacCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATaa                                                                                                      >  2:240550/1‑140 (MQ=255)
           tGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCACAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCACAAtctc                                                                                                  <  2:196380/140‑1 (MQ=255)
                  tCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTAc                                                                                           <  1:37346/140‑1 (MQ=255)
                    tAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTActc                                                                                         <  1:200119/140‑1 (MQ=255)
                    tAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTActc                                                                                         <  1:254030/140‑1 (MQ=255)
                         gAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATCATCTCTAGGTACTCTAATa                                                                                    <  1:213580/140‑1 (MQ=255)
                                      tGCAAAGTATTCCGGATGACCTGCAGAAGCATTATACCAATCGACTGGATGTACGTTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAatata                                                                       >  2:162535/1‑140 (MQ=255)
                                          aaGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTc                                                                   <  2:321078/140‑1 (MQ=255)
                                          aaGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTc                                                                   <  2:97326/140‑1 (MQ=255)
                                          aaGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTc                                                                   <  2:240710/140‑1 (MQ=255)
                                            gTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTg                                                                 <  2:279920/140‑1 (MQ=255)
                                              aTTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCt                                                               <  1:314681/140‑1 (MQ=255)
                                               ttCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCtt                                                              <  2:22003/140‑1 (MQ=255)
                                                        aCCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAAc                                                     <  2:196108/140‑1 (MQ=255)
                                                             cAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATcc                                                >  2:52819/1‑140 (MQ=255)
                                                                taGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAAt                                             <  1:103883/139‑1 (MQ=255)
                                                                   cgttATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATTGt                                          <  2:5183/138‑1 (MQ=255)
                                                                   cTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATATATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGt                                          <  2:195523/140‑1 (MQ=255)
                                                                        taCCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACGCTAATAGATACTAAATATATGTCTGCTTTTCTAAAACTATCCACTAGTTCAtt                                     >  2:188335/1‑140 (MQ=255)
                                                                                           cATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGTTCATTTAATTGATCTTTCGTAAAt                  >  1:185439/1‑140 (MQ=255)
                                                                                                  gcgcTTTTCAGCTGCCTCATGAATTAATTAATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGTTCATTAAATTGATCTTTCGTAAATGCCCCAt           >  1:215846/1‑140 (MQ=255)
                                                                                                           aGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAATTCATTTAATTGATCTTTCGTAAATCCCCCATTTGTACCaa  >  1:23757/1‑140 (MQ=255)
                                                                                                                               |                                                                                                                       
GCTTTACTACCTGCATATTCTAAGTGAATACCGTCATATGCAAAGTATTCCGGATGACCTGCAGAAGCTTTATACCAATCGACTAGATGTACATTAGAGCGCTTTTCAGCTGCCTCATAAATTAATTTATTTATTCTACCTTCATAATCTCTAGGTACTCTAATAGAAACTAAATATATGTCTGCTTTTCCAAAACTATCCAATAGTTCATTTAATTGATCTTTCGTAAATGCCCCATTTGTACCAA  >  CP000730/2709170‑2709416

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