Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 445,241 A→G 9.7% L207P (CTT→CCT)  USA300HOU_0421 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730445,2410AG 9.7% 105.8 / 3.0 41L207P (CTT→CCT) USA300HOU_0421hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/21);  new base G (2/2);  total (18/23)
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.72e-01

GTCCTCTTGCTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTGGA  >  CP000730/445101‑445375
                                                                                                                                            |                                                                                                                                      
gTCCTCTTGCTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCaa                                                                                                                                        <  2:240723/141‑1 (MQ=255)
      ttGCTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCtt                                                                                                                                  >  1:98412/1‑141 (MQ=255)
         cTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTtca                                                                                                                               >  1:331355/1‑141 (MQ=255)
                     aCATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCGGTTAAAATTGCAGCCAGTTCTTTCATATCATAACTGt                                                                                                                   >  1:24729/1‑141 (MQ=255)
                       aTTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAg                                                                                                                 >  1:145917/1‑141 (MQ=255)
                          gATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATACCTGTAGCCt                                                                                                              >  2:270848/1‑141 (MQ=255)
                                   acggtaCCAAAGCCTTGAGGTTCATAATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAg                                                                                                     <  2:299333/136‑1 (MQ=255)
                                    atGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGt                                                                                                    <  1:255142/141‑1 (MQ=255)
                                    atGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGt                                                                                                    <  1:255195/141‑1 (MQ=255)
                                        aCCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTTCAGGAAGTTCTTTCATATCATTACTGTAGCCTGATAATAAGTAgag                                                                                                >  2:94767/1‑139 (MQ=255)
                                         ccAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGt                                                                                               <  2:58506/141‑1 (MQ=255)
                                                   tAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGTCTGATAATAAGGAGCGTTTGCCCCAAg                                                                                     <  2:62647/141‑1 (MQ=255)
                                                   tAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAAGTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAg                                                                                     <  2:62583/141‑1 (MQ=255)
                                                   tAGGTTCATCATACATTTCGTCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAg                                                                                     >  2:187430/1‑141 (MQ=255)
                                                        tcatcaTACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCt                                                                                <  1:56510/141‑1 (MQ=255)
                                                         catcatACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTg                                                                               >  1:26442/1‑141 (MQ=255)
                                                             atACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGAtt                                                                           <  2:317017/141‑1 (MQ=255)
                                                                          aaTGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAGGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCa                                                              >  1:10434/1‑141 (MQ=255)
                                                                          aaTGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAGGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCa                                                              >  1:10444/1‑141 (MQ=255)
                                                                          aaTGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCa                                                              >  2:122396/1‑141 (MQ=255)
                                                                              tctcTAATGAAACGGGCTCATATTTAATCTCTGTGCCGGATGCATAAAGGGATATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGa                                                          <  1:269761/141‑1 (MQ=255)
                                                                               ctctAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGac                                                         >  1:361096/1‑141 (MQ=255)
                                                                               ctctAATGAAACGGGCTCAGATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTACTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGac                                                         >  1:278319/1‑141 (MQ=255)
                                                                                         aCGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCAAAAGTGTCTGGATTTTTAATAATAGCAATGACACCTCTAGCa                                               >  1:68563/1‑141 (MQ=255)
                                                                                          cGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCaa                                              >  1:362148/1‑141 (MQ=255)
                                                                                          cGGGCTCATATTTAATTCCCATGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCaa                                              >  2:361436/1‑141 (MQ=255)
                                                                                             gCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAAAGACACCTCTAGCAatat                                           >  1:259817/1‑141 (MQ=255)
                                                                                               tCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGTTTTTTAATAATAGCAATGACACCTCTAGCAATATCa                                         <  2:186474/141‑1 (MQ=255)
                                                                                                  tatTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCAttt                                      <  2:336837/141‑1 (MQ=255)
                                                                                                         tttCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTaa                               <  2:230137/141‑1 (MQ=255)
                                                                                                           tCTGTGCCTGATGCCTCAGATAAAATTGCAGCAGGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAata                             <  2:10434/141‑1 (MQ=255)
                                                                                                           tCTGTGCCTGATGCCTCAGATAAAATTGCAGCAGGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAata                             <  2:10444/141‑1 (MQ=255)
                                                                                                            cTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAatat                            <  1:343647/141‑1 (MQ=255)
                                                                                                               tgCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATAtaat                         <  1:187430/141‑1 (MQ=255)
                                                                                                                    gATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCCTAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATa                    <  2:280837/141‑1 (MQ=255)
                                                                                                                      tGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACg                  <  1:140077/141‑1 (MQ=255)
                                                                                                                             gATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCg           <  2:78188/141‑1 (MQ=255)
                                                                                                                              aTAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGc          <  1:122396/141‑1 (MQ=255)
                                                                                                                                   actGCATCAAGTTCTTGCATATCATAACTGTAGCCTGATCATAGTTAGCGTGTGCCTCAAGTATCTGCATTTTTAATAATAGCAATGTCACCTCTAGCTTTATCATTTCTAGTAGTATGATTAATACGACCATTGCCAggt     <  2:141141/139‑3 (MQ=255)
                                                                                                                                    ttGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTg    <  1:287259/141‑1 (MQ=255)
                                                                                                                                     tGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTgg   <  1:360592/141‑1 (MQ=255)
                                                                                                                                      gcagcaAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGGTTGCCCCAAGTATCTGGATTTTTAATCATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTGGa  >  1:312154/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                      
GTCCTCTTGCTCCTGCGTCGTACATTGATGCCAATAATGCACCAAAGCCTTTAGGTTCATCATACATTTCTGCAAATGTCTCTAATGAAACGGGCTCATATTTAATTTCTGTGCCTGATGCCTCAGATAAAATTGCAGCAAGTTCTTTCATATCATAACTGTAGCCTGATAATAAGTAGCGTTTGCCCCAAGTATCTGGATTTTTAATAATAGCAATGACACCTCTAGCAATATCATTTCTAGTAATATAATTAATACGACCATCGCCAGCTGGA  >  CP000730/445101‑445375

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