Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 76,969 G→A 100% intergenic (+509/‑107) USA300HOU_RS00140 → / → USA300HOU_RS00350 hypothetical protein/Monoacylglycerol lipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE76,9690GA95.9% 135.9 / ‑1.2 49intergenic (+509/‑107)USA300HOU_RS00140/USA300HOU_RS00350hypothetical protein/Monoacylglycerol lipase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (26/21);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAG  >  USA300TCH1516_ALE/76833‑77099
                                                                                                                                         |                                                                                                                                  
ttCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGt                                                                                                                                 <  2:13772/140‑1 (MQ=255)
     aTCTACGGTTTTTCGCCATGCCAAGAAATTAGCATCATGCTAGAAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAACACATGGGCTTAAAAGGTTAAtt                                                                                                                            <  2:277061/140‑1 (MQ=17)
      gcTACGGATTTACGCCATGCAACGAAACTAGCATCATGCTAGTAAGTTAAGGGAACACTTAAATGATAAATTAGTTGTTAGATATCTTTTTTTACTTTGCAGCAGAACCTATAAAAAACATGTGCTTAAAAGGTTAAttt                                                                                                                           <  2:375031/139‑1 (MQ=255)
         aCGGATGTTCGGCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTtata                                                                                                                        <  2:302127/140‑1 (MQ=255)
                ttCGCCATGCCACGAAATTAGCATCACGCTAGCAAGTTAAGCGCTATTTGGTCAGATAAATTAGTGGTTAGTTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTAtt                                                                                                                 <  1:152791/140‑1 (MQ=18)
                           aCGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTATTGAAAACAt                                                                                                      >  2:458237/1‑140 (MQ=255)
                                      atcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTGTTACGGCAAACATTGATTTGGGTc                                                                                           <  1:544610/140‑1 (MQ=255)
                                      agcatGCTAGCAAGTTAAGCGAACACTTACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTc                                                                                           <  2:371882/138‑1 (MQ=255)
                                          tGCTAGAAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATTTTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATAAAGTATTACTGCAAACATTGATTTGGGTCaaaa                                                                                       <  1:64335/140‑1 (MQ=255)
                                            cTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAAt                                                                                     >  1:526241/1‑140 (MQ=255)
                                                    ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTATTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCATAAATTAGGGGTa                                                                             >  1:89936/1‑140 (MQ=255)
                                                    ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTa                                                                             >  2:484707/1‑140 (MQ=255)
                                                           aaCACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaa                                                                      >  1:389796/1‑140 (MQ=255)
                                                             cacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaaaa                                                                    <  1:452590/140‑1 (MQ=255)
                                                              actgaactgaAAAATTAGTAGTTAGTTATATTTTTTTACTTTGCAACAGAACCTATAAATAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTAAAGTGAAAAATTAGGGGTATTAAAAAaat                                                                   <  2:277942/140‑1 (MQ=37)
                                                               cTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATTAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTGAAAAATTAGGGGTATTAAAAAaatg                                                                  >  1:152924/1‑140 (MQ=255)
                                                                tGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatga                                                                 >  2:388935/1‑140 (MQ=255)
                                                                  aCATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGAAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaat                                                               >  1:432218/1‑140 (MQ=255)
                                                                    aTGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAatgaatga                                                             >  1:533474/1‑140 (MQ=255)
                                                                               gTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACTTGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGACAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAtttt                                                  >  1:8345/1‑140 (MQ=255)
                                                                               gTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAtttt                                                  >  1:396142/1‑140 (MQ=255)
                                                                                tGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATAAGAGGTATTAAAAAAATAAATGAATCATAAtgtaa                                                 >  1:73442/1‑136 (MQ=39)
                                                                                           atTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATGAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGGTgttg                                      >  2:435760/1‑137 (MQ=255)
                                                                                               ttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCGAGCTGATTGgag                                  >  2:495383/1‑140 (MQ=255)
                                                                                               ttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgag                                  >  1:231467/1‑140 (MQ=255)
                                                                                               ttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgag                                  >  2:492899/1‑140 (MQ=255)
                                                                                               ttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCAGAATGgag                                  >  2:46840/1‑140 (MQ=255)
                                                                                                tttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAAGTTAATTTATATAGTATTACTACAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTACTTGgaaa                                 >  2:76871/1‑138 (MQ=255)
                                                                                                  tattttGTAACAGAACCTATAAAAAACATGGGCATAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAgg                               <  2:128638/137‑1 (MQ=255)
                                                                                                  tACTTTGAAACAGAACCTATAAAAAACATGGGCTAAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGAGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAgg                               <  2:1862/140‑1 (MQ=255)
                                                                                                  tACGTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTAGTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAgg                               <  2:160605/140‑1 (MQ=255)
                                                                                                   attatGCAACAATACCTGTAAAAATCATGGGCATAAAATGTAAATTTATATAATATAACTCCAAACAAAGATTTGTGTTATAAATTATGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGt                              <  2:46820/136‑1 (MQ=255)
                                                                                                   aCTTTGCAACAGAACCTATAAAAGACATGGGCTTAAAAGGTTAATTTATATAGTATTGCTGCAAACATTGATTTGGGTCAAAAATTTGGGGTATTAAAAAAATGAATGAGTCATAATTTTATCAAGCTGATTGGAGAGGt                              <  2:252549/140‑1 (MQ=255)
                                                                                                    cTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGACAAAATAATAATTTTATAAAGCTGAGTGTAGAGGta                             >  1:118267/1‑139 (MQ=255)
                                                                                                    cTTTGCAAAAGAACCTATAAAAAACATGGGCTTAAAAGTTTAATTTATATAGTATTACTGCAAACTTTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGtt                             <  1:282901/140‑1 (MQ=255)
                                                                                                      ttGCAAAAGAACCTATAAAAAACATGGGCTTAAAAAGGTAAATTTATATAGTATTACTGCAAACATTAAAAAAAAAAAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTa                            <  2:432561/140‑1 (MQ=255)
                                                                                                             aGAACCTATAAAAAACATAGGCTTAAAAGGATAATTTATATAGTATTACTGCAAACATTTATTTGGGACAAAAATTAGGAGTAATAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCAt                    <  1:200078/140‑1 (MQ=255)
                                                                                                             aGAACCTACAAAAAGCATTGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGAGTTGGTTCAAAAATTAGGGCTACTAAAAAAATGAATGAATCATAAGTTTAGCAAGCTGATTGGAGAGGTTAAAATGCAt                    <  2:171474/140‑1 (MQ=255)
                                                                                                                aCCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGAATtat                 >  2:548749/1‑140 (MQ=255)
                                                                                                                 ccTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGTACAGGTTAAAATGCATtata                >  2:211501/1‑140 (MQ=255)
                                                                                                                 ccTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtgta                >  2:369875/1‑140 (MQ=255)
                                                                                                                 ccTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtata                >  2:250173/1‑140 (MQ=255)
                                                                                                                 ccTACAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAGAAGCATtgta                >  2:64240/1‑140 (MQ=255)
                                                                                                                   tataAAAAACGTGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTGGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGAGATGAGGGGGTAAAATGCATtatata              >  1:85499/1‑140 (MQ=255)
                                                                                                                   tataAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGTGGTTAAAATGAATtatata              >  2:252255/1‑140 (MQ=255)
                                                                                                                        aaaaCAGCGGCTTAAAAGGTTAATTTATATATTATTACCGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAAtt         <  1:55623/140‑1 (MQ=255)
                                                                                                                              tGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGa   >  1:515674/1‑140 (MQ=255)
                                                                                                                              tGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGa   >  1:279020/1‑140 (MQ=255)
                                                                                                                               gggCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATATGAGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAg  <  2:113286/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
TTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAG  >  USA300TCH1516_ALE/76833‑77099

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: