Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,324,256 (A)9→8 coding (1028/2010 nt) mutL → DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,324,2480A.88.6% 155.9 / 15.9 44coding (1020/2010 nt)mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base . (15/24);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGTAAATGTTCATCCAACAAAACTAGAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACGATGA  >  USA300TCH1516_ALE/1324110‑1324391
                                                                                                                                          |                                                                                                                                               
agACGTAAATGTTCATCCAACAAAACTAGAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATGAATTCCTAAATATAACTTGGATTATGTGCCaaaaaaaata                                                                                                                                        <  3:373592/148‑3 (MQ=255)
   cGTAAATGTTCATCCAACAAAACTAGAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAg                                                                                                                                    >  3:84666/1‑148 (MQ=255)
          gTTCATCCAACAAAACTAGAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGAATTATGTGCCAAAAAAAATAAAg                                                                                                                                    >  1:3306/1‑141 (MQ=255)
            tCATCCAACAAAACTAGAAGTGCGTGTAGGAGATGAAGAGCAACTATACCAATAGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTAcatt                                                                                                                           >  4:190691/1‑148 (MQ=255)
                  aacaaAACTAGAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCTAAAAAAATAAAGTATTAcat                                                                                                                            <  1:48452/141‑1 (MQ=255)
                       aaCTAGAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGaacaa                                                                                                                >  3:41050/1‑148 (MQ=255)
                            gAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACaaaa                                                                                                           >  3:396908/1‑148 (MQ=255)
                                   gTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACaaaa                                                                                                           >  2:27614/1‑141 (MQ=255)
                                     ttATCAAAAGAAGAGCAACTATACCAATTGATTGTAAGCAAAATACAAAAAGCATTTAAAGACCGAATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTAATACATTCATTCGAACAACAAAAAATCGAAtt                                                                                                  <  3:391307/148‑1 (MQ=255)
                                       aTCAAAAGAAGAGCAACTTTACCAATTGATAGGGAGCAAAATACAAGAAGCATTTAAAGACAGTATATTAATTTCTAAAAATAACTTGGAAAAAAAGAAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTg                                                                                                <  4:161152/148‑1 (MQ=255)
                                           aaaGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGTATTATGTACCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACa                                                                                            <  3:389199/148‑1 (MQ=255)
                                             agaagaGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACaaa                                                                                          <  4:174731/148‑1 (MQ=255)
                                               aagaGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACATCAAAAAATCGAATTTGAACAAAGa                                                                                        <  4:165546/148‑1 (MQ=255)
                                               aagaGCAAATATAACAATTGATTGTGACAAAAATACAAGAAGCATTTAAAGACCGTATATTAATACATAAAAATAACTTGGATTATATGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGa                                                                                        <  4:79300/148‑1 (MQ=255)
                                               aagaGAAACAATAACAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACAGTATATTAATTCCTAAAAATAACTTGGAAAAAAAAAAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGa                                                                                        <  4:160173/148‑1 (MQ=255)
                                                  agCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACaa                                                                                     <  3:60575/148‑1 (MQ=255)
                                                  agCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACaa                                                                                     <  3:58260/148‑1 (MQ=255)
                                                        tataCCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAAcaca                                                                               <  3:286173/148‑1 (MQ=255)
                                                         ataCCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACag                                                                              <  4:143640/148‑1 (MQ=255)
                                                         ataCCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCATATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACag                                                                              <  4:269953/148‑1 (MQ=255)
                                                                ttgattgTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATaa                                                                       >  3:164570/1‑148 (MQ=255)
                                                                ttgattgTGAGAAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATAAAAAAAAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATaa                                                                       <  4:27870/148‑1 (MQ=255)
                                                                   attgTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACCGAGAATAACCa                                                                    <  4:333843/148‑1 (MQ=255)
                                                                        gAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAgaga                                                               >  3:154726/1‑148 (MQ=255)
                                                                                 aCAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGtttt                                                      >  4:337733/1‑148 (MQ=255)
                                                                                     gaagCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTAAGAGCAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATc                                                  <  4:41820/148‑1 (MQ=255)
                                                                                        gCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCAAAAAAAAATAAAGTAATACATTAATTCGAACAACAAAAAATCGAATATGAACAAAGACAAAACACAGAGAAAAACCAAGAAAAGACGattt                                                      >  2:53862/1‑141 (MQ=255)
                                                                                           tttAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTgaaga                                            <  3:246268/148‑1 (MQ=255)
                                                                                           tttAAAGACCGTATATCAATGCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTgaaga                                            <  3:127817/148‑1 (MQ=255)
                                                                                               aaGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGt                                        <  3:190691/148‑1 (MQ=255)
                                                                                                      tataTTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGt                                        <  1:27614/141‑1 (MQ=255)
                                                                                                      tataTTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTa                                 <  4:396908/148‑1 (MQ=255)
                                                                                                           tAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCa                            <  4:59827/148‑1 (MQ=255)
                                                                                                              ttCCTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCAttt                         <  3:284519/148‑1 (MQ=255)
                                                                                                              atCCTAAAAATAACTAGGAATATGAGCCAAAAAAAATAAAGTATTACATTCATTCGAACAAAAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCAttt                         <  4:2264/147‑1 (MQ=255)
                                                                                                                 cTAAAAATAACTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATg                      <  4:154354/148‑1 (MQ=255)
                                                                                                                   aaaaaTAACTTGGATTATGTGCCAAAAAAAATAAAGTATTGCATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAACACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGga                    >  3:216915/1‑148 (MQ=255)
                                                                                                                           cTTGGATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATc            >  3:81988/1‑148 (MQ=255)
                                                                                                                              ggATTATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCaaa         >  3:152162/1‑148 (MQ=255)
                                                                                                                                 ttATGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACg      <  4:364009/148‑1 (MQ=255)
                                                                                                                                   aTGTGCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCa           <  2:51984/141‑1 (MQ=255)
                                                                                                                                     gtgCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCATAACGAAga  >  3:105874/1‑148 (MQ=255)
                                                                                                                                     gtgCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACGATga  >  3:153487/1‑148 (MQ=255)
                                                                                                                                     gtgCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACGATga  >  4:47705/1‑148 (MQ=255)
                                                                                                                                     gtgCCAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACGATga  <  3:114620/148‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
AGACGTAAATGTTCATCCAACAAAACTAGAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCCAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACGATGA  >  USA300TCH1516_ALE/1324110‑1324391

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: