Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,460,199 T→A 15.2% intergenic (‑231/+169) USA300HOU_1356 ← / ← ctpA possible acetyltransferase/S41A family carboxy‑terminal peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,460,1990TA15.2% 78.0 / 8.2 33intergenic (‑231/+169)USA300HOU_1356/ctpApossible acetyltransferase/S41A family carboxy‑terminal peptidase
Reads supporting (aligned to +/- strand):  ref base T (11/17);  new base A (2/3);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

GAAAAATAGTAATTATCACACTATGGGATAAAAATAACCAATAAAAGTATCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATA  >  CP000730/1460062‑1460336
                                                                                                                                         |                                                                                                                                         
gaaAAATAGTAATTATCACACTATGGGATAAAAATAACCAATAAAAGTATCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAg                                                                                                                                        <  2:359100/141‑1 (MQ=255)
  aaaaTAGTAATTATCACACTATGGGATAAAAATAACCAATAAAAGTATCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCg                                                                                                                                      <  1:153405/141‑1 (MQ=255)
   aaaTAGTAATTATCACACTATGGGATAAAAATAACCAATAAAAGTATCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCga                                                                                                                                     >  1:241941/1‑141 (MQ=255)
        gTAATTATCACACTATGGGATAAAAATAACCAATAAAAGTATCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGCATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATa                                                                                                                                >  1:94372/1‑141 (MQ=255)
              aTCACACTATGGGATAAAAATAACCAATAAAAGTATCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAAtt                                                                                                                          >  2:179217/1‑141 (MQ=255)
                           aTAAAAATAACCAATAAAAGTATCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTAcat                                                                                                             >  2:1869/1‑141 (MQ=255)
                                    aCCAATAAAAGTATCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTaa                                                                                                    <  2:400840/141‑1 (MQ=255)
                                            aaGTATCATATACCAATTAATGGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATTTGCTCATGGCATAGTAAACCAGGAATTGGTATGAGCGATTATAGAAATTGCGGTATTTACATCATTTGGAAATGTCGaa                                                                                            >  1:456195/1‑141 (MQ=255)
                                                aTCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTta                                                                                        <  2:241941/141‑1 (MQ=255)
                                                 tCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTtat                                                                                       >  2:235763/1‑141 (MQ=255)
                                                         cAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTa                                                                               >  2:72537/1‑141 (MQ=255)
                                                         cAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTAATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTAAATCATTTGTAAATGTCTAATTTATATCttttt                                                                               >  2:413938/1‑140 (MQ=255)
                                                           attaattGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTaga                                                                             <  1:49797/141‑1 (MQ=255)
                                                             taattGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGt                                                                           <  2:150187/141‑1 (MQ=255)
                                                                     tgTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAATCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTc                                                                   <  2:180375/141‑1 (MQ=255)
                                                                         aaaaCTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAACTGTTATCAGCGATGATAGAAATTGCGTTATTTATATCTTTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCaata                                                               >  1:316883/1‑141 (MQ=255)
                                                                                      aataCTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGc                                                  <  1:87628/141‑1 (MQ=255)
                                                                                           tATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGcctcct                                             <  2:408241/141‑1 (MQ=255)
                                                                                                aaaaaGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAAtt                                        <  1:39594/141‑1 (MQ=255)
                                                                                                        tatGATCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCATCAATTGCAATCtt                                <  2:91317/141‑1 (MQ=255)
                                                                                                         atGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCttt                               <  2:180810/141‑1 (MQ=255)
                                                                                                         atGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATAttt                               >  2:414139/1‑141 (MQ=255)
                                                                                                            cTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAAtttgcaaagatagc                            <  1:5537/141‑13 (MQ=255)
                                                                                                                       ttAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTaggagg                 <  2:101318/141‑1 (MQ=255)
                                                                                                                           agcATGAATTGTTAACAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGGGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAAGCgtg             <  1:203873/139‑1 (MQ=255)
                                                                                                                                 aaTTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTa       <  1:169292/141‑1 (MQ=255)
                                                                                                                                  aTTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTaa      <  1:166767/141‑1 (MQ=255)
                                                                                                                                  aTTGTTAACAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAAt                                                                                                  >  2:365803/1‑49 (MQ=255)
                                                                                                                                  aTTGTTAACAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAAt                                                                                                  <  1:365802/49‑1 (MQ=255)
                                                                                                                                  aTTGTTAACAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAAt                                                                                                  >  2:365802/1‑49 (MQ=255)
                                                                                                                                  aTTGTTAACAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAAt                                                                                                  <  1:365803/49‑1 (MQ=255)
                                                                                                                                    tGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACa    <  2:164846/141‑1 (MQ=255)
                                                                                                                                      ttATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACata  >  1:49119/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                         
GAAAAATAGTAATTATCACACTATGGGATAAAAATAACCAATAAAAGTATCATATACCAATTAATTGATTGTTAAAACTAGTGAATAATACTATTTAAAAAGTATATGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATA  >  CP000730/1460062‑1460336

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