Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 849,986 G→T E68D (GAG→GAT)  USA300HOU_RS04150 → Gluconeogenesis factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE849,9860GT97.9% 124.7 / NA 47E68D (GAG→GAT) USA300HOU_RS04150Gluconeogenesis factor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (12/34);  total (12/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGAGTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGACTAATATTACGAATGATTTCGGACATGCCATTAAAGC  >  USA300TCH1516_ALE/849855‑850111
                                                                                                                                   |                                                                                                                             
ttCCCAATTGATATTACGGCGTTTGTAACGTTTGCTGATAATGGTTGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATGGCAGCTTTAAGTGATTCTGATTCa                                                                                                                            >  2:460695/1‑135 (MQ=255)
ttCCCAATTGATATTACGACAAATATAACAATTGCTGAAAATAATGGAAATACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCCGCTTTAAGTGATTATGATTcg                                                                                                                            >  1:222812/1‑134 (MQ=255)
        tGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAg                                                                                                                    <  2:405634/135‑1 (MQ=255)
                  gCGATTGTAACAGTTGCTGATAATGGTGGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACttttt                                                                                                          >  2:38876/1‑135 (MQ=255)
                  gCGATTGTAACAGTTGCTGATAATGGTGGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTAATTCTGATTCAGTTTTAAGCCAACttttt                                                                                                          >  2:425908/1‑135 (MQ=255)
                   cGATTGTAACAGTTGCTGATAATGGTGGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTTAGTTTTAAGCCAACttttgc                                                                                                         <  1:229127/135‑3 (MQ=255)
                      ttGTAACAGTTGCTGATAATGGTGGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGAACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCCACTGTTTCAGt                                                                                                      >  2:561025/1‑135 (MQ=255)
                          aaCAGTTGCTGATAATGGTGGAAGCACATGTAAAATCCGAGCACATATATTTACACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCg                                                                                                  <  1:655360/135‑1 (MQ=255)
                             ggTTGCTGATAGTCGTGGGAGAAAAGGGAAAATTAGAGATGAAATGGATATACCAGCACCATGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCtt                                                                                               <  2:87689/134‑1 (MQ=255)
                              gTTGCTGATAATGGTGGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTTAGTTTTAAGCCAACTTTTTCAGTATCGCttt                                                                                              <  2:593858/135‑1 (MQ=255)
                              gTTGCTGATAATGGTGGGAGAACAAGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCttt                                                                                              >  1:291125/1‑135 (MQ=255)
                                       aaTGGTGGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAAt                                                                                     >  2:317157/1‑135 (MQ=255)
                                         tggtggGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATACAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCa                                                                                   >  2:288566/1‑135 (MQ=255)
                                            tggGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCATAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAATTCAACTTTTTCATTATCGCTTTGAAGAAAAGCAAAt                                                                                <  1:31061/135‑1 (MQ=255)
                                              ggaggaCAGGGAAAATCAGAGATGAAATGTATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAATTTTTTCAGTATCGCTTTGAAGAAAATCaaagga                                                                              <  1:277211/135‑4 (MQ=255)
                                                  tACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCgg                                                                          >  1:431405/1‑135 (MQ=255)
                                                   aCAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGt                                                                         >  2:282622/1‑135 (MQ=255)
                                                     aGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCACTTCTAAGCCAATTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCa                                                                       <  1:295183/135‑1 (MQ=255)
                                                     aGGGAAAATCAGAGAGTAAATCGATATTGCAGCACCAGGAGACATCAGAAATGTTATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAACTCAAATTAGCGGTCa                                                                       <  2:542030/135‑1 (MQ=255)
                                                                 agaTGAAATGGATATTCCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAGGTATTTGGGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTaa                                                           <  1:452897/135‑1 (MQ=255)
                                                                      aaaTGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATgtat                                                      <  1:116349/135‑1 (MQ=255)
                                                                       aaTGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAATCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATttatt                                                     <  2:542764/135‑1 (MQ=255)
                                                                        aTGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATttatta                                                    <  1:589627/135‑1 (MQ=255)
                                                                        aTGGAGATCCCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTAGCGCTTTTAAGAATATCAAATTAGCGGTCACTCATTAGGGAATttatta                                                    <  2:90018/135‑1 (MQ=255)
                                                                            ttataCCAGCACCAGGAGACATCAGTAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATTAAATTAGCGGTTACTCATTAGGTAATTTATTCAtag                                                <  1:459677/134‑3 (MQ=255)
                                                                            ctataCCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGTTCACTCATTAGGTAATTTATTCAttg                                                <  1:459838/134‑3 (MQ=255)
                                                                            atatACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTACCTGCGCATATAATTTATTAATCg                                                <  2:642514/135‑1 (MQ=255)
                                                                            atacACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCATTCATTAGGTAATTTATTAATCg                                                <  1:88461/135‑1 (MQ=255)
                                                                                  cAGCACCATGAGACATCATAAAAGTGATCTGAGCTTGAAGTGATTCTGATCCAGTTCTAAGCCTACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGTTCACTCATTAGGTAATTTATTAATCGCAGGTa                                          <  2:214047/135‑1 (MQ=255)
                                                                                         aGGAGAAATCACAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTTAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGACTaa                                   <  2:109681/135‑1 (MQ=255)
                                                                                              aCATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGACTAATATTa                              >  1:234970/1‑135 (MQ=255)
                                                                                               cATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATTACTAATCTTAc                             >  2:453157/1‑135 (MQ=255)
                                                                                                    gAAATGTGATTGCAGCTTTAACTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTATCAGTCACTCATTAGCTAATCTAATAATCGCAGGTATGACTAATATTACGAATg                        <  2:264859/135‑1 (MQ=255)
                                                                                                     aaaTGTGATTTCAGCTTTAAGTGACTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGACTAATATTACGAAtaa                       <  2:173630/135‑3 (MQ=255)
                                                                                                     aaaTGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGACTAATATTACGAATGa                       <  1:604531/135‑1 (MQ=255)
                                                                                                     aaaTGTGATTGCAGCTTTAAGTGATGCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGTATAATATTACGAATGa                       <  1:347821/135‑1 (MQ=255)
                                                                                                      aaTGTGATTGCTATTTTTAGTGCTTCTGATTCAGTTTTCAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGTGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGACTAATATTACGAATGAt                      <  2:333503/135‑1 (MQ=255)
                                                                                                      aaTGTGATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATTGCAGGTATGACTAATATTACGAATGAt                      <  1:81309/135‑1 (MQ=255)
                                                                                                          tgATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGCTAATTTATTAATCGCAGGTATGACTAATATTACGAATGATTTCg                  <  2:144711/135‑1 (MQ=255)
                                                                                                          tgATTGCAGCTTTAAGTGATTCTGAGTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATCAGCGGTCACTCATTAGGTAATTTATCAATCGCTGCTTCCACCAATATTACGAATGATTTct                  <  1:460695/135‑2 (MQ=255)
                                                                                                           gATTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGACTAATATTACGAATGATTTCgg                 <  1:282622/135‑1 (MQ=255)
                                                                                                            aTTGCAGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGACTAATATTACGAATGATTTCGGa                <  1:657250/135‑1 (MQ=255)
                                                                                                             ttGCAGCTTTAAGTGATTCTGATTCAGATTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGATTAATATTACGAATGATTTCGGAc               <  1:317157/135‑1 (MQ=255)
                                                                                                                 aGCTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGATTAATATTACGAATGATTTCGGACATGc           <  1:38876/135‑1 (MQ=255)
                                                                                                                   cTTTAAGTGATTCTGATTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCTCTCATTATGTAATTTATTAATCGCAGGTATGATTAATATTACGAATGATTTCGGACATGCCa         <  1:30507/135‑1 (MQ=255)
                                                                                                                     ttAAGTGATTCTGATTCAGTTTTAAGCCAACTTTCTCAGGATCGGTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGTAGGTATGACTAATATTACGAATGATTTCGGACATGCCAtt       <  1:231098/135‑1 (MQ=255)
                                                                                                                          tGATTCTGCTTCAGTTTTAAGCCAACTTTTTCAGTTTCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGAATAATATTACGACTGATTTTGTACAGGCCATTAAAGc  <  2:115349/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                             
TTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGAGTACAGGGAAAATCAGAGATGAAATGGATATACCAGCACCAGGAGACATCAGAAATGTGATTGCAGCTTTAAGTGATTCTGAGTCAGTTTTAAGCCAACTTTTTCAGTATCGCTTTGAAGAAAATCAAATTAGCGGTCACTCATTAGGTAATTTATTAATCGCAGGTATGACTAATATTACGAATGATTTCGGACATGCCATTAAAGC  >  USA300TCH1516_ALE/849855‑850111

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: ‑