Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,231,396 G→T 19.0% R14I (AGA→ATA)  stp → Serine/threonine phosphatase stp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,231,3960GT19.0% 118.1 / 16.0 55R14I (AGA→ATA) stpSerine/threonine phosphatase stp
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base T (5/6);  total (27/28)
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

ATTAATGCCACAATACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGT  >  USA300TCH1516_ALE/1231257‑1231532
                                                                                                                                           |                                                                                                                                        
aTTAATGCCACAATACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATaga                                                                                                                                         <  2:425005/141‑1 (MQ=255)
  tAATGCCACAATACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAtaga                                                                                                                                       <  1:323956/141‑1 (MQ=255)
       ccACAATACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGa                                                                                                                                  >  2:725243/1‑141 (MQ=255)
           aaTACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAAtaa                                                                                                                              <  2:281888/141‑3 (MQ=255)
              aCGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGa                                                                                                                           <  1:528025/141‑1 (MQ=255)
                gtcgtGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGTTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATg                                                                                                                         <  1:525858/141‑1 (MQ=255)
                     gAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGggt                                                                                                                    <  1:72650/141‑1 (MQ=255)
                       aCAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGggtgg                                                                                                                  <  1:517935/141‑1 (MQ=255)
                       aCAAGGTTCGGATAGTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGggtgg                                                                                                                  <  2:639788/141‑1 (MQ=255)
                          aGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTAt                                                                                                               >  1:37037/1‑141 (MQ=255)
                           ggTTCGGATATTGACGCAGCTTGTGGTCAATTAAGCGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGGATGAAGATGCGGGTGGTAtt                                                                                                              <  1:212958/141‑1 (MQ=255)
                                gggatgtGACGCAGCTTGGGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTta                                                                                                         <  2:60553/135‑1 (MQ=255)
                                gggagttGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTta                                                                                                         <  2:229655/136‑1 (MQ=255)
                                ggctatTGACGCAGCTTGTGGTCAATTAAGAGGAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTta                                                                                                         <  2:509099/138‑1 (MQ=255)
                                ggATATTGACGCTGCTTGTGGTCGATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACTATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTta                                                                                                         <  2:229927/141‑1 (MQ=255)
                                 gATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTtat                                                                                                        <  1:643734/141‑1 (MQ=255)
                                    atTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAAt                                                                                                     >  1:236510/1‑141 (MQ=255)
                                        aCGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCaaa                                                                                                 >  1:726982/1‑141 (MQ=255)
                                         cGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAAc                                                                                                >  1:378390/1‑141 (MQ=255)
                                           cAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTa                                                                                              <  1:677789/141‑1 (MQ=255)
                                            aGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTaa                                                                                             >  1:77991/1‑141 (MQ=255)
                                             ccTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAAt                                                                                            <  1:390791/140‑1 (MQ=255)
                                                tgtgGTCAATTAAGAGCAAAGGGACGACAAGTAGACACGAGGTCAAGATAAATGCTAGCGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATCATCAAAGTAATcaa                                                                                         >  2:188302/1‑141 (MQ=255)
                                                 gtgGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaac                                                                                        >  1:404951/1‑141 (MQ=255)
                                                 ctgGTCAATGAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAAGGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaac                                                                                        <  2:656793/140‑1 (MQ=255)
                                                   ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaacaa                                                                                      >  2:706768/1‑141 (MQ=255)
                                                   ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATACTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaacaa                                                                                      >  1:521562/1‑141 (MQ=255)
                                                             gagCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTc                                                                            >  1:147868/1‑141 (MQ=255)
                                                              agCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCt                                                                           >  1:164699/1‑141 (MQ=255)
                                                                      aaCGACAAGTAGAAACACGCAAGAGGTTAATGTTTGCGGCACAATATTTCGTTGAGTCTGTCCACCATAGAGTCATGTATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATg                                                                   <  1:372310/141‑1 (MQ=255)
                                                                       aCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTTTGAttg                                                                  <  2:682528/141‑3 (MQ=255)
                                                                       aCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATgg                                                                  >  1:198958/1‑141 (MQ=255)
                                                                       aCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATgg                                                                  >  1:449959/1‑141 (MQ=255)
                                                                        cGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGt                                                                 >  1:709379/1‑141 (MQ=255)
                                                                        cGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGt                                                                 >  1:348248/1‑141 (MQ=255)
                                                                         caCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCCGTGTGATGGTa                                                                <  1:733744/140‑1 (MQ=255)
                                                                              gTAGAAACGAGGTAAAAACCAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGGATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATgggg                                                           >  1:441947/1‑140 (MQ=255)
                                                                                aGAAATGAGGTAAAGACAAATGCTAGAGGCCAAAATTTTATGGGGTTCCGGACATCATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTgg                                                         <  2:754922/141‑1 (MQ=255)
                                                                                   aaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAAGTTTTAGTTCTGTGTGATGGTATGGGCGGGCa                                                      >  1:292675/1‑141 (MQ=255)
                                                                                   aaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCa                                                      >  1:39283/1‑141 (MQ=255)
                                                                                    aCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCAt                                                     <  2:688366/141‑1 (MQ=255)
                                                                                    aCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTGAGTTCTGTGTGATGGTATGGGTGGCCAt                                                     <  1:567333/141‑1 (MQ=255)
                                                                                            aaGACAAATGCAGTAGGCACAATTTTTGACCGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATCTGTTATAATCAAACTAATCAACAACTTTTAGATCTGTGTGATGGTATGGGTGGCCATAAAGCacg                                             <  2:217281/141‑3 (MQ=255)
                                                                                               aCAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgaga                                          >  1:135133/1‑141 (MQ=255)
                                                                                                    tGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAGCGAGTCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTg                                     <  1:738975/141‑1 (MQ=255)
                                                                                                     gCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAACCAGGATAAGTTGc                                    >  2:287437/1‑141 (MQ=255)
                                                                                                      cTAGAGGCACAATTTTTTTCTTTTACTGGACAACAAATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                   >  2:503824/1‑141 (MQ=255)
                                                                                                        agagGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCTTAAATCAGGAGAAGTTGCAAg                                 <  1:405285/141‑1 (MQ=255)
                                                                                                         gcggCACAATTTTTTCCTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGt                                <  1:92872/139‑1 (MQ=255)
                                                                                                          agGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTa                               >  2:139033/1‑141 (MQ=255)
                                                                                                           ggCACAATTTGTTACTGATACGGAACAACATATAAATAAGAATGAAGATGCGGGTAGTATTTTTTATAAGCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTaa                              <  2:772516/141‑1 (MQ=255)
                                                                                                                  ttttttACTAATACTGGACAACATAGAGATAAGAATAAAGAAGCGGGTGGTATTTTTTCTAATCAAACTAATCAACAACTTTTAGTTCGGTGTGATGGTATGAGTGGCCATAAAGAGGAAGAAGTTGCAGATTAATTTGtt                       >  2:339339/1‑141 (MQ=255)
                                                                                                                     tttACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACa                    >  1:191975/1‑141 (MQ=255)
                                                                                                                          tGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATgg               >  1:320551/1‑140 (MQ=255)
                                                                                                                           gATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAg              >  1:183831/1‑141 (MQ=255)
                                                                                                                                 ggACAACAAAGGGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTTTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGaaa        <  2:275265/141‑1 (MQ=255)
                                                                                                                                    caacaaAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATcc     <  2:649914/141‑1 (MQ=255)
                                                                                                                                       aatatAGATTAGGATGAAGGTGCTGGTTGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCCTAAAGCGGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGt  <  2:539923/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
ATTAATGCCACAATACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGT  >  USA300TCH1516_ALE/1231257‑1231532

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: