Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 880,093 G→T intergenic (+271/‑70) emp → / → nuc Extracellular matrix protein‑binding protein emp/Thermonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE880,0930GT96.4% 181.4 / 0.6 57intergenic (+271/‑70)emp/nucExtracellular matrix protein‑binding protein emp/Thermonuclease
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (19/35);  total (20/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

AACGATTCATGTTGTAGGCTTTTTAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTGTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGTTTCAATATTACTTATA  >  USA300TCH1516_ALE/879959‑880222
                                                                                                                                       |                                                                                                                                 
aaCGATTCATGTTGTAGGCTTTTTAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTaaaaa                                                                                                                             >  1:451882/1‑141 (MQ=255)
  cGATTCATGTTGTAGGCTTTTTAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTtaaaaat                                                                                                                           >  3:228559/1‑139 (MQ=255)
   gATTCATGTTGTAGGCTTTTTAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGt                                                                                                                          <  1:92915/141‑1 (MQ=255)
        aTGTTGTATACTTTTTAATTTTACAAATAAGGCTAAATATATAAGTTCTTGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAAt                                                                                                                     <  4:63042/141‑1 (MQ=255)
         tgttgtAGGCTTTTTAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAAtt                                                                                                                    <  1:170438/141‑1 (MQ=255)
         tgttgtAGGCTTTTTAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAAtt                                                                                                                    >  2:613165/1‑141 (MQ=255)
               aGGCTTTTTAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTAtaatt                                                                                                              <  1:199053/141‑1 (MQ=255)
                    ttttAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTATAAAAATGTAGAATTATAattaatta                                                                                                         <  2:56404/141‑1 (MQ=255)
                       tAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCAACTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACg                                                                                                      >  1:372714/1‑141 (MQ=255)
                          ttttACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGttt                                                                                                   <  2:408649/141‑1 (MQ=255)
                             tACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTGTAATTAATTAACGTTTAat                                                                                                <  2:73234/141‑1 (MQ=255)
                             tACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAat                                                                                                <  2:271478/141‑1 (MQ=255)
                                     aGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAAt                                                                                        >  3:236142/1‑141 (MQ=255)
                                                 taAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAaaaga                                                                            >  3:46684/1‑141 (MQ=255)
                                                 taAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAaaaga                                                                            >  3:46686/1‑141 (MQ=255)
                                                 taAGTTCTGGCACCTAAAATATAGAAAATACATAAAAATAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAAAAAAAAAAaaaga                                                                            <  2:611102/141‑1 (MQ=255)
                                                  aaGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATAAACGTTTAATATTAAAATTAACTAAaaagaa                                                                           >  1:15295/1‑141 (MQ=255)
                                                     ttCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAaaagaaag                                                                         <  4:228559/141‑1 (MQ=255)
                                                      tCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAagag                                                                       <  2:483141/141‑1 (MQ=255)
                                                          gCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGgtgt                                                                   <  1:29672/141‑1 (MQ=255)
                                                          gCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGgtgt                                                                   <  1:465901/141‑1 (MQ=255)
                                                          gCAACTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGgtgt                                                                   <  4:219247/141‑1 (MQ=255)
                                                              cTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTgttagt                                                               <  3:155516/141‑1 (MQ=255)
                                                               aaaaaTATAGAAAATACATAAAAGTAGGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTgttagtt                                                              <  2:556893/140‑1 (MQ=255)
                                                                 aaaTATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTAt                                                            >  4:222418/1‑141 (MQ=255)
                                                                 aaaTATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTAt                                                            <  4:169444/141‑1 (MQ=255)
                                                                  aaTATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATg                                                           <  3:279064/141‑1 (MQ=255)
                                                                  aaTATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATg                                                           <  3:299196/141‑1 (MQ=255)
                                                                   atatAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGa                                                          <  4:268463/141‑1 (MQ=255)
                                                                       aGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGa                                                      >  1:179158/1‑141 (MQ=255)
                                                                       aGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGa                                                      <  3:262552/141‑1 (MQ=255)
                                                                       aGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAAGTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGa                                                      <  3:83290/141‑1 (MQ=255)
                                                                        gAAAATACATAAAAGTACGCATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGaa                                                     <  2:140536/141‑1 (MQ=255)
                                                                        gAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAATAGAATGTGGTGTTAGTTATGACAGaa                                                     <  2:451882/141‑1 (MQ=255)
                                                                         aaaaTACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAAt                                                    <  3:100087/141‑1 (MQ=255)
                                                                           aaTACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTTTTTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAACAGGTGTTAGTTATGAAAGAATaa                                                  >  2:224204/1‑140 (MQ=255)
                                                                               cATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACttat                                              <  1:178026/141‑1 (MQ=255)
                                                                                  aaaaGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTaa                                           >  1:144363/1‑141 (MQ=255)
                                                                                       taagtaTAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCt                                      <  1:376463/141‑1 (MQ=255)
                                                                                         agtaTAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTTCAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTgg                                    >  3:252660/1‑141 (MQ=255)
                                                                                          gtaTAGTTATTTTATTATAATTATTAAATTTTTATTAATTGATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGc                                   >  2:606169/1‑141 (MQ=255)
                                                                                                     ttattaTAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGc                        <  2:182/141‑1 (MQ=255)
                                                                                                     ttattaTAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATGAACTAAAAGGAGAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGc                        <  2:241883/141‑1 (MQ=255)
                                                                                                     ttattaTAATTATTAAATTTTTATTAATTAATTGTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGc                        <  2:55967/141‑1 (MQ=255)
                                                                                                        ttatAATTATTAAATTTTTATTAATTAATTTTACAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAAAGGAGTTAGTTATGACAGAATACTTATTAATTGCTGGCATATGTATGGCaaa                     >  2:591519/1‑140 (MQ=255)
                                                                                                        ttatAATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAAt                     >  1:322720/1‑141 (MQ=255)
                                                                                                           taATTATTAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGt                  <  1:64020/141‑1 (MQ=255)
                                                                                                             attattAAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGttt                >  2:401425/1‑141 (MQ=255)
                                                                                                                  taAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGATGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGTTTCAata           <  3:155659/141‑1 (MQ=255)
                                                                                                                  taAATTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGTTTCAata           <  1:19771/141‑1 (MQ=255)
                                                                                                                     aTTTTTATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGTTGGCATATGTATGGCAATCGTTTCAATATTa        <  2:446908/141‑1 (MQ=255)
                                                                                                                      tttttATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGTTTCAATATTAc       >  2:508114/1‑141 (MQ=255)
                                                                                                                       ttttATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGTTTCAATATTACt      <  3:222418/141‑1 (MQ=255)
                                                                                                                        tttATTAATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGTTTCAATATTACtt     <  2:35579/141‑1 (MQ=255)
                                                                                                                         ttattaATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATCGTCATCGCTTCAATATTACTta    <  2:179158/141‑1 (MQ=255)
                                                                                                                          tattaATTAATTTTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGTTTCAATATTACTtat   >  3:104334/1‑141 (MQ=255)
                                                                                                                           attaattaATTGTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGTTTCAATATTACTtata  >  1:510911/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                 
AACGATTCATGTTGTAGGCTTTTTAATTTTACAAATAAGGCTAAATATATAAGTTCTGGCACCTAAAATATAGAAAATACATAAAAGTAAGTATAGTTATTTTATTATAATTATTAAATTTTTATTAATTAATTGTAAAAATGTAGAATTATAATTAATTAACGTTTAATATTAAAATTAACTAAAAAGAAAGAGGTGTTAGTTATGACAGAATACTTATTAAGTGCTGGCATATGTATGGCAATCGTTTCAATATTACTTATA  >  USA300TCH1516_ALE/879959‑880222

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