Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,081,153 T→C 6.8% N295N (AAT→AAC purM → Phosphoribosylformylglycinamidine cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,081,1530TC 6.8% 169.3 / 3.6 59N295N (AAT→AACpurMPhosphoribosylformylglycinamidine cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base T (34/21);  new base C (2/2);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

TATGAAAATATTCCACGTGCATTGCCAGCCGGATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCAACTGCAATTGAATTGTTGGGGGT  >  USA300TCH1516_ALE/1081016‑1081293
                                                                                                                                         |                                                                                                                                            
tATGAAAATATTCCACGTGCATTGCCAGCCGGATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATg                                                                                                                                           <  1:95578/141‑1 (MQ=255)
tATGAAAATATTCCACGTGCATTGCCAGCCGGATATGCTGCTAGAATTGATACTACAGCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATg                                                                                                                                           <  1:301321/141‑1 (MQ=255)
tATGAAAATATTCCACCTGACTTGCCAGCCGGATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATg                                                                                                                                           <  2:48338/141‑1 (MQ=255)
    aaaaTATTCCACGTGCATTGCCAGCCGGATATGCTGCTAGAATTGCTACGCCGTCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGCCACAACTGAAATGTATAACATTTTTAATATgggat                                                                                                                                       >  1:148255/1‑139 (MQ=255)
     aaaTATTCCACGTGCATTGCCAGCCGGATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTAt                                                                                                                                      <  1:645001/141‑1 (MQ=255)
           tCCACGTGCATTGCCAGCCGGATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCta                                                                                                                                >  1:106498/1‑141 (MQ=255)
            ccACGTGCATTGCCAGCCGGATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCtat                                                                                                                               >  1:690295/1‑141 (MQ=255)
             cACGTGCATTGCCAGCCGGATATGCTGCTAATATTGATACCACATCATCTCCAACACCATACATAGTTGATTGGTTACCACAACAAGGCAATATAGACACAAATTAAAGGGATAACATTTTGAATATGCGGATTGGCgatt                                                                                                                              <  2:188216/141‑2 (MQ=255)
               cGTGCATTGCCAGCCGGATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACg                                                                                                                            <  2:65060/141‑1 (MQ=255)
                          aGCCGGATATGCTGCTAGAATTGAAACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTtgg                                                                                                                 >  1:491115/1‑140 (MQ=255)
                           gCCGGATATGCTGCTAGAATTGATACTACACCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGCTATTGGCTACACGGTTATCGTtgat                                                                                                                <  1:385042/141‑1 (MQ=255)
                            ccGGATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTtgatg                                                                                                               <  1:640372/141‑1 (MQ=255)
                              ggATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGaa                                                                                                             >  1:442412/1‑141 (MQ=255)
                              ggATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAACATGGGTATTGGCTATACGGTTATCGTTGATGaa                                                                                                             >  1:180359/1‑141 (MQ=255)
                              ggATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAACATGGGTATTGGCTATACGGTTATCGTTGATGaa                                                                                                             >  1:180287/1‑141 (MQ=255)
                                atatGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGaaaa                                                                                                           >  1:668414/1‑141 (MQ=255)
                                      tgctAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGt                                                                                                     <  2:745318/141‑1 (MQ=255)
                                            atttGATACTACATCATTTCCAACTCCAGAAATATTTGATTGGTGACAGTAACGGGGCGTTATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCAcg                                                                                               <  2:79964/139‑1 (MQ=255)
                                              ttGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCAcgcc                                                                                             >  1:651462/1‑140 (MQ=255)
                                                 aTACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCCCgcgcgt                                                                                          >  2:254919/1‑139 (MQ=255)
                                                  tactacATCATTTCCAACACCAAAAATATTTGATTGGTTACAAGAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCttt                                                                                         <  1:141977/141‑1 (MQ=255)
                                                  tactacATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCttt                                                                                         >  2:119906/1‑141 (MQ=255)
                                                  tactacATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCttt                                                                                         >  2:119901/1‑141 (MQ=255)
                                                       catcatTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAAATGTATCACGGGCTTTGAAGa                                                                                    >  2:239013/1‑141 (MQ=255)
                                                             ttCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTGTGAAGATTTTAg                                                                              >  1:526068/1‑141 (MQ=255)
                                                                cAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATTAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGcca                                                                           >  2:730960/1‑139 (MQ=255)
                                                                cAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAg                                                                           >  1:285045/1‑141 (MQ=255)
                                                                 aaCACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGa                                                                          >  1:703072/1‑141 (MQ=255)
                                                                    acCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAACATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACa                                                                       <  2:180359/141‑1 (MQ=255)
                                                                    acCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAACATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACa                                                                       <  2:180287/141‑1 (MQ=255)
                                                                     ccAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTCAGCGTGGAGGAAAAAGATGTAGCACGCGCTTTGAAGATTTTAGCAGAACaa                                                                      >  1:629793/1‑141 (MQ=255)
                                                                        aaaaTATATAATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAAt                                                                   <  2:557623/141‑1 (MQ=255)
                                                                                tGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGc                                                           >  1:562841/1‑141 (MQ=255)
                                                                                tGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGc                                                           >  1:14733/1‑141 (MQ=255)
                                                                                tGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGATTTGAAGGTTTTCGCGGAACGAAATGTGTAAGc                                                           >  1:203756/1‑141 (MQ=255)
                                                                                 gATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGcc                                                          >  2:529284/1‑141 (MQ=255)
                                                                                 gATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGcc                                                          >  2:509310/1‑141 (MQ=255)
                                                                                   ttGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTa                                                        <  1:321100/141‑1 (MQ=255)
                                                                                     ggTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATc                                                      <  1:275234/141‑1 (MQ=255)
                                                                                     ggTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATc                                                      <  2:690295/141‑1 (MQ=255)
                                                                                      gTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCa                                                     <  1:662570/141‑1 (MQ=255)
                                                                                            acaacaaGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTgg                                               >  1:103883/1‑141 (MQ=255)
                                                                                            acaacaaGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTgg                                               <  1:301062/141‑1 (MQ=255)
                                                                                                caaGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCat                                           <  2:372748/141‑1 (MQ=255)
                                                                                                caaGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATAACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCat                                           >  2:527815/1‑141 (MQ=255)
                                                                                                    gCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGGATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATGTTAGCAGAACAAAATGTGGAGACCTATCAAATTGGTCatatat                                       >  1:158542/1‑139 (MQ=255)
                                                                                                         ataGACACAAATGAAATGTATAACATTTTTAATATGAGTATTGGCTATACGGTTATCGTTGATTAAAAAGATGTATCAAGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCATACCAAATTGGTAATATcgttcaa                                  >  2:668361/1‑137 (MQ=255)
                                                                                                            gACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCGTATCAAATTGGTCATATTGTGAAAAAt                               >  1:585417/1‑141 (MQ=255)
                                                                                                              cacaAATGAAAGGTATAACATTTTTAATATGGGTATTGCCTATACGGTTATGGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGCTCATATTGTGAAAAATGa                             >  2:351477/1‑141 (MQ=255)
                                                                                                                   atgaaatgTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCaa                        >  2:325685/1‑141 (MQ=255)
                                                                                                                    tgaaatgTATAACATTTTTAATATGGGACCTGGCGATACGGTACACAACCGTGATAAAGGTGTATCACGCGCTTTGAAGATTTTAGCCGGACAAAAAATGGCAACCAACCCCATGGGTCATATGGCGACAAATggcgcaac                       >  2:658451/1‑134 (MQ=255)
                                                                                                                     gaaatgTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCaagt                      >  2:97792/1‑139 (MQ=255)
                                                                                                                      aaatgTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAAGGAGTCAACTg                     >  2:473511/1‑141 (MQ=255)
                                                                                                                          atataAGATTTTTTATATGCGTATTGGTTATACGGTTATCGTTGATTAAAAAGATGTAGCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCAACTGCaat                 <  1:310792/140‑1 (MQ=255)
                                                                                                                           tataACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCAACTGCaatt                <  2:285045/141‑1 (MQ=255)
                                                                                                                              aaCATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTGAACTGCaattgaa             >  2:705473/1‑141 (MQ=255)
                                                                                                                              aaCATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCAACTGCaattgaa             <  2:651462/141‑1 (MQ=255)
                                                                                                                               aCATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTCTCAACTTGTTCATATTGTGAAAAATGAGTCAACTGCaattgaat            >  1:590622/1‑141 (MQ=255)
                                                                                                                                       aaTATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCAACTGCAGTTGAATTGTTgggg    >  1:232871/1‑141 (MQ=255)
                                                                                                                                        atatGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCAACTGCAATTGAATTGTTggggg   <  2:153913/141‑1 (MQ=255)
                                                                                                                                         tatGGGTATTGGCTATACGGTTATCGTTGATAAAAAAGATGTATCACGCGCTTCGAAGCTTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCAACTGCAATTGAATTGTTGGGGGt  >  1:523459/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                            
TATGAAAATATTCCACGTGCATTGCCAGCCGGATATGCTGCTAGAATTGATACTACATCATTTCCAACACCAAAAATATTTGATTGGTTACAACAACAAGGCAATATAGACACAAATGAAATGTATAACATTTTTAATATGGGTATTGGCTATACGGTTATCGTTGATGAAAAAGATGTATCACGCGCTTTGAAGATTTTAGCAGAACAAAATGTGGAAGCCTATCAAATTGGTCATATTGTGAAAAATGAGTCAACTGCAATTGAATTGTTGGGGGT  >  USA300TCH1516_ALE/1081016‑1081293

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