Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,565,939 A→T 13.0% I533N (ATT→AAT)  mutL ← DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,565,9390AT13.0% 80.5 / 3.2 30I533N (ATT→AAT) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (16/10);  new base T (2/2);  total (18/12)
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

ACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGC  >  NC_002947/5565792‑5566083
                                                                                                                                                   |                                                                                                                                                
aCGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATc                                                                                                                                                <  2:307844/150‑1 (MQ=255)
      ccGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGATTCGCCAga                                                                                                                                          <  1:435607/150‑2 (MQ=255)
                             cagcTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCaa                                                                                                                   <  1:386948/150‑1 (MQ=255)
                                   aTTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGaa                                                                                                             >  1:225739/1‑150 (MQ=255)
                                              ccTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCt                                                                                                  <  2:139182/150‑1 (MQ=255)
                                                       gATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAg                                                                                                                                           >  2:586299/1‑100 (MQ=255)
                                                       gATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAg                                                                                                                                           <  1:586299/100‑1 (MQ=255)
                                                       gATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACt                                                                                         >  2:478869/1‑150 (MQ=255)
                                                             tGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCAt                                                                                   >  1:116274/1‑150 (MQ=255)
                                                                   cGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTTCTTGAGCAAGACCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCCCTGCAATTCGAAGCCCAGGCGCTGAAACTACTGTGCATGCTCtt                                                                             >  1:556999/1‑150 (MQ=255)
                                                                        tCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGca                                                                        >  1:322436/1‑150 (MQ=255)
                                                                          cATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGcaca                                                                      >  1:286054/1‑150 (MQ=255)
                                                                           aTAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCAc                                                                                          <  2:446506/129‑1 (MQ=255)
                                                                           aTAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCAc                                                                                          >  1:446506/1‑129 (MQ=255)
                                                                                        gCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACgct                                                        >  1:528404/1‑150 (MQ=255)
                                                                                                tGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAgg                                                <  2:225739/150‑1 (MQ=255)
                                                                                                            ttGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTcc                                    >  2:280875/1‑150 (MQ=255)
                                                                                                              ggCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCgg                                  >  2:531978/1‑150 (MQ=255)
                                                                                                                    aGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGACCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCCCTTCGGCACAGTCGGCTTCACGCTGGCTGATGGCGAGGGCCTACGGCCCcag                            >  2:493191/1‑150 (MQ=255)
                                                                                                                      ccTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACcagca                          >  2:79910/1‑150 (MQ=255)
                                                                                                                        tGCTTGAGCAAGGCCAGGATATGGCGAAATGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGTGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCCTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCagcagc                        <  2:556999/150‑1 (MQ=255)
                                                                                                                                       gggATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCCAGCCGCtgctg         <  1:480428/150‑5 (MQ=255)
                                                                                                                                        ggATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCCAGCCGCtgctgg        <  1:79910/150‑6 (MQ=255)
                                                                                                                                          cccTGGCGATTCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTGCACACCTGATGAGTCTGCACTCTGGCTGCGTGCCCGCGTCACCCGCAATAGCGGGGGCTGGTCGCTCAGGccc      >  1:318776/3‑150 (MQ=255)
                                                                                                                                           ccTGGCGATTCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAgg                                          >  2:168358/2‑113 (MQ=255)
                                                                                                                                           ccTGGCGATTCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAgg                                          <  1:168358/112‑1 (MQ=255)
                                                                                                                                             tGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCg   >  1:62317/1‑150 (MQ=255)
                                                                                                                                              ggCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGc  >  2:453467/1‑150 (MQ=255)
                                                                                                                                               gCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAgg                                          <  2:203609/109‑1 (MQ=255)
                                                                                                                                               gCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAgg                                          >  1:203609/1‑109 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGC  >  NC_002947/5565792‑5566083

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