Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,709,919 G→T 32.9% Y336* (TAC→TAA)  USA300HOU_2562 ← possible acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,709,9190GT32.9% 44.5 / 14.4 32Y336* (TAC→TAA) USA300HOU_2562possible acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (12/9);  new base T (2/9);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 6.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACTGTTCTAGCAAATTGCCCCTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAAATTTAGCAAATAAACTAGAA  >  CP000730/2709780‑2710057
                                                                                                                                           |                                                                                                                                          
acacTTTTCTAGCAAATTGCCCCTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAAGAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACg                                                                                                                                            <  1:100502/140‑1 (MQ=255)
    ctttCTAACAAATGGCCCCTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTtata                                                                                                                                        <  1:161410/138‑1 (MQ=255)
                ttGCCCCTACTTTTTAAGTAAAAAAGCAAAACCTTAAAATCCTTTTTTACTTTTATGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTc                                                                                                                            <  1:9925/140‑1 (MQ=255)
                     ccTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGt                                                                                                                       <  2:144304/140‑1 (MQ=255)
                            tttAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGtaat                                                                                                                >  2:301540/1‑140 (MQ=255)
                                    aaaaTGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTc                                                                                                        >  1:292953/1‑140 (MQ=255)
                                              aGCTTTAAATCCTTTTTTACGTATAGGTTTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACaa                                                                                              <  1:307178/140‑1 (MQ=255)
                                                 tttAAATCCTTTTTTACCAAACGGTTTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGa                                                                                           <  2:324589/140‑1 (MQ=255)
                                                    aaaTCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATACCTTTTCACAAAAACAATGAtcg                                                                                        >  2:190600/1‑138 (MQ=255)
                                                     aaTCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATACAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAgg                                                                                       <  1:53189/140‑1 (MQ=255)
                                                         cTTTTTTATGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATaa                                                                                   <  1:267318/140‑1 (MQ=255)
                                                             tttACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTCACAAAATGTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTATTAACTGTTAACAAAAACAATGATAGGATAATGcc                                                                               >  2:24611/1‑140 (MQ=255)
                                                               tACGTATAGTTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACTGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCAt                                                                             >  2:62756/1‑140 (MQ=255)
                                                               tACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCAt                                                                             <  1:203259/140‑1 (MQ=255)
                                                               tACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGATAACAAAAACAATGATAGGATAATGCCAt                                                                             <  2:177479/140‑1 (MQ=255)
                                                                    ataGGTGTTTGAATAAGGGGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATAGCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGta                                                                        <  2:25449/140‑1 (MQ=255)
                                                                          gtTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCt                                                                  >  2:95933/1‑140 (MQ=255)
                                                                              cAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCAATAAGTATAAGCTaaat                                                              >  2:123171/1‑140 (MQ=255)
                                                                                aTAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGCTAATGCCATAAGTATAAGCTATATGa                                                            >  1:21655/1‑140 (MQ=255)
                                                                                aTAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGa                                                            <  1:42153/140‑1 (MQ=255)
                                                                                   aaGCGATACGAAATTTCAGCCATTAACGCTTTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCg                                                         <  2:150287/140‑1 (MQ=255)
                                                                                   aaGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCg                                                         <  1:301540/140‑1 (MQ=255)
                                                                                      cGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGttt                                                      <  2:292953/140‑1 (MQ=255)
                                                                                         tACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATACCTGTTCACAAAAACAATGATAGGATAATGCCATAATTATCAGCTATATGATCGTTTAcc                                                   >  2:13011/1‑140 (MQ=255)
                                                                                          accaaaTTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCt                                                  <  1:93348/137‑1 (MQ=255)
                                                                                          aCGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCt                                                  >  1:113722/1‑140 (MQ=255)
                                                                                                   catctatcatcgtggTTAACAAACTTTCTATAATATACACTTATACCGGCATTTGTCCTTGTACGTAATAACTTTTATCATATATAAGGGTAGTATAATGCCATAAGTGTAAGCTATATGATCGTTTATCTATAATTAGt                                         <  2:75623/126‑1 (MQ=255)
                                                                                                    aGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTa                                        <  2:5150/140‑1 (MQ=255)
                                                                                                            aCGCTGTTAACAAAATTTCTCTAACATAAACTTATACCGGTATTTGTCCTTGGACCTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTc                                <  2:57168/140‑1 (MQ=255)
                                                                                                                             tCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAATTTTAGc               >  2:238268/1‑140 (MQ=255)
                                                                                                                                     ataaACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCTTTTACCTATAATTAGTAAATGTTTCATACTTAAAAATTTAGCAAATAAAc       >  1:49914/1‑140 (MQ=255)
                                                                                                                                          cGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAAATTTAGCAAATAAACTAGaa  >  2:107119/1‑140 (MQ=255)
                                                                                                                                          cGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAAATTTAGCAAATAAACTAGaa  >  1:335593/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                          
ACACTGTTCTAGCAAATTGCCCCTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAAATTTAGCAAATAAACTAGAA  >  CP000730/2709780‑2710057

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: ‑