Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 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
*CP0007301,324,2480A.89.5% 139.5 / 12.1 38coding (1020/2010 nt)mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base . (13/21);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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