Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,307,497 Δ1 bp 16.7% coding (751/828 nt) USA300HOU_RS06455 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,307,4970A.16.7% 49.4 / 9.5 24coding (751/828 nt)USA300HOU_RS06455hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (6/14);  new base . (2/2);  total (8/16)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

GTTATGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAA  >  USA300TCH1516_ALE/1307359‑1307632
                                                                                                                                           |                                                                                                                                       
gTTATGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTGTTTGTTGAAACAAGTACTTGTACAATGTAAGATTTTGACTATAGTAGCGAATTGATAGAAAATAAAAATCAAGAAATACATTACGAAGAAATgaa                                                                                                                                        >  2:132385/1‑140 (MQ=255)
   aTGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATgaagaa                                                                                                                                     >  2:23506/1‑140 (MQ=255)
        gTTACCATTAAACAAAGTCGTTCAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCTAATTGAGAGAAAATGAAAATCAAGAAATACATTACCAAGAAATGAAGAAacaga                                                                                                                                <  2:9500/140‑1 (MQ=255)
        gTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACCAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAACAAATACATTACGAAGGAATGAAGAAtcaga                                                                                                                                <  2:176756/140‑1 (MQ=255)
            ccATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAATTGAAGAAACAGACAGt                                                                                                                            <  1:145141/140‑1 (MQ=255)
             cATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTc                                                                                                                           >  2:85786/1‑140 (MQ=255)
                    caaAGTCGTTGAACATCCAGATTTCTTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTa                                                                                                                    <  2:67182/140‑1 (MQ=255)
                                     cAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTg                                                                                                  >  2:3555/1‑140 (MQ=255)
                                     cAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTg                                                                                                  >  2:236233/1‑140 (MQ=255)
                                        atttatttGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGga                                                                                                >  1:203144/1‑140 (MQ=255)
                                                        aaGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAACTCAAGAAATACATTACGAAGAAATGAAGAAACAGACAATCCAAAGTAAGAATAATTTAGATTTTGGAGATCAGTTtggtggtg                                                                                >  1:111810/1‑140 (MQ=255)
                                                               tGTACAATGGAAGATTTTGAGTATATTATCGAGTTGAGAGAAAATGAAAATCAAGAAATTCATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTGAGATTTTGGAGATCCGTTTGGTGGtggttttgg                                                                         <  1:105824/140‑1 (MQ=255)
                                                                                tGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAgtgt                                                        <  2:211742/140‑1 (MQ=255)
                                                                                       attatCGAATGCAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAaggagg                                                 <  2:163341/140‑1 (MQ=255)
                                                                                         tatCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCt                                               <  1:34244/140‑1 (MQ=255)
                                                                                          atCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAATAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTa                                              <  2:284409/140‑1 (MQ=255)
                                                                                          atCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGGCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTa                                              <  1:23506/140‑1 (MQ=255)
                                                                                            cGAATTGAGAGAAAATGAAAATAAAGAAATAAAATACGATATGATTAAAGAAATGACAGAAGGAAATAATTATAACATAAATGTGGGATATCCGTTGGGTGATGGTTTTTGTTTCTTCGAGTAAGTGTAAGGAGGCTATg                                            <  2:68560/140‑1 (MQ=255)
                                                                                                      gaaaatgaaaatCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGAAAGAATAATTTAGATTTTGGAGATCCGGTTGGTGGAGGTGTTGGTTTCTTCGAGCAAGTGTAAGGAGGCAATGAATAGAacac                                  >  2:275712/1‑140 (MQ=255)
                                                                                                                  cAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTCGGTGGTGGTTTTGGTTTCATCGAGTATGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAGgc                      >  2:158912/1‑140 (MQ=255)
                                                                                                                        aTACATTACAAAGAAAAGAAGAAACAGACAGTCCAAAGTAAGAATAATGTAGATATTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAGGTGTAAGAAGGCTATGAATTGAAAACGCTCGGTGAAGCGCTaaa                >  2:61812/1‑140 (MQ=255)
                                                                                                                           cagtACGAAGAAATGAAGAAACAGACTGTCCAAAGTAATAATGATTTAGATGTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTAGAGTAAGTGTAAGGAGGCTATGTATTGAAAACGGTCGGTGAAGCGCTAAAAgg             <  2:177005/137‑1 (MQ=255)
                                                                                                                            aTTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGGGTAATTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGt            <  1:184029/140‑1 (MQ=255)
                                                                                                                            aTTACGAAGAAATGATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTa           <  1:3555/140‑1 (MQ=255)
                                                                                                                            aTTACGAAGAAATGATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTa           <  1:236233/140‑1 (MQ=255)
                                                                                                                                       aaTGAAGAAACCGCCAGTCCAAAGTAATAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAGCGGTCGGTGAAGCGCTAAAAGGTAGACGTGaaa  <  1:158912/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
GTTATGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAA  >  USA300TCH1516_ALE/1307359‑1307632

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: