Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,279,705 G→A A943V (GCA→GTA)  ebh_3 ← Extracellular matrix‑binding protein ebh

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,279,7050GA87.5% 35.1 / ‑4.9 16A943V (GCA→GTA) ebh_3Extracellular matrix‑binding protein ebh
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (8/6);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGTTATTAGTTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTGCTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATT  >  USA300TCH1516_ALE/2279570‑2279823
                                                                                                                                       |                                                                                                                      
gtgCCGTTATTAGTTGATTCTACTTCTGGTTTACTAGTTTCATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACtgtt                                                                                                                   <  1:205379/141‑1 (MQ=255)
      ttattaGTTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTGCTGTTGTTTCa                                                                                                             >  2:150939/1‑141 (MQ=255)
       tattaGTTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATTCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTTATGTTACTGTTGTTTCAc                                                                                                            >  2:80056/1‑141 (MQ=255)
               gATTCTACTTCTGGTTTACTAGTTACATCGTTATCTATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCGGATGTTATTGTTGTTTCACCTGTTGcc                                                                                                    <  2:221617/141‑1 (MQ=255)
               gATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGcc                                                                                                    >  1:190160/1‑141 (MQ=255)
                      cTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTCTGCCCTGTTGCCGCATccc                                                                                             >  1:133282/1‑139 (MQ=255)
                                                  cATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCACTTGCATGATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTc                                                                 >  1:34182/1‑141 (MQ=255)
                                                         ggACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCCTTTGCATCATCAGTTGCTGATGTTGCTGTTGTTTCACCTGTTGCCGCATCACTATTCTTTGGTGTTGTCGGAGAAGCGTCTGCTTTg                                                          >  1:201506/1‑141 (MQ=255)
                                                              gtttgttGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCAtt                                                     <  1:217286/141‑1 (MQ=255)
                                                                             tACACTAGAATTGTTACTAGCCTGCTGTTTTTCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGa                                      <  1:32404/141‑1 (MQ=255)
                                                                                            attaGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTcc                       <  1:70263/141‑1 (MQ=255)
                                                                                                     tggTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTc              <  1:182992/140‑1 (MQ=255)
                                                                                                           aTCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTGCGGTATTTTATGgtgt        >  1:189035/1‑141 (MQ=255)
                                                                                                           aTCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGgtgt        >  1:183726/1‑141 (MQ=255)
                                                                                                               tttGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCa    >  1:147894/1‑141 (MQ=255)
                                                                                                                 tGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATCCGTTAGGTTTTCCAGTATTTTCTGGTCTTGCAtt  >  2:228322/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                      
GTGCCGTTATTAGTTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTGCTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATT  >  USA300TCH1516_ALE/2279570‑2279823

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: