Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 76,969 G→A 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,9690GA100.0% 68.6 / NA 21intergenic (+509/‑107)USA300HOU_RS00140/USA300HOU_RS00350hypothetical protein/Monoacylglycerol lipase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/16);  total (5/16)

TCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAAG  >  USA300TCH1516_ALE/76831‑77106
                                                                                                                                          |                                                                                                                                         
tctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAagg                                                                                                                                         >  1:450610/1‑141 (MQ=255)
 ctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGt                                                                                                                                        <  1:157514/141‑1 (MQ=255)
    cAGATCTCCGGATTTCCGCCATGCCCCGAAATTAGCAGCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTaa                                                                                                                                     <  1:265144/141‑1 (MQ=17)
       aTCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAAttt                                                                                                                                  <  1:15538/141‑1 (MQ=255)
         cTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTtat                                                                                                                                <  1:333702/141‑1 (MQ=255)
               aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTa                                                                                                                          >  1:100371/1‑141 (MQ=255)
                                       catcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTc                                                                                                  <  1:292310/141‑1 (MQ=255)
                                           aagCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCaaaa                                                                                              <  1:163516/139‑1 (MQ=255)
                                                          gCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaa                                                                               <  1:240573/141‑1 (MQ=255)
                                                             aaCACTGACATGATAAATTAGTGGTTAGCTATATTGTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaaa                                                                            <  1:195010/141‑1 (MQ=255)
                                                               cacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaat                                                                          <  1:345400/141‑1 (MQ=255)
                                                               cacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaat                                                                          <  1:81873/141‑1 (MQ=255)
                                                                             aTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTTTATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAt                                                            <  1:285989/141‑1 (MQ=255)
                                                                             aTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAt                                                            <  1:82837/141‑1 (MQ=255)
                                                                                   ggTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATc                                                      <  1:260375/141‑1 (MQ=255)
                                                                                        gCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCt                                                 <  1:310523/141‑1 (MQ=255)
                                                                                                           cAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATg                              >  1:210246/1‑141 (MQ=255)
                                                                                                             aCAGAACCTATAAAAACAATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCa                            <  1:185513/141‑1 (MQ=255)
                                                                                                                   ccTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtatat                      <  1:185977/141‑1 (MQ=255)
                                                                                                                       taAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATaaaa                  >  1:383824/1‑141 (MQ=255)
                                                                                                                             aCATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGCCCAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTAtt            >  1:416592/1‑141 (MQ=255)
                                                                                                                                  ggCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTc       >  1:95260/1‑141 (MQ=255)
                                                                                                                                       aaaaGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAAg  <  1:418101/141‑1 (MQ=255)
                                                                                                                                       aaaaGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAAg  <  1:418155/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
TCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAAG  >  USA300TCH1516_ALE/76831‑77106

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: