Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,375,713 T→A 7.8% I280I (ATT→ATA)  sbcC → Nuclease SbcCD subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,375,7130TA 7.8% 161.2 / 2.9 51I280I (ATT→ATA) sbcCNuclease SbcCD subunit C
Reads supporting (aligned to +/- strand):  ref base T (21/26);  new base A (2/2);  total (23/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 = 9.97e-01

GATTGAAAATATACCACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATATTGAAAATAGA  >  USA300TCH1516_ALE/1375575‑1375842
                                                                                                                                          |                                                                                                                                 
gattgaAAATATACCACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATAAAAAGATAATATTaa                                                                                                                                 >  1:89590/1‑141 (MQ=255)
   tgaAAATATACCACTTTTACAAGCTAGGTCGAAAGAAAGACCAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAAtca                                                                                                                              <  1:513293/141‑1 (MQ=255)
     aaaaTATACCACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAAtcaac                                                                                                                            >  1:549728/1‑141 (MQ=255)
          ataCCACTTTTACAAGCTTGTGCTAAATAAATACTAGCATCTGGAAATGAAAGTAAAGAAACGGCTATTAAAGAATTTGAAATAATAGAAGAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAAtcaactcaac                                                                                                                       <  2:298917/141‑1 (MQ=255)
          ataCCACATTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAAtcaactcaac                                                                                                                       <  1:282790/141‑1 (MQ=255)
             ccACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACaaa                                                                                                                    <  1:597379/141‑1 (MQ=255)
                     aCAAACGAGGTAGAAAGAAATAATAACATGTGTAAATGCAAGTAAAGAAACTGATATGAACGAACTTGAAATATTAGAAAAAACAACAGTAGAAAATAATATATTAAAAGATAATATTAATCAACACAAAAaaaataaaat                                                                                                            >  2:232688/1‑141 (MQ=255)
                      cAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATGATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAAtt                                                                                                           >  1:74577/1‑141 (MQ=255)
                       aaGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATCTTAAAAGATAATATAAATCAACTCAACAAAAATAAAATTg                                                                                                          >  1:512963/1‑141 (MQ=255)
                       aaGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATAAATCAACTCAACAAAAATAAAATTg                                                                                                          >  1:513026/1‑141 (MQ=255)
                         gCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGAt                                                                                                        >  1:502329/1‑141 (MQ=255)
                             ggTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCg                                                                                                    >  1:591689/1‑141 (MQ=255)
                             ggTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCg                                                                                                    >  2:582079/1‑141 (MQ=255)
                                cGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTc                                                                                                 >  1:201631/1‑141 (MQ=255)
                                  aaagaaaTACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCaa                                                                                               <  2:326911/141‑1 (MQ=255)
                                     gaacTATTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAAGATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTg                                                                                            <  1:41661/141‑1 (MQ=255)
                                      aaaTACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGa                                                                                           >  2:168674/1‑141 (MQ=255)
                                         tACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAg                                                                                        <  1:613694/141‑1 (MQ=255)
                                         tACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAg                                                                                        >  1:280181/1‑141 (MQ=255)
                                          aCTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGa                                                                                       >  1:309414/1‑141 (MQ=255)
                                                  cTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGaacaacaac                                                                               <  1:555651/141‑1 (MQ=255)
                                                   tGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAATATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGaacaacaaca                                                                              >  1:10120/1‑140 (MQ=255)
                                                   tGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAAcc                                                                              <  2:318576/141‑1 (MQ=255)
                                                      aaaTGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGa                                                                           >  2:287778/1‑141 (MQ=255)
                                                        aTGAAAGTAAAGAAACGTGTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGaaa                                                                         <  1:143367/141‑1 (MQ=255)
                                                           aaaGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATag                                                                      <  1:449022/141‑1 (MQ=255)
                                                                aaaGAAACTGCTATTAAAGGATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGa                                                                 >  1:398004/1‑141 (MQ=255)
                                                                 aaGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAACGAACAACAACCTGAAATAGAGGGaa                                                                >  2:16733/1‑141 (MQ=255)
                                                                               aaaGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTa                                                  <  1:410100/141‑1 (MQ=255)
                                                                                  aaaattGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAAATGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAg                                               <  2:123988/137‑1 (MQ=255)
                                                                                   aaTTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAgt                                              <  2:74577/141‑1 (MQ=255)
                                                                                    aTTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAgtt                                             <  2:385398/141‑1 (MQ=255)
                                                                                      ttGAAATAATATAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAgttgt                                           <  2:264487/141‑1 (MQ=255)
                                                                                         aaaTAATAGAAAACAAAACATTAGAAAATAATATATTAAAAGATAAGATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTAc                                        <  1:606725/141‑1 (MQ=255)
                                                                                               taGAAAAGAAATCAATAGAAAATAATAGATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGata                                  <  2:304231/141‑1 (MQ=255)
                                                                                                   aaagaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTAc                              <  1:171803/141‑1 (MQ=255)
                                                                                                    aagaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACa                             >  2:296829/1‑141 (MQ=255)
                                                                                                    aagaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACa                             >  2:189718/1‑141 (MQ=255)
                                                                                                    aagaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACa                             <  1:615806/141‑1 (MQ=255)
                                                                                                      gaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACaaa                           <  1:287778/141‑1 (MQ=255)
                                                                                                      gaaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACaaa                           >  1:429676/1‑141 (MQ=255)
                                                                                                       aaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAAc                          <  1:234305/141‑1 (MQ=255)
                                                                                                       aaaaCATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAAc                          <  1:163744/141‑1 (MQ=255)
                                                                                                                  aaaTAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATta               >  2:583175/1‑141 (MQ=255)
                                                                                                                   aataatATATTAAAAGATAATATAAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATtat              <  2:512963/141‑1 (MQ=255)
                                                                                                                   aataatATATTAAAAGATAATATAAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATtat              <  2:513026/141‑1 (MQ=255)
                                                                                                                   aataatATATAAAAAGATAATCTTATGCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATtat              >  1:661074/1‑141 (MQ=255)
                                                                                                                      aataTATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATAtt           <  2:251367/141‑1 (MQ=255)
                                                                                                                      aataTATTAAAAGATAATATTAATCAACTCAACAAAAATAAAAATGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATAtt           <  2:658243/141‑1 (MQ=255)
                                                                                                                            ttAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATATTGAAAAt     <  1:168674/141‑1 (MQ=255)
                                                                                                                               aaaGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATATTGAAAATaga  >  2:588704/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                 
GATTGAAAATATACCACTTTTACAAGCTAGGTCGAAAGAAATACTAGCATCTGTAAATGAAAGTAAAGAAACTGCTATTAAAGAATTTGAAATAATAGAAAAGAAAACATTAGAAAATAATATATTAAAAGATAATATTAATCAACTCAACAAAAATAAAATTGATTTCGTTCAATTGAAAGAACAACAACCTGAAATAGAGGGAATTGAAGCTAAGTTAAAGTTGTTACAAGATATTACAAACCTATTGAATTATATTGAAAATAGA  >  USA300TCH1516_ALE/1375575‑1375842

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