Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,018,190 C→T 22.5% intergenic (+115/‑27) acp → / → xynZ Sodium/proton‑dependent alanine carrier protein/Endo‑1,4‑beta‑xylanase Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,018,1900CT22.5% 34.4 / 4.7 23intergenic (+115/‑27)acp/xynZSodium/proton‑dependent alanine carrier protein/Endo‑1,4‑beta‑xylanase Z
Reads supporting (aligned to +/- strand):  ref base C (10/7);  new base T (2/3);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AAAACAAGCACGTGCCAAAAAGTAATCTACTTTTGTTTATAGTATATGTAGTGATCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATTAATAAACATGTTTTATATAGTAATATTTTAAATAGAGATGTCACGTTAAGTATTTATTTACCAGAATCTTATAATCAACTTG  >  USA300TCH1516_ALE/1018051‑1018325
                                                                                                                                           |                                                                                                                                       
aaaaCAAGCACGTGCCAAAAAGTAATCTACTTTTGTTTATAGTATATGTAGTGATCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTacaataaca                                                                                                                                        >  2:77410/1‑141 (MQ=255)
        cACGTGCCAAAAAGTAATCTACTTTTGTTTATAGTATATGTAGTGATCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTAAACTGCAAAAATATAGTGTCTCTTGGTACAATAACAATACAact                                                                                                                                >  2:61736/1‑141 (MQ=255)
                      tAAGCTACTTTTGTTTATAGTATATGTAGTGCACAATTAACTCTTCGGATAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAATAATACAACTACTAGGGGCACttt                                                                                                                  <  1:125614/141‑1 (MQ=255)
                             cTTTTGTTTATAGTATATGTAGTGATCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGt                                                                                                           >  2:265551/1‑141 (MQ=255)
                                 tGTTTATAGTATATGTAGTGATCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTACACACAGAGTTTCTCTTGGTACAATAATAATACAACTACTAGAGGCACTTTTTTATGTCAGa                                                                                                       >  2:270635/1‑141 (MQ=255)
                                       taGTATATGTAGTGATCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTaa                                                                                                 >  2:185915/1‑141 (MQ=255)
                                              tGTAGTGATCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGATCTTTGCGTTTGAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGt                                                                                          <  1:138811/141‑1 (MQ=255)
                                              agaaatGATCAATTAATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAACAATATAGTTAGTCTTCGTCTTACACGTTGACAACTACTAGGGGCACTTTACTATGTCAGAATTTAAAACTGGt                                                                                          <  2:134045/136‑1 (MQ=255)
                                                  gTGATCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATGTGAAATTTAAAAATATAGTGTATCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGa                                                                                      >  2:272538/1‑141 (MQ=255)
                                                    taTCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGAtt                                                                                    <  2:11167/140‑1 (MQ=255)
                                                       cATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATtaat                                                                                 >  2:295112/1‑141 (MQ=255)
                                                         tttGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATtaataa                                                                               <  1:260134/141‑1 (MQ=255)
                                                         tttGATAAAAAAGAAAAGTAGTGAGAATTGTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATtaataa                                                                               <  1:2477/141‑1 (MQ=255)
                                                                 aaaagaaaagTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTATACCGGGCAAGATTAATAAACATGttt                                                                       >  2:219865/1‑141 (MQ=255)
                                                                 aaaagaaaagTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACAACTAAGAGCACTATTCTATGTCAGAATTTAAAACTGGTAAGTTTAATAATCCTGttt                                                                       >  2:81997/1‑141 (MQ=255)
                                                                 aaaagaaaagTATTGAGAATTTTAGGTGCTCAGAAAAATGAATATTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATCAATAAACATGttt                                                                       >  2:281226/1‑141 (MQ=255)
                                                                        aagTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAATAATACAACTACTAGCGGCACTTTTTTATGACAGAATTTAAAACTGGTAAGATTAATAAACATGTTTtatataa                                                                >  2:157328/1‑140 (MQ=255)
                                                                        aagTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACTACTTCTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGCAATGTTAATTAACATGTTTTATATAg                                                                >  1:8265/1‑141 (MQ=255)
                                                                         agTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATTAATAAACATGTTTTATATAGt                                                               >  2:66139/1‑141 (MQ=255)
                                                                                agaATTTTAGGTGCGCAGAAATTTGAATTTTAAAAAGATAGTGTCTTTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATTAATAAACATGTTTTATATAGTAATAttt                                                        <  2:104849/141‑1 (MQ=255)
                                                                                  aaTTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATTAATAAACATGTTTTATATAGTAATATTTTa                                                      <  2:53218/141‑1 (MQ=255)
                                                                                                   aaTTTGAATTTTAAAAATACAGTGTCTCTTGGTACAACAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATTAATAAACATGTTTTATATAGTAATATTTTAAATAGAGATGTCACGtt                                     <  1:253366/141‑1 (MQ=255)
                                                                                                                             tctTGGTATAATAATAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTTGTAAGATTAATCACTGTGTTGTATATAGTAATATTTTAAATAGAGATGTCACGTTAAGTATTTATTTACCAGAATCTtata           <  2:90070/141‑1 (MQ=255)
                                                                                                                             tctTGGTACAATAATAATACAACTACTAGGGGCTCTTTTTTATGTCAGAATTTAAAACTGGTAAGATTAATAAACATGTTTTATATAGTAATATTTTAAATAGAGATGTCACGTTAAGTATTTATTTACCAGAATCTtata           <  2:254811/141‑1 (MQ=255)
                                                                                                                                      aataacaataCAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATTAATAAACATGTTTTATATAGTAATATTTTAAATAGAGATGTCACGTTAAGTATTTATTTACCAGAATCTTATAATCAACTTg  <  1:185915/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
AAAACAAGCACGTGCCAAAAAGTAATCTACTTTTGTTTATAGTATATGTAGTGATCATTTGATAAAAAAGAAAAGTATTGAGAATTTTAGGTGCTCAGAAATTTGAATTTTAAAAATATAGTGTCTCTTGGTACAATAACAATACAACTACTAGGGGCACTTTTTTATGTCAGAATTTAAAACTGGTAAGATTAATAAACATGTTTTATATAGTAATATTTTAAATAGAGATGTCACGTTAAGTATTTATTTACCAGAATCTTATAATCAACTTG  >  USA300TCH1516_ALE/1018051‑1018325

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: