Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,994,589 (T)6→7 coding (118/1833 nt) uvrC ← excision nuclease subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,994,5831.T90.0% 57.4 / 0.5 20coding (124/1833 nt)uvrCexcision nuclease subunit C
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (8/10);  total (9/11)
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.57e-01

GATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCTTTTGCGTCAAACTGATCACTCACTGAC  >  NC_000913/1994452‑1994710
                                                                                                                                    |                                                                                                                                
gATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGt                                                                                                                                 <  2:428702/133‑1 (MQ=255)
 aTGTAGTTGTGTTCTAGCAACAGCGCTTCGGTTTCTGTGTGAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGtt                                                                                                                                <  1:350760/133‑1 (MQ=255)
      gTTGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTTCTACGGTAATAGCTCGAAAGCGGttttttt                                                                                                                          >  1:194660/1‑132 (MQ=255)
       ttGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGtttttttc                                                                                                                         <  2:286552/133‑3 (MQ=255)
                       gcgcTTCGGTTTCTGTGTGAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGc                                                                                                          <  2:16785/133‑1 (MQ=255)
                          cTTCGGTTTCTGTGTGAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGa                                                                                                       <  1:15956/133‑1 (MQ=255)
                                     gtgtgAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAAcc                                                                                            >  1:503787/1‑133 (MQ=255)
                                        tgAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTa                                                                                         >  2:108869/1‑133 (MQ=255)
                                            tAACCGTTACGTCAACTTGCTGGATCTGTGCGTCCAGCGCTGCGGTTTTGCGCGAAGCGAGGTTGCTACGGTAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTaccac                                                                                     <  2:194660/133‑1 (MQ=255)
                                              aCCGTTACATCAATTTGCTGGATCCGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTCTCAGGTCTATCGATTTACCGACATATATAATCGTAccacca                                                                                   >  2:137397/1‑133 (MQ=255)
                                                cGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGc                                                                                 >  2:266325/1‑133 (MQ=255)
                                                          aTTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCTCTTTGCCGCCATAGATAACCGTACCACCAGGATCGTACATg                                                                       >  2:142123/1‑133 (MQ=255)
                                                                        ggCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCctgg                                                         <  2:339884/133‑1 (MQ=255)
                                                                        ggCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCctgg                                                         <  1:505049/133‑1 (MQ=255)
                                                                            aCCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCctggctgg                                                     <  1:360437/133‑1 (MQ=255)
                                                                                       ggTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGt                                          <  2:393288/133‑1 (MQ=255)
                                                                                                       aGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCtttt                          <  1:266325/133‑1 (MQ=255)
                                                                                                             cTACGGAAATAGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCTTTTGCGTCa                    >  1:17361/1‑133 (MQ=255)
                                                                                                                       aGCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCTTTTGCGTCAAACTGAtcac          >  2:464762/1‑133 (MQ=255)
                                                                                                                        gCTGGAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCTTTTGCGTCAAACTGAtcact         >  2:439769/1‑133 (MQ=255)
                                                                                                                          tGGAAAGCTGTTTTTCTCAGGTTTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCTTTTGCGTCAAACTGAtcactca       <  1:142123/133‑1 (MQ=255)
                                                                                                                            gAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCTTTTGCGTCAAACTGATcactcact     <  1:108486/133‑1 (MQ=255)
                                                                                                                            gAAAGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCTTATGCTTCAAACTGATcactcact     >  1:187843/1‑133 (MQ=255)
                                                                                                                               aGCCGTTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCTTTTGCGTCAAACTGATCACTCACTGAc  <  1:464762/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                
GATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGAGTAACCGTTACATCAATTTGCTGGATCTGGGCGACCAGCGCTTCGGTTTTGCGCGAAGCGAGGTTGCTACGGAAATAGCTGGAAAGCCGTTTTTTCAGGTCTTTCGCTTTGCCGACATAGATAACCGTACCACCAGCATCGTACATGCGATAAACGCCTGGCTGGCTGGTTACGGTTTTTAAAAACGCTTTTGCGTCAAACTGATCACTCACTGAC  >  NC_000913/1994452‑1994710

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: