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% 137.9 / ‑4.3 43intergenic (+509/‑107)USA300HOU_RS00140/USA300HOU_RS00350hypothetical protein/Monoacylglycerol lipase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/18);  minor base C (1/0);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-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.

GATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCA  >  USA300TCH1516_ALE/76837‑77102
                                                                                                                                    |                                                                                                                                      
gATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAAtt                                                                                                                                <  1:670169/141‑1 (MQ=255)
      cGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCACCTGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTtatata                                                                                                                          >  1:373763/1‑141 (MQ=255)
          ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTAt                                                                                                                      <  1:870081/141‑1 (MQ=255)
          ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTAt                                                                                                                      <  1:276715/141‑1 (MQ=255)
                      cACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACAt                                                                                                          <  1:152798/141‑1 (MQ=255)
                      cACGAAACTAGCATCATGCTANCNAGTNAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTCATTTATATAGGATTACTGCAAACAt                                                                                                          <  1:477467/141‑1 (MQ=18)
                          aaaTTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTCGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGAt                                                                                                      >  1:848558/1‑141 (MQ=255)
                             ttAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTg                                                                                                   <  1:5696/141‑1 (MQ=255)
                              tAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTgg                                                                                                  >  1:507483/1‑141 (MQ=255)
                                          aGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAg                                                                                      >  1:730020/1‑141 (MQ=255)
                                                    gCGAACACTGCCATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTTCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaa                                                                           <  1:382042/141‑1 (MQ=255)
                                                        acacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaaaa                                                                       >  1:800284/1‑141 (MQ=255)
                                                             gACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGATTAAAAATTCGGGGTATTAAAAAaatgaat                                                                  >  1:328709/1‑141 (MQ=255)
                                                              aCATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAGCAGAACCTATAAAAAACAAGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACAATGATTTGGGTCAAAAAATAGGAGTTTTAAAGAaatgcata                                                                 >  1:62545/1‑140 (MQ=37)
                                                              aCATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaatg                                                                 >  1:729087/1‑141 (MQ=255)
                                                                  gATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATc                                                             <  1:290175/141‑1 (MQ=255)
                                                                   aTAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCa                                                            <  1:477170/141‑1 (MQ=255)
                                                                   aTAAATGAGTGGTGAGTGATGTTTTTTTACTTTGCAAGCGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTAAAAAATTAGGGGTATTAAAAAAATGAATGAATCa                                                            <  1:826618/141‑1 (MQ=255)
                                                                     aaaTTAGCGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATa                                                          >  1:790730/1‑141 (MQ=255)
                                                                      aaTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATaa                                                         >  1:48584/1‑141 (MQ=255)
                                                                            tGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTAt                                                   >  1:642810/1‑141 (MQ=255)
                                                                               ttAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCaa                                                >  1:66629/1‑141 (MQ=255)
                                                                                tAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTACAAAAATGAATGAATCATAATTTTATCAAg                                               >  1:748619/1‑141 (MQ=255)
                                                                                tAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAg                                               >  1:748678/1‑141 (MQ=255)
                                                                                   cTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAACGGTTACTTTATATAGTATTACTGCAACCATTGATTTGGGTCAAACATTCGGGGTATTAAAAAACTGAATGAATCATACATTTCTCAAGCTg                                            >  1:831291/1‑141 (MQ=255)
                                                                                   cTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTg                                            <  1:149109/141‑1 (MQ=255)
                                                                                   cTATATCTTTTTACTTTGCAAAAGAACCTAAAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTg                                            <  1:579623/141‑1 (MQ=255)
                                                                                        tttttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGCGACCACGAGTATTGATAAAAACACTAAATGAATCATAATTTTATCAAGCTGATTgg                                       >  1:196618/1‑141 (MQ=255)
                                                                                           ttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgaga                                    >  1:14331/1‑141 (MQ=255)
                                                                                             ttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAgg                                  >  1:212996/1‑141 (MQ=255)
                                                                                             ttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATCTTATGAAGCTGATTGGAGAgg                                  <  1:445590/141‑1 (MQ=255)
                                                                                                cTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTa                               >  1:751081/1‑141 (MQ=255)
                                                                                                    gCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAAt                           >  1:67492/1‑141 (MQ=255)
                                                                                                       aCAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATCTAGTGTTACGGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCa                        <  1:495802/141‑1 (MQ=255)
                                                                                                           aaCCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtat                    >  1:327669/1‑141 (MQ=255)
                                                                                                           aaCCGATAAAAAGCATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGGCAAAAATTAGGGGTATTAAAAAAATGAATGAGTCATAGTTTTATCAAGCTGAATGGAGACGTTAAAATGCATtat                    >  1:614491/1‑141 (MQ=255)
                                                                                                            aCCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtata                   >  1:763158/1‑141 (MQ=255)
                                                                                                                 taAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATaaaa              >  1:509421/1‑141 (MQ=255)
                                                                                                                       aCATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTAtt        <  1:213274/141‑1 (MQ=255)
                                                                                                                           gggCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGt    >  1:712411/1‑141 (MQ=255)
                                                                                                                            ggCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTc   <  1:23971/141‑1 (MQ=255)
                                                                                                                             gCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAATTGCATTATATAAAATTTATTGAGTCa  <  1:730394/141‑1 (MQ=255)
                                                                                                                             gCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCa  <  1:22708/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
GATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCA  >  USA300TCH1516_ALE/76837‑77102

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: