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)
   cGTAAATGTTCATCCAACAAAACTAGAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAg                                                                                                                                    >  1:84666/1‑148 (MQ=255)
            tCATCCAACAAAACTAGAAGTGCGTGTAGGAGATGAAGAGCAACTATACCAATAGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTAcatt                                                                                                                           >  2:190691/1‑148 (MQ=255)
                       aaCTAGAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGaacaa                                                                                                                >  1:41050/1‑148 (MQ=255)
                            gAAGTGCGTTTATCAAAAGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACaaaa                                                                                                           >  1:396908/1‑148 (MQ=255)
                                     ttATCAAAAGAAGAGCAACTATACCAATTGATTGTAAGCAAAATACAAAAAGCATTTAAAGACCGAATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTAATACATTCATTCGAACAACAAAAAATCGAAtt                                                                                                  <  1:391307/148‑1 (MQ=255)
                                       aTCAAAAGAAGAGCAACTTTACCAATTGATAGGGAGCAAAATACAAGAAGCATTTAAAGACAGTATATTAATTTCTAAAAATAACTTGGAAAAAAAG‑AAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTg                                                                                                <  2:161152/148‑1 (MQ=255)
                                           aaaGAAGAGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGTATTATGTACC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACa                                                                                            <  1:389199/148‑1 (MQ=255)
                                             agaagaGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACaaa                                                                                          <  2:174731/148‑1 (MQ=255)
                                               aagaGCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACATCAAAAAATCGAATTTGAACAAAGa                                                                                        <  2:165546/148‑1 (MQ=255)
                                               aagaGCAAATATAACAATTGATTGTGACAAAAATACAAGAAGCATTTAAAGACCGTATATTAATACATAAAAATAACTTGGATTATATGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGa                                                                                        <  2:79300/148‑1 (MQ=255)
                                               aagaGAAACAATAACAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACAGTATATTAATTCCTAAAAATAACTTGGA‑AAAAAAAAAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGa                                                                                        <  2:160173/148‑1 (MQ=255)
                                                  agCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACaa                                                                                     <  1:60575/148‑1 (MQ=255)
                                                  agCAACTATACCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACaa                                                                                     <  1:58260/148‑1 (MQ=255)
                                                        tataCCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAAcaca                                                                               <  1:286173/148‑1 (MQ=255)
                                                         ataCCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACag                                                                              <  2:143640/148‑1 (MQ=255)
                                                         ataCCAATTGATTGTGAGCAAAATACAAGAAGCATTTAAAGACCATATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACag                                                                              <  2:269953/148‑1 (MQ=255)
                                                                ttgattgTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATaa                                                                       >  1:164570/1‑148 (MQ=255)
                                                                ttgattgTGAGAAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGA‑TAAAAAAAAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATaa                                                                       <  2:27870/148‑1 (MQ=255)
                                                                   attgTGAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACCGAGAATAACCa                                                                    <  2:333843/148‑1 (MQ=255)
                                                                        gAGCAAAATACAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAgaga                                                               >  1:154726/1‑148 (MQ=255)
                                                                                 aCAAGAAGCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGtttt                                                      >  2:337733/1‑148 (MQ=255)
                                                                                     gaagCATTTAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTAAGAG‑CAAAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATc                                                  <  2:41820/148‑1 (MQ=255)
                                                                                           tttAAAGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTgaaga                                            <  1:246268/148‑1 (MQ=255)
                                                                                           tttAAAGACCGTATATCAATGCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTgaaga                                            <  1:127817/148‑1 (MQ=255)
                                                                                               aaGACCGTATATTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGt                                        <  1:190691/148‑1 (MQ=255)
                                                                                                      tataTTAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTa                                 <  2:396908/148‑1 (MQ=255)
                                                                                                           tAATTCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCa                            <  2:59827/148‑1 (MQ=255)
                                                                                                              ttCCTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCAttt                         <  1:284519/148‑1 (MQ=255)
                                                                                                              atCCTAAAAATAACTAGGAATATGAGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAAAAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCAttt                         <  2:2264/147‑1 (MQ=255)
                                                                                                                 cTAAAAATAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATg                      <  2:154354/148‑1 (MQ=255)
                                                                                                                   aaaaaTAACTTGGATTATGTGCC‑AAAAAAAATAAAGTATTGCATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAACACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGga                    >  1:216915/1‑148 (MQ=255)
                                                                                                                           cTTGGATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATc            >  1:81988/1‑148 (MQ=255)
                                                                                                                              ggATTATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCaaa         >  1:152162/1‑148 (MQ=255)
                                                                                                                                 ttATGTGCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACg      <  2:364009/148‑1 (MQ=255)
                                                                                                                                     gtgCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCATAACGAAga  >  1:105874/1‑148 (MQ=255)
                                                                                                                                     gtgCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACGATga  >  1:153487/1‑148 (MQ=255)
                                                                                                                                     gtgCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACGATga  >  2:47705/1‑148 (MQ=255)
                                                                                                                                     gtgCC‑AAAAAAAATAAAGTATTACATTCATTCGAACAACAAAAAATCGAATTTGAACAAAGACAAAACACAGAGAATAACCAAGAGAAGACGTTTTCATCTGAAGAAAGTAACAGTAAGCCATTTATGGAAGAAAATCAAAACGATga  <  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: ‑