Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,397,834 G→A V141V (GTG→GTA)  mutL → DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,397,8340GA97.1% 111.6 / NA 34V141V (GTG→GTA) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (14/19);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGG  >  NC_000913/4397714‑4397952
                                                                                                                        |                                                                                                                      
tATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTATAACCGGCCGCGCATCCCGTAGGGACGACGCTg                                                                                                            >  2:197024/1‑133 (MQ=255)
  tCAGTTCGTTTTCCCGCCCGACTCTCACTTTACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGa                                                                                                          <  1:197024/133‑1 (MQ=255)
  tCAGTTCGGTTTCCCTCCTGACGCTCACTTCAGGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGa                                                                                                          <  2:246776/133‑1 (MQ=255)
  tCAGTTCGGTCTCCCGCCTGACGCTCACGTCATGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGa                                                                                                          <  2:186120/133‑1 (MQ=255)
           tttCCCGCCTGACGCTCATTTCTCGCACCGTAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGa                                                                                                 <  1:279805/133‑1 (MQ=255)
           tttCCCGCCTGACGCTCAATTCACGCACCGAAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGa                                                                                                 <  1:271206/133‑1 (MQ=255)
            ttCCCGCCTGACGCTGACTTCACGCACCCCAGAACAGCAGGCAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGAt                                                                                                <  2:211678/133‑1 (MQ=255)
            ttCCCGCCTGACGCTCGCTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGACTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGAt                                                                                                <  1:526914/133‑1 (MQ=255)
            ttCCCGCCTGACGCTCACTTCACGCATCTCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGAt                                                                                                <  2:268122/133‑1 (MQ=255)
             tCCCGCCTTACGCTCGCTTCACCCACCGCACAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATc                                                                                               <  1:302305/133‑1 (MQ=255)
               ccGCCTGACGCTCACTTCGCGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTg                                                                                             <  2:18539/133‑1 (MQ=255)
                cGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGAGGCGCATCGTGTGGGGACGACGCTGGAGGTGCTGGATCTGt                                                                                            >  2:160804/1‑133 (MQ=255)
                      aCGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCTGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTaca                                                                                      <  2:449309/133‑1 (MQ=255)
                                 aCGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGTCGACGCTGGAGGAGCTCGATCTGTTCTACACCACCCcggcg                                                                           >  2:249864/1‑133 (MQ=255)
                                      ccGCAGAACAGCAGGAAGCCTGGTAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCcggcgcggcg                                                                      >  1:462500/1‑133 (MQ=255)
                                      ccGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCCGTTCTACAACACCCcggcgcggcg                                                                      >  2:223765/1‑133 (MQ=255)
                                        gCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCa                                                                    >  1:545922/1‑133 (MQ=255)
                                                 cagGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTgc                                                           >  1:346733/1‑133 (MQ=255)
                                                      aGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCAcc                                                      >  2:417338/1‑133 (MQ=255)
                                                      aGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCAcc                                                      >  2:346416/1‑133 (MQ=255)
                                                      aGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGACGGAAATTCATGCGCAcc                                                      >  1:262701/1‑133 (MQ=255)
                                                          tGGTAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCgaga                                                  <  2:462500/133‑1 (MQ=255)
                                                           ggCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGTGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGaa                                                 >  2:574653/1‑133 (MQ=255)
                                                             cAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGaaaa                                               <  1:44738/133‑1 (MQ=255)
                                                                  cTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGaa                                          <  1:417338/133‑1 (MQ=255)
                                                                  cTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGaa                                          <  2:442724/133‑1 (MQ=255)
                                                                     tGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAAttt                                       >  1:343455/1‑133 (MQ=255)
                                                                          aaGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACca                                  <  1:603148/133‑1 (MQ=255)
                                                                           aGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACcac                                 <  1:412714/133‑1 (MQ=255)
                                                                           aGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACcac                                 >  2:373195/1‑133 (MQ=255)
                                                                                gcgATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATtga                            <  1:346416/133‑1 (MQ=255)
                                                                                              aCGGTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATccgcc              >  1:590772/1‑133 (MQ=255)
                                                                                                 gTAAAACCGGCGGCGCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCa           <  1:130354/133‑1 (MQ=255)
                                                                                                          gcggcgCATCCTGTAGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTgg  >  2:403124/1‑133 (MQ=255)
                                                                                                                        |                                                                                                                      
TATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAACACCCCGGCGCGGCGCAAATTCCTGCGCACCGAGAAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGG  >  NC_000913/4397714‑4397952

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