Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,342,380 A→C 100% L78V (TTG→GTG)  USA300HOU_2199 ← possible HAD superfamily hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,342,3800AC94.3% 38.5 / ‑5.9 12L78V (TTG→GTG) USA300HOU_2199possible HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (6/5);  minor base G (1/0);  total (7/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCAAATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAAC  >  CP000730/2342259‑2342507
                                                                                                                         |                                                                                                                               
gAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAAt                                                                                                              <  2:57371/141‑1 (MQ=255)
                             ttcaaATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCAt                                                                                 <  2:113273/141‑1 (MQ=255)
                                               gCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACGCAGAGCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTg                                                               >  1:83222/1‑141 (MQ=255)
                                                 tCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGaa                                                             <  2:30788/141‑1 (MQ=255)
                                                             ttttATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTg                                                 <  1:53894/141‑1 (MQ=255)
                                                                              tCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATc                                >  2:71058/1‑141 (MQ=255)
                                                                               cAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCa                               <  1:113546/141‑1 (MQ=255)
                                                                                          ccATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCattatta                    >  1:137998/1‑141 (MQ=255)
                                                                                                 tCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTAttagtta             >  1:23707/1‑141 (MQ=255)
                                                                                                  ccAAGTGAACATTTTAGAACTCACATCAGAAAAGAGATGTAAGGGTTTAAACGTCATATTAACTTGACCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTGttagttag            >  1:119552/1‑141 (MQ=255)
                                                                                                         aaCATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTcc     >  2:54568/1‑141 (MQ=255)
                                                                                                            aTTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCcaac  >  1:131404/1‑141 (MQ=255)
                                                                                                                         |                                                                                                                               
GAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCAAATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAAC  >  CP000730/2342259‑2342507

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