Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,345,220 A→G 18.8% intergenic (+32/‑18) USA300HOU_1256 → / → USA300HOU_1257 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,345,2200AG18.8% 66.6 / 9.7 32intergenic (+32/‑18)USA300HOU_1256/USA300HOU_1257hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/13);  new base G (3/3);  total (16/16)
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.44e-01

ATGAATTCGAGAGCGAAACAGGAAATCTTGTTGTAGATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAAATAATTCATCTTTATTACCAAGTTATACTAT  >  CP000730/1345082‑1345359
                                                                                                                                          |                                                                                                                                           
aTGAATTCGAGAGCGAAACAGGAAATCTTGATGTAGATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGATaat                                                                                                                                           >  2:174293/1‑140 (MQ=255)
     ttCGAGAGCGAAACAGGAAATCTTGTTGTAGATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGa                                                                                                                                      >  1:95341/1‑141 (MQ=255)
         agagCGAAACAGGAAATCTTGTTGTAGATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATACTTTCTACAAAAAGCTAAAAAAGAAATg                                                                                                                                  >  1:146002/1‑141 (MQ=255)
                aaCAGGAAATCTTGTTGTAGATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTGAAAAAGAAATGCTTGaaa                                                                                                                           >  1:220266/1‑141 (MQ=255)
                aaCAGGAAATCTTGTTGTAGATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTGAAAAAGAAATGCTTGaaa                                                                                                                           >  1:25135/1‑141 (MQ=255)
                aaCAGGAAATCTTGTTGTAGATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTGAAAAAGAAATGCTTGaaa                                                                                                                           >  1:226513/1‑141 (MQ=255)
                   aGGATATCTTGTTGTAGATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGa                                                                                                                        <  2:181045/141‑1 (MQ=255)
                      aaaTCTTTTTGTAGATATACCGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAAGAATTACAAAATAGCATCCTTTCTACAAAAAGCTACAAAAGAAATGCTTGAAATGATaa                                                                                                                     >  1:180794/1‑141 (MQ=255)
                                   gATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTaaaaa                                                                                                        >  2:279511/1‑140 (MQ=255)
                                                 ttttGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGa                                                                                          >  2:114339/1‑141 (MQ=255)
                                                  tttGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGaa                                                                                         <  2:252024/141‑1 (MQ=255)
                                                   ttGCTGTTTATTCGTTTGATGAAACTGAGAAAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGaat                                                                                        <  1:225871/141‑1 (MQ=255)
                                                                      tgaAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTc                                                                     <  1:233951/141‑1 (MQ=255)
                                                                      tgaAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTc                                                                     <  1:170289/141‑1 (MQ=255)
                                                                          aCTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAACGGCagat                                                                 >  1:103006/1‑141 (MQ=255)
                                                                          aCTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCGTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCatat                                                                 >  1:151048/1‑141 (MQ=255)
                                                                                aTAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATg                                                           <  2:80372/141‑1 (MQ=255)
                                                                                 aaaaaaGTCTTGTATTATTAAAATAATTAAAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGc                                                          <  2:190999/140‑1 (MQ=255)
                                                                                        tATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATg                                                   >  2:158301/1‑141 (MQ=255)
                                                                                         aTTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGa                                                  <  1:270974/141‑1 (MQ=255)
                                                                                         aTTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGa                                                  <  1:158301/141‑1 (MQ=255)
                                                                                                 attaAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACACATGAATataata                                          >  2:209747/1‑141 (MQ=255)
                                                                                                        taatTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATAAATATAACAATACGaa                                   >  2:147853/1‑141 (MQ=255)
                                                                                                             aCAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAaata                              >  2:143841/1‑141 (MQ=255)
                                                                                                             aCAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATGAGATGCACAAATGAATATAATAATAAGAAAaata                              >  1:193762/1‑141 (MQ=255)
                                                                                                                aaaTAGCATCCTTTCTACAAAAAGCTGAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAAATAAtt                           <  2:226513/141‑1 (MQ=255)
                                                                                                                aaaTAGCATCCTTTCTACAAAAAGCTGAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAAATAAtt                           <  2:25135/141‑1 (MQ=255)
                                                                                                                aaaTAGCATCCTTTCTACAAAAAGCTGAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATAAAAAAATAAGAAAAATAAtt                           <  2:220266/141‑1 (MQ=255)
                                                                                                                 aaTAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAAATAAtta                          >  2:35301/1‑140 (MQ=255)
                                                                                                                            ttCTACAAAAAGTTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAAATAATTCATCTTTATTAc               <  1:120334/141‑1 (MQ=255)
                                                                                                                                        cTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAAATAATTCATCTTTATTACCAAGTTATACta   <  1:286816/141‑1 (MQ=255)
                                                                                                                                         tAAAAAAGAAATGCTTGAAATGATAAGTAATAGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAAATAATTCATCTTTATTACCAAGTTATACtat  <  1:30070/141‑1 (MQ=255)
                                                                                                                                         tAAAAAAGAAATGCTTGAAATGATAAGGAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAAATAATTCATCTTTATTACCAAGTTATACtat  <  2:95341/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
ATGAATTCGAGAGCGAAACAGGAAATCTTGTTGTAGATATACAGACTGATTTTGCTGTTTATTCGTTTGATGAAACTGAGATAAAAAGTATTAGATTATTAAAATAATTACAAAATAGCATCCTTTCTACAAAAAGCTAAAAAAGAAATGCTTGAAATGATAAGTAATGGTAAAATGACGATGAAATTGAATAATGTTAAGCAAAAGCGTCATATTTTATGCACAAATGAATATAATAATAAGAAAAATAATTCATCTTTATTACCAAGTTATACTAT  >  CP000730/1345082‑1345359

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: