Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 51,334:1 (G)7→8 100% coding (36/1092 nt) USA300HOU_0043 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073051,3271.G53.8% 34.3 / 45.7 26coding (29/1092 nt)USA300HOU_0043hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base . (5/7);  new base G (7/7);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

ATTTCTTATTCATTTTAAAATTCCTCCAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTGATGTAACAGTTGAGGTTTATGAAGATTTGTTGTTATTAAGAAGTTTAAG  >  CP000730/51190‑51453
                                                                                                                                           |                                                                                                                               
aTTTCTTATTCATTTTAAAATTCCTCCAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAggg                                                                                                                              >  1:371795/1‑141 (MQ=255)
aTTTCTGATTCATTTTAAAATTCCGCCAATAATATTTAAGGAGTCAATCAATATCATATAGGTTTAACTATCTTCAATGGATTGCTATAATAAGCATAGATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAggg                                                                                                                              <  2:181645/141‑1 (MQ=255)
    cTTATTCATTTTAAAATTCCTCCAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAggggggg                                                                                                                         >  1:69770/1‑141 (MQ=255)
    cTTATTCATTTTAAAATTCCTCCAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATGAGACGCGGAACAAAAAATGGGGAggggggg                                                                                                                         >  2:251156/1‑141 (MQ=255)
           aTTTTAAACTTCATCCTATAAGTTTAAAGGAGGCAATCAATATCAAATAGGTTAAACTATCTTCAATTGATTGCTATATCAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGTatatat                                                                                                                  <  1:158445/141‑1 (MQ=255)
               tAAAATTCCTCCAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAGAATGGGGAGGGGGGGGTATATATTGa                                                                                                               <  2:241672/141‑1 (MQ=255)
                aaaaTTCCTCCAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGAt                                                                                                              >  1:338317/1‑141 (MQ=255)
                    ttCCTCCAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGTATATATTGATACCTa                                                                                                         <  1:56905/141‑1 (MQ=255)
                     tcctccAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAACATGGGGAGCGGGGGCATAGACTGATACCTAc                                                                                                        >  2:243727/1‑141 (MQ=255)
                        tccAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGTATATATTGATACCTACGCt                                                                                                     >  1:86875/1‑141 (MQ=255)
                             taatTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTagaa                                                                                                 <  2:283703/141‑1 (MQ=255)
                             aaatTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAAAAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTagaa                                                                                                 <  2:9509/140‑1 (MQ=255)
                              aatTTTTAAGGAGTCAATCAATATCATATAGGGTAAACTAACTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTagaag                                                                                                <  1:40705/141‑1 (MQ=255)
                                                       atatAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACtt                                                                       >  1:363305/1‑141 (MQ=255)
                                                        tataGGTTAACCTATCTTCAATTGATTGCTAGAATAAGTATAAGAAACACTTTATGCATAAATATACAACAATAAATGCGGAGGGGGGGTATATATTGATACCTACGCTATAAGAAATTATTGTTAAGTATGGAAACTTAg                                                                     <  2:22352/141‑1 (MQ=255)
                                                            ggTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTCTATATGGATACCTAAGCTAGAAGAAAGTATTGAAAAGAATGGAAACATAGTTg                                                                  >  2:353979/1‑141 (MQ=255)
                                                             gTTAAACTATCTTCACTTGGTTGAAATAATAAGCATAAATGAAGCATGATGATTAGACGCGGAACAAAAAATGGGGCGGGGGGGGTATATATTGATCCCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGa                                                                 <  1:8567/141‑1 (MQ=255)
                                                                     atcccaagttgtgtgCTATAATAAGCATAAGTGAAGCAGTATGTTTATACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGa                                                         <  2:355677/128‑1 (MQ=255)
                                                                     aTCTTCAATTGATTGCTATAGTAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTCAGTTGATTATTTGaa                                                        >  1:247245/1‑141 (MQ=255)
                                                                          cAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTg                                                   <  1:80212/141‑1 (MQ=255)
                                                                          cAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGTATATATTGATACCGACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTg                                                   <  1:54290/141‑1 (MQ=255)
                                                                           aaTTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTAGAAGAGATTATTGATAAGTATGGAAACTTGGTTGATTATTTGAAATTTg                                                   >  2:114585/1‑141 (MQ=255)
                                                                           aaTTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTg                                                   >  2:253406/1‑141 (MQ=255)
                                                                             ttgattgCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGGTAAGTTTGGAAACTTAGTTGGATATTTGAAAtttattt                                                >  2:23495/1‑137 (MQ=255)
                                                                                attgCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTGATGTa                                              >  1:270982/1‑141 (MQ=255)
                                                                                                       aGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTGATGTAACAGTTGAGGTTTATGAAGAttt                       >  1:383707/1‑141 (MQ=255)
                                                                                                          aTTATGACTAGAAGAAGAAAGAAAAATGGGGAGGGTGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTGATGTAACAGTTGAGGTTTATGAAGATttgtt                    <  2:212090/141‑1 (MQ=255)
                                                                                                              tGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTGATGTAACAGTTGAGGTTTATGAAGATTTGTTGtta                <  1:355663/141‑1 (MQ=255)
                                                                                                                 ttAGACGCGGAACAAAAAATGGGGAGGGGGGGTATATATTGATACCTACGATAGAAGAAATTAGTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTGATGTAACAGTTGAGGTTTATGAAGATTTGTTGTTATTaa            >  2:290586/1‑141 (MQ=255)
                                                                                                                            aCAAAAAATGGGGAGGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAACTTTGATGTAACAGTTGAGGTTTATGAAGATTTGTTGTTATTAAGAAGTTTAAg  <  1:264973/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                               
ATTTCTTATTCATTTTAAAATTCCTCCAATAATTTTTAAGGAGTCAATCAATATCATATAGGTTAAACTATCTTCAATTGATTGCTATAATAAGCATAAATGAAGCATTATGATTAGACGCGGAACAAAAAATGGGGAGGGGGGGTATATATTGATACCTACGCTAGAAGAAATTATTGATAAGTATGGAAACTTAGTTGATTATTTGAAATTTGATGTAACAGTTGAGGTTTATGAAGATTTGTTGTTATTAAGAAGTTTAAG  >  CP000730/51190‑51453

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: