Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 629,681 C→A L84I (CTA→ATA)  USA300HOU_RS02990 → 3‑hexulose‑6‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE629,6810CA83.3% 118.7 / 23.8 54L84I (CTA→ATA) USA300HOU_RS029903‑hexulose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (7/2);  new base A (16/29);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 2.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

ACAACGAAGGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATT  >  USA300TCH1516_ALE/629541‑629819
                                                                                                                                            |                                                                                                                                          
aCAACGAAGGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAc                                                                                                                                            >  1:800576/1‑141 (MQ=255)
        ggTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATATGGCGCGGATGTAATTACAATAATAGGTGtt                                                                                                                                    >  1:620167/1‑141 (MQ=255)
          tttACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGc                                                                                                                                  >  1:330797/1‑141 (MQ=255)
          tttACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGc                                                                                                                                  <  1:442023/141‑1 (MQ=255)
            tACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCag                                                                                                                                <  1:581453/141‑1 (MQ=255)
             aCCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAACTTGGCGCGGATGTAATTACAATAATAGGTGTTGCaga                                                                                                                               >  1:944185/1‑141 (MQ=255)
              ccAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCagaa                                                                                                                              >  1:647857/1‑141 (MQ=255)
                agcagTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCagaaga                                                                                                                            >  1:340279/1‑141 (MQ=255)
                      tAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGAAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATc                                                                                                                      >  1:417285/1‑141 (MQ=255)
                             atGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTaaa                                                                                                               <  1:924606/141‑1 (MQ=255)
                             atGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATGACAATAATAGGTGTTGCAGAAGATGCATCAATTaaa                                                                                                               <  2:258733/141‑1 (MQ=255)
                                   gACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCt                                                                                                         <  2:294237/141‑1 (MQ=255)
                                       acaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTg                                                                                                     >  2:353680/1‑141 (MQ=255)
                                             gTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGaagaag                                                                                               >  1:65377/1‑141 (MQ=255)
                                                aTGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTc                                                                                            <  2:309318/141‑1 (MQ=255)
                                                aTGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATCGAAGAAGCTc                                                                                            <  1:239120/141‑1 (MQ=255)
                                                  gTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCAt                                                                                          <  2:5391/141‑1 (MQ=255)
                                                   tAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATa                                                                                         <  2:87883/141‑1 (MQ=255)
                                                    aaaaGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATaa                                                                                        <  2:32314/141‑1 (MQ=255)
                                                               cAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAAc                                                                             >  2:92567/1‑141 (MQ=255)
                                                                aGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACa                                                                            <  2:277381/141‑1 (MQ=255)
                                                                    aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTa                                                                        >  1:5754/1‑141 (MQ=255)
                                                                    aTGAAAATTATAGATGCAGCTAATTATGAAGTTAGACAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGAGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTa                                                                        <  1:75143/141‑1 (MQ=255)
                                                                     tGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTAc                                                                       <  1:202537/141‑1 (MQ=255)
                                                                          aTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGtt                                                                  <  1:685394/141‑1 (MQ=255)
                                                                           ttATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAAATTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTg                                                                 <  2:365118/141‑1 (MQ=255)
                                                                             ataGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGat                                                               >  1:577372/1‑141 (MQ=255)
                                                                             ataGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGat                                                               >  1:362179/1‑141 (MQ=255)
                                                                               aGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGatat                                                             <  2:254775/141‑1 (MQ=255)
                                                                                gATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATg                                                            >  2:236654/1‑141 (MQ=255)
                                                                                  tGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGAt                                                          >  2:411345/1‑141 (MQ=255)
                                                                                   gCAGCTGATTATGAAGTTAGCCAAACAATTAAATTTAGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGAtt                                                         <  1:902179/141‑1 (MQ=255)
                                                                                       cTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTg                                                     <  1:481065/141‑1 (MQ=255)
                                                                                        tGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGt                                                    >  2:254039/1‑141 (MQ=255)
                                                                                         gATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGtt                                                   <  2:119015/141‑1 (MQ=255)
                                                                                          aTTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTc                                                  >  1:728906/1‑141 (MQ=255)
                                                                                             aTGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAg                                               <  2:250570/141‑1 (MQ=255)
                                                                                             aTGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAg                                               <  2:235/141‑1 (MQ=255)
                                                                                                 aGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGAttt                                           >  1:754074/1‑141 (MQ=255)
                                                                                                  gTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTa                                          >  1:783765/1‑141 (MQ=255)
                                                                                                  gTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTa                                          >  1:886740/1‑141 (MQ=255)
                                                                                                     aGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGaa                                       >  2:213472/1‑141 (MQ=255)
                                                                                                         aaGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAAc                                   >  1:60237/1‑141 (MQ=255)
                                                                                                            cAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTg                                >  1:770069/1‑141 (MQ=255)
                                                                                                            cAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTg                                <  1:798458/141‑1 (MQ=255)
                                                                                                              aTTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCa                              <  1:236957/141‑1 (MQ=255)
                                                                                                               ttAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCaa                             <  2:14796/141‑1 (MQ=255)
                                                                                                                 aaaTTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCaaaa                           <  2:41948/141‑1 (MQ=255)
                                                                                                                     ttGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAAc                       <  2:3937/141‑1 (MQ=255)
                                                                                                                                  aaTTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATggg          <  2:6606/141‑1 (MQ=255)
                                                                                                                                    ttACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTg        <  1:538469/141‑1 (MQ=255)
                                                                                                                                    ttACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTg        <  1:636717/141‑1 (MQ=255)
                                                                                                                                      aCAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCt      <  2:413/141‑1 (MQ=255)
                                                                                                                                       cAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTg     <  1:770306/141‑1 (MQ=255)
                                                                                                                                          tACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAGAAGAACTAGATGAAATGGGTGCTGAtt  >  1:35833/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                          
ACAACGAAGGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATT  >  USA300TCH1516_ALE/629541‑629819

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: