Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,230,806 A→G T14A (ACT→GCT)  umuD → DNA polymerase V protein UmuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,230,8060AG96.1% 171.6 / ‑4.5 51T14A (ACT→GCT) umuDDNA polymerase V protein UmuD
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (32/17);  total (33/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGT  >  NC_000913/1230670‑1230940
                                                                                                                                        |                                                                                                                                      
tATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCtt                                                                                                                                     <  2:560848/140‑1 (MQ=255)
     aGCGACAAATATTGATAGCCTGCATCAGTTTTGATCTACTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCg                                                                                                                                <  1:314452/140‑1 (MQ=255)
       cGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCt                                                                                                                              <  1:49272/140‑1 (MQ=255)
        gACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTa                                                                                                                             >  2:436999/1‑140 (MQ=255)
            aaTATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTa                                                                                                                         >  2:523823/1‑140 (MQ=255)
             atatTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAg                                                                                                                        <  2:615257/140‑1 (MQ=255)
                     aGCCTGAATCAGTATTGATCTGCTGGCAAGATCAGTCTAATGTATATATAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTg                                                                                                                <  2:207331/140‑1 (MQ=255)
                      gCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGt                                                                                                               >  2:82740/1‑140 (MQ=255)
                                 tATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCt                                                                                                    >  2:338504/1‑140 (MQ=255)
                                           ctgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCAcc                                                                                          >  2:93040/1‑140 (MQ=255)
                                                             tGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGa                                                                        >  1:518807/1‑140 (MQ=255)
                                                                  ataAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCGCCAGCAGCACATTACGTTGAACAgc                                                                   >  1:314806/1‑140 (MQ=255)
                                                                    aaaaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAgcgc                                                                 >  1:344122/1‑140 (MQ=255)
                                                                    aaaaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAgcgc                                                                 >  2:491907/1‑140 (MQ=255)
                                                                     aaaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCa                                                                <  1:436999/140‑1 (MQ=255)
                                                                      aaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCAt                                                               >  2:63358/1‑140 (MQ=255)
                                                                      aaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCAt                                                               <  2:518807/140‑1 (MQ=255)
                                                                         cAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGa                                                            >  1:505128/1‑140 (MQ=255)
                                                                         cAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGa                                                            >  1:605230/1‑140 (MQ=255)
                                                                         cAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGa                                                            >  2:408760/1‑140 (MQ=255)
                                                                           gTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATc                                                          <  1:338504/140‑1 (MQ=255)
                                                                           gTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATc                                                          >  2:577313/1‑140 (MQ=255)
                                                                           gTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATc                                                          >  2:460867/1‑140 (MQ=255)
                                                                            tataACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCt                                                         <  1:460867/140‑1 (MQ=255)
                                                                                  ttCAGGCAGATTATTATTTTGTCTATCAAGCCTGCTGATATCCGCGAACTTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCa                                                   <  2:199048/140‑1 (MQ=255)
                                                                                      ggCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGTAGATTACATTGTACGGCGCATTGATCTGAATCAACTg                                               >  1:192862/1‑140 (MQ=255)
                                                                                      ggCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTg                                               >  2:165598/1‑140 (MQ=255)
                                                                                      ggCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTg                                               >  2:137325/1‑140 (MQ=255)
                                                                                      ggCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTg                                               >  2:663547/1‑140 (MQ=255)
                                                                                      ggCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCCGAATCAACTg                                               >  1:606603/1‑140 (MQ=255)
                                                                                      ggCAGATTATTATGTTGCTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTg                                               >  2:47568/1‑140 (MQ=255)
                                                                                           attattATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCa                                                   >  2:619977/1‑131 (MQ=255)
                                                                                           attattATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCa                                                   <  1:619977/131‑1 (MQ=255)
                                                                                                     tgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATccc                                >  1:541077/1‑140 (MQ=255)
                                                                                                     tgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATccc                                >  1:335048/1‑140 (MQ=255)
                                                                                                     tgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATccc                                >  1:565802/1‑140 (MQ=255)
                                                                                                       tttATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAg                              >  1:577228/1‑140 (MQ=255)
                                                                                                        ttATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAgc                             <  2:442804/140‑1 (MQ=255)
                                                                                                         tATCAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGAGCAACGATCTGAATCAGCTGTCGATCCAGCATCCCAgcg                            >  2:271808/1‑140 (MQ=255)
                                                                                                            cAAGCCTGCGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGa                         <  1:165598/140‑1 (MQ=255)
                                                                                                                    cGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTc                 >  1:528424/1‑140 (MQ=255)
                                                                                                                    cGGATCTCCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTCCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTc                 >  2:131116/1‑140 (MQ=255)
                                                                                                                    cGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTCCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTc                 >  2:501899/1‑140 (MQ=255)
                                                                                                                       aTCTCCGCGAAATTGTGACTTTACCGCTATTTAGCAATCATGTTCAGTATGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTc              <  2:169973/140‑1 (MQ=255)
                                                                                                                        tctcCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCa             <  2:528424/140‑1 (MQ=255)
                                                                                                                         ctcCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTc                 >  1:515153/1‑135 (MQ=255)
                                                                                                                         ctcCGCGAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTc                 <  2:515153/135‑1 (MQ=255)
                                                                                                                              cgAAATTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCaa       >  2:372225/1‑140 (MQ=255)
                                                                                                                                aaaTTGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgt     >  2:348321/1‑140 (MQ=255)
                                                                                                                                   ttGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgtggt  <  2:605230/140‑1 (MQ=255)
                                                                                                                                   ttGTGGCTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgtggt  <  1:491907/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGT  >  NC_000913/1230670‑1230940

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