Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,454,626 A→T L48L (CTA→CTT)  USA300HOU_2315 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,454,6260AT90.3% 86.6 / 3.4 31L48L (CTA→CTT) USA300HOU_2315hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (12/16);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

GAACAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTC  >  CP000730/2454488‑2454762
                                                                                                                                          |                                                                                                                                        
gaaCAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGc                                                                                                                                        >  2:198742/1‑141 (MQ=255)
 accattacaaTTTCTGGTTTTCAACGGACAATGACGATTTTCGACTTCTTTGTCATTACAATGGCGATATCGATGATGCTCAGAGATTTCCAGTCTACGAGTAGTGTCAAACACTTTTTATTTTAAGTTACAGTTCTTGct                                                                                                                                       >  2:226778/9‑140 (MQ=255)
  aCAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCAc                                                                                                                                      <  2:107648/141‑1 (MQ=255)
     aTAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACTCTTTATATTTGAAGTTACAGATCTTGCACCAt                                                                                                                                   <  1:327015/141‑1 (MQ=255)
     aTAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGACATTTCCAGTCTATAATTGGTCTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCAt                                                                                                                                   >  2:394835/1‑141 (MQ=255)
       aaaaTTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTa                                                                                                                                 >  2:267571/1‑141 (MQ=255)
             tCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATtgct                                                                                                                           <  2:278362/141‑1 (MQ=255)
              cTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATtgctg                                                                                                                          <  1:130740/141‑1 (MQ=255)
               tGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCGAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATtgctgc                                                                                                                         <  2:215555/141‑1 (MQ=255)
                ggTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCa                                                                                                                        <  2:288516/141‑1 (MQ=255)
                 gTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCaa                                                                                                                       <  1:43486/141‑1 (MQ=255)
                          gggCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTTCAATCATTTGta                                                                                                              >  2:138967/1‑141 (MQ=255)
                               aTGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTc                                                                                                         <  1:337051/141‑1 (MQ=255)
                               aTGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTc                                                                                                         <  1:19209/141‑1 (MQ=255)
                                   cGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGttt                                                                                                     >  1:325168/1‑141 (MQ=255)
                                   cGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGttt                                                                                                     >  1:302560/1‑141 (MQ=255)
                                        ttCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCaaa                                                                                                <  1:254126/141‑1 (MQ=255)
                                                  ttGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAgg                                                                                      <  2:353764/141‑1 (MQ=255)
                                                     tCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTcc                                                                                   >  2:400471/1‑141 (MQ=255)
                                                      cATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGGCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCa                                                                                  >  1:11228/1‑141 (MQ=255)
                                                              tGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAg                                                                          >  1:352785/1‑141 (MQ=255)
                                                                                  agagATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGc                                                      <  1:255931/141‑1 (MQ=255)
                                                                                          ccAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGt                                              <  1:341953/141‑1 (MQ=255)
                                                                                                     ttGGTTTCAAGCACTTTATATTTGAAGTTAAAGATTTTGCTCCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATTAAAAGGTATAACTATCTTTTTTAAAATTATAAATCCAACTGATCGTAATTCaaagac                                   >  2:301939/1‑137 (MQ=255)
                                                                                                           tCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATTAGCCTGAAAGTAATTGAACGTCTATTGc                             >  2:383207/1‑141 (MQ=255)
                                                                                                            cAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCt                            <  1:137740/141‑1 (MQ=255)
                                                                                                            cAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCt                            >  1:96309/1‑141 (MQ=255)
                                                                                                                cacTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCt                        <  2:2712/141‑1 (MQ=255)
                                                                                                                   tttATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTa                     >  1:220346/1‑141 (MQ=255)
                                                                                                                          ttGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTAc              <  1:322824/141‑1 (MQ=255)
                                                                                                                                      aTCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTc  <  2:412465/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                        
GAACAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTC  >  CP000730/2454488‑2454762

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: ‑