Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,460,315 T→C 13.5% intergenic (‑347/+53) 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,3150TC13.5% 67.0 / 5.9 30intergenic (‑347/+53)USA300HOU_1356/ctpApossible acetyltransferase/S41A family carboxy‑terminal peptidase
Reads supporting (aligned to +/- strand):  ref base T (13/12);  new base C (2/2);  total (16/14)
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.20e-01

TGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTAAATTTATTATTCGTTTCTTTA  >  CP000730/1460168‑1460461
                                                                                                                                                   |                                                                                                                                                  
tGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAGTTTa                                                                                                                                                   >  2:30537/1‑149 (MQ=255)
tGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTAGCTTTGCAAATTTa                                                                                                                                                   >  1:290512/1‑149 (MQ=255)
 gCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTag                                                                                                                                                  <  1:162944/149‑1 (MQ=255)
  cTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTagg                                                                                                                                                 >  1:221152/1‑149 (MQ=255)
  cTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTagg                                                                                                                                                 >  2:201034/1‑149 (MQ=255)
       tGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCg                                                                                                                                            >  1:9423/1‑149 (MQ=255)
       tGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCg                                                                                                                                            >  1:168583/1‑149 (MQ=255)
          atatTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCgtgt                                                                                                                                         >  1:147127/1‑149 (MQ=255)
                    aTGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACata                                                                                                                               <  1:235409/149‑1 (MQ=255)
                       aaTTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACataata                                                                                                                            >  1:6746/1‑149 (MQ=255)
                        aTTGTTATCAGCGATGATAGAAATTGCGTTATTTACATAATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACataataa                                                                                                                           >  1:85055/1‑149 (MQ=255)
                           gTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACataataataa                                                                                                                        >  2:48555/1‑149 (MQ=255)
                                   cGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAAttt                                                                                                                >  2:66673/1‑149 (MQ=255)
                                    gatgatAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAAtttt                                                                                                               <  1:56179/149‑1 (MQ=255)
                                            aaaTTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAgtgt                                                                                                       <  2:9423/149‑1 (MQ=255)
                                            aaaTTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCGTGTCATTAACATAATAATAATACAAATTTAGTAgtgt                                                                                                       <  2:168583/149‑1 (MQ=255)
                                                        ttACATTATTTGTAAATGTCTAATTTATATCTTTTAGAGTTATACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTAtt                                                                                           <  2:288991/149‑1 (MQ=255)
                                                             tcatTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTaaa                                                                                      >  1:27860/1‑149 (MQ=255)
                                                                       aTGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCa                                                                            <  1:201034/149‑1 (MQ=255)
                                                                                ttATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCga                                                                   <  1:32911/149‑1 (MQ=255)
                                                                                    atCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAAc                                                               >  2:216511/1‑149 (MQ=255)
                                                                                            gagTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAAcatcatcat                                                       >  2:239611/1‑149 (MQ=255)
                                                                                            gagTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAAcatcatcat                                                       >  1:187033/1‑149 (MQ=255)
                                                                                                 aGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTtta                                                  <  2:145290/149‑1 (MQ=255)
                                                                                                      tCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGc                                             <  2:211313/149‑1 (MQ=255)
                                                                                                          taaaaataaaCACGCCTCCTAAATTGCTATCTTTGCAAATTGAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCtttt                                         >  1:157093/1‑149 (MQ=255)
                                                                                                           aaaaataaaCACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCttttt                                        <  2:2252/149‑1 (MQ=255)
                                                                                                                           cctAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGt                        <  2:159953/149‑1 (MQ=255)
                                                                                                                              aaaTTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTaaa                     <  2:264622/149‑1 (MQ=255)
                                                                                                                              aaaTTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATAGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTaaa                     <  1:108746/149‑1 (MQ=255)
                                                                                                                                                 tttAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTAAATTTATTATTCGTTTCTtta  <  2:6746/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTAAATTTATTATTCGTTTCTTTA  >  CP000730/1460168‑1460461

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