Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,143,110 G→A 40.3% V163V (GTG→GTA USA300HOU_RS10935 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,143,1100GA40.3% 2.8 / 18.8 20V163V (GTG→GTAUSA300HOU_RS10935hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (7/5);  new base A (4/4);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

CTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCTTGCT  >  USA300TCH1516_ALE/2142976‑2143242
                                                                                                                                      |                                                                                                                                    
cTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCAtt                                                                                                                                 >  2:277158/1‑140 (MQ=255)
               cGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTGAGAAAAGTAATCTTTGAAGAATTTTTTAATGTGATTAAAGATAATAGTATCGTGGTATTTATTATTGCTgcaat                                                                                                                  <  2:187436/140‑2 (MQ=255)
                  aaaaTACCTATGTTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAATGTAATCGTATCGTAGCATTTATTATTGCTACAACAGt                                                                                                               <  2:135041/140‑1 (MQ=255)
                  aaaaTACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATCATTGCTACAACAGt                                                                                                               <  1:176205/140‑1 (MQ=255)
                          tataTTTATCTTCTTTGTATCTCTTATTGGTCCTTTATTAGAAGAATATGTCTTCAGAAAAGTAATATTTGAAGAATTATTTAATGCGATTCAAGGTAATCGTATCGTGGCTTTTATCTTTGGTTAAACAGTAAGTTCAt                                                                                                       >  2:167413/1‑140 (MQ=255)
                          tataTTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCAt                                                                                                       >  2:264077/1‑140 (MQ=255)
                                   cttcttTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTg                                                                                              >  1:38917/1‑140 (MQ=255)
                                     tcttTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCa                                                                                            <  2:293195/140‑1 (MQ=255)
                                                              attaGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTa                                                                   >  2:78880/1‑140 (MQ=255)
                                                                 agaagaATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATCATTTCAAATTTATTc                                                                >  1:253516/1‑140 (MQ=255)
                                                                    agaaTATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAg                                                             <  1:277158/140‑1 (MQ=255)
                                                                    agaaTATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAg                                                             <  1:110086/140‑1 (MQ=255)
                                                                            tagtCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATGGCTACAACAGTAAGGTCATTAATAGTTGAATTAGCACCTAATGATTTCAAATTTATTCCATTTGAtttt                                                     >  2:151994/4‑140 (MQ=255)
                                                                                              ttGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCAttttt                                   >  2:245199/1‑140 (MQ=255)
                                                                                                    aaTTATTTAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTCCTATAACTGTATGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTa                             >  2:36458/1‑140 (MQ=255)
                                                                                                         tttAATGCGATTAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCata                        <  1:287411/140‑1 (MQ=255)
                                                                                                               gCGATGAAAGGTAATCGTATCGTAGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTAATTAGCATATGtttt                  >  2:179753/1‑139 (MQ=255)
                                                                                                                   ttAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACa              >  2:275926/1‑140 (MQ=255)
                                                                                                                               gtatcgtGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCTTGCt  <  1:264077/140‑1 (MQ=255)
                                                                                                                               gtatcgtGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCTTGCt  <  1:167413/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAATCTTTGGAGAATTATTTAATGCGATTAAAGGTAATCGTATCGTGGCATTTATTATTGCTACAACAGTAAGTTCATTAATATTTGCATTAGCACATAATGATTTCAAATTTATTCCAGTTTATTTTGGTATGGGTGTCATTTTTTCATTAGCATATGTTTGGACAAAACGGCTTGCT  >  USA300TCH1516_ALE/2142976‑2143242

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