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,9690GA97.7% 143.4 / ‑3.0 44intergenic (+509/‑107)USA300HOU_RS00140/USA300HOU_RS00350hypothetical protein/Monoacylglycerol lipase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/25);  minor base T (1/0);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-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.

CTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGT  >  USA300TCH1516_ALE/76832‑77100
                                                                                                                                         |                                                                                                                                   
ctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTATCAAGTTATGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAAGATGGGCTTAAAAGGt                                                                                                                                  <  1:601371/141‑1 (MQ=255)
                  tCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGGTAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTat                                                                                                                >  1:699847/1‑140 (MQ=255)
                                 aTTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGAttt                                                                                                 >  1:519045/1‑141 (MQ=255)
                                   tAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTgg                                                                                               <  1:178903/141‑1 (MQ=255)
                                      catcatGCCAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTc                                                                                            <  1:254637/141‑1 (MQ=255)
                                                 cAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAggg                                                                                 >  1:686318/1‑141 (MQ=255)
                                                    gTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTa                                                                              <  1:381352/141‑1 (MQ=255)
                                                    gTTAAGCGAACACTGACATGATAAATTAGTGGGTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTa                                                                              <  1:518573/141‑1 (MQ=255)
                                                     ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCCACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTCGGCCAAAAATTAGGGGTAt                                                                             <  1:456674/141‑1 (MQ=255)
                                                     ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCTAAAATTAGGGGTAt                                                                             <  1:446596/141‑1 (MQ=255)
                                                     ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTAt                                                                             <  1:653174/141‑1 (MQ=255)
                                                     ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTAt                                                                             <  1:690039/141‑1 (MQ=255)
                                                     ttAACCGAACACTGACACTATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACAATGATTTGGGTCAAAAATTAGGGGTAt                                                                             <  1:246394/141‑1 (MQ=37)
                                                          cGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTGAAAAAGTAGGGGTATTaaaa                                                                        >  1:676442/1‑141 (MQ=255)
                                                          cGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaa                                                                        >  1:544037/1‑141 (MQ=255)
                                                          cGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaa                                                                        >  1:269265/1‑141 (MQ=255)
                                                             acacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTTTTaaaaaaa                                                                     >  1:650597/1‑141 (MQ=255)
                                                             acacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTGGGGTCACAAATTAGGGGTATTaaaaaaa                                                                     >  1:89483/1‑141 (MQ=255)
                                                                               aGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAtttt                                                   >  1:63604/1‑141 (MQ=255)
                                                                                        cTATATTGTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTg                                          <  1:419142/141‑1 (MQ=255)
                                                                                          atatTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGAt                                        <  1:474689/141‑1 (MQ=255)
                                                                                          atatTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAGTTTATCAAGCTGAt                                        <  1:437418/141‑1 (MQ=255)
                                                                                                ttttACTTTGCAACAGACCCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgaga                                  <  1:116967/141‑1 (MQ=255)
                                                                                                ttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgaga                                  <  1:225779/141‑1 (MQ=255)
                                                                                                ttttACTTTGCAACAGAACCTATAAAAAACATGGGCGTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgaga                                  <  1:148999/141‑1 (MQ=255)
                                                                                                   tACTTTGCAACAGAACCTATAAAAAACATGGGCTTGAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGTGGTAGTAAAAAAATGAATGAATCATAATTTTATCAAGCTGCTGGGAAAGGt                               >  1:696074/1‑141 (MQ=255)
                                                                                                      tttGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTaa                            >  1:276546/1‑141 (MQ=255)
                                                                                                         gCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAAt                         >  1:11102/1‑141 (MQ=255)
                                                                                                          cAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATg                        >  1:464266/1‑141 (MQ=255)
                                                                                                              aGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCAtt                    >  1:62850/1‑141 (MQ=255)
                                                                                                              aGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCAtt                    >  1:170395/1‑141 (MQ=255)
                                                                                                              aGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATGTGGGGCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTAGCAAGCTGATTGGAGAGGTTAAAATGCAtt                    >  1:11716/1‑141 (MQ=255)
                                                                                                              aGAACCTATAAAAAACAAGGACTTTATTAAAAAATGTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCAtt                    >  1:195290/1‑141 (MQ=255)
                                                                                                               aaaCATATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGGCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATta                   <  1:621048/140‑1 (MQ=255)
                                                                                                                   cTATAAAAAACATGGGCTTAGAAGGTTAGTTTATATAGTATTACTGCAAACATTTATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtatata               >  1:527787/1‑141 (MQ=255)
                                                                                                                    tataCAAAACATTGGCTTAAAAGCATAATTGGTATAGTTTTCTTGTCATTACTAAGTTTGGTCAAAAATTAGGTGTATTAAAAAAATGAATAAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATaa              <  1:794785/141‑1 (MQ=255)
                                                                                                                      taAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATTATAATTTGATCAAGCTGATTGGAGAGGTTAAAATGCATTATATaaaa            <  1:157643/141‑1 (MQ=255)
                                                                                                                      taAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATaaaa            <  1:659889/141‑1 (MQ=255)
                                                                                                                            aCATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTAtt      <  1:763678/141‑1 (MQ=255)
                                                                                                                               tGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAg   >  1:165949/1‑141 (MQ=255)
                                                                                                                               tGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAg   <  1:403471/141‑1 (MQ=255)
                                                                                                                               tGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAg   <  1:428829/141‑1 (MQ=255)
                                                                                                                               ggggCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAg   <  1:327008/140‑1 (MQ=255)
                                                                                                                                gggCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGt  <  1:251290/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                   
CTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGT  >  USA300TCH1516_ALE/76832‑77100

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