Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,232,002 C→T R273W (CGG→TGG)  umuC → translesion error‑prone DNA polymerase V subunit; DNA polymerase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,232,0020CT98.0% 182.5 / ‑4.3 49R273W (CGG→TGG) umuCtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (22/26);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATA  >  NC_000913/1231859‑1232149
                                                                                                                                               |                                                                                                                                                   
cGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGggcggc                                                                                                                                              >  1:168949/1‑151 (MQ=255)
  gTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGcggcgg                                                                                                                                            >  1:692067/1‑151 (MQ=255)
      gTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGaaa                                                                                                                                         >  1:668979/1‑150 (MQ=255)
       tGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGaaaa                                                                                                                                        >  2:599983/1‑150 (MQ=255)
             gcgcGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCg                                                                                                                                  <  1:561083/150‑1 (MQ=255)
             gcgcGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGc                                                                                                                                 <  2:324325/151‑1 (MQ=255)
              cgcgGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGc                                                                                                                                 >  2:40578/1‑150 (MQ=255)
                cggcgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGc                                                                                                                              <  1:576532/151‑1 (MQ=255)
                 ggcgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGggcgg                                                                                                                                               <  2:639918/133‑1 (MQ=255)
                 ggcgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGggcgg                                                                                                                                               >  1:639861/1‑133 (MQ=255)
                  gcgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCg                                                                                                                             <  2:119798/150‑1 (MQ=255)
                   cgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAg                                                                                                                           <  2:453183/151‑1 (MQ=255)
                      aCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCAt                                                                                                                        >  2:552823/1‑151 (MQ=255)
                      aCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCAt                                                                                                                        <  2:4279/151‑1 (MQ=255)
                           gTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAata                                                                                                                   <  1:598291/151‑1 (MQ=255)
                              tGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATAtt                                                                                                                 <  2:122073/150‑1 (MQ=255)
                                  aCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCgg                                                                                                            <  1:574019/151‑1 (MQ=255)
                                     ggAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGttt                                                                                                         <  2:361425/151‑1 (MQ=255)
                                       aaGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAt                                                                                                       <  1:403242/151‑1 (MQ=255)
                                        agagTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTTGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAtc                                                                                                      <  1:596744/151‑1 (MQ=255)
                                           gTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTcc                                                                                                   <  2:479392/151‑1 (MQ=255)
                                           gTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTcc                                                                                                   <  2:705958/151‑1 (MQ=255)
                                                  ccGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTtta                                                                                            <  2:440823/151‑1 (MQ=255)
                                                   cgacgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTttat                                                                                           <  2:160505/151‑1 (MQ=255)
                                                        aaGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTttatta                                                                                         >  1:371915/1‑148 (MQ=255)
                                                                   tATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCAttt                                                                           >  2:545771/1‑151 (MQ=255)
                                                                    aTCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTg                                                                          >  2:408099/1‑151 (MQ=255)
                                                                        gTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCt                                                                      >  2:283072/1‑151 (MQ=255)
                                                                           cccGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCa                                                                    <  1:318194/150‑1 (MQ=255)
                                                                             cGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATg                                                                 >  1:680832/1‑151 (MQ=255)
                                                                              gcccGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCCGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAGCAACTTCGCAGCGACCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAAt                                                                  >  1:627624/4‑149 (MQ=255)
                                                                                 cGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAcc                                                             <  2:761146/151‑1 (MQ=255)
                                                                                    ttGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCtt                                                           >  2:693239/1‑150 (MQ=255)
                                                                                         gAACGCATCACGGATTAGCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAGGACGTCACCATTTGCGCTCAATGAACCTTATTAc                                                      <  1:320910/150‑1 (MQ=255)
                                                                                          aaCGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACgg                                                    <  1:511166/151‑1 (MQ=255)
                                                                                            cGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGc                                                   >  1:388626/1‑150 (MQ=255)
                                                                                               aTCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATa                                               <  2:163200/151‑1 (MQ=255)
                                                                                                  aCGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTACCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAgcg                                            <  1:704978/151‑1 (MQ=255)
                                                                                                         aTCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGt                                     >  1:597775/1‑151 (MQ=255)
                                                                                                            cGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTaa                                   >  2:146165/1‑150 (MQ=255)
                                                                                                             gTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTa                                    >  1:559681/1‑148 (MQ=255)
                                                                                                                        aGGCCATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGcc                      <  1:411525/151‑1 (MQ=255)
                                                                                                                           ccATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCACTATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGccc                     >  1:292595/1‑149 (MQ=255)
                                                                                                                           ccATTTGTAGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCAc                   >  1:267783/1‑151 (MQ=255)
                                                                                                                                   aGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGAc            >  2:126860/1‑150 (MQ=255)
                                                                                                                                   aGTTACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGAc            >  2:724398/1‑150 (MQ=255)
                                                                                                                                      tACGCTGCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGAcagc         <  1:633779/150‑1 (MQ=255)
                                                                                                                                          cTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAggg     >  1:416498/1‑150 (MQ=255)
                                                                                                                                            gCCTGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGata  <  1:552771/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATA  >  NC_000913/1231859‑1232149

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