Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,481 C→T 14.2% intergenic (+27/+19) umuC → / ← dsbB translesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,232,4810CT14.2% 56.0 / 3.9 30intergenic (+27/+19)umuC/dsbBtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base C (12/13);  new base T (2/2);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

AAGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAA  >  NC_000913/1232333‑1232621
                                                                                                                                                    |                                                                                                                                            
aaGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                                                              <  2:35637/149‑1 (MQ=255)
     ggCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGGGCTCCCggcgg                                                                                                                                         >  2:139240/1‑145 (MQ=255)
     ggCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAgg                                                                                                                                         >  2:532109/1‑149 (MQ=255)
     ggCAGAGGAACACTCTATTTTGCCGGGCAGAGGATCAAGCAACAGTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAAGCCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAgg                                                                                                                                         >  1:599453/1‑149 (MQ=255)
         gagGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                                                                     >  2:184365/1‑149 (MQ=255)
           cgAACCCTCTATGTTGCCGGGCAGGGGATCCAGCAACCATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTACACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                                                                   <  2:410697/148‑1 (MQ=255)
            gAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCg                                                                                                                                  >  1:85966/1‑149 (MQ=255)
             aaCACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGACGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                                                                 <  1:123111/149‑4 (MQ=255)
             aaCACTCTATTTCGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                                                                 <  1:678715/149‑4 (MQ=255)
                             gggCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGa                                                                                                                 >  2:449837/1‑149 (MQ=255)
                               gCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                                                               <  2:638391/149‑1 (MQ=255)
                                   gggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGGAGGCGCGCCGAATGGATTAGCGACCGAACAGAt                                                                                                           >  1:168704/1‑149 (MQ=255)
                                     gATCCAGCAACAATGGCAGATGAAGAGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCa                                                                                                         >  1:447325/1‑149 (MQ=255)
                                                 aTGGAAGAGGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGTGTCAAATAAACATATCGGCAGGAAAAAAGCGTTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCt                                                                                             <  2:603129/149‑1 (MQ=255)
                                                               gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGGCCGACCAGACACCGGTTTTTCGCTTTAAACGGCTGGGa                                                                               >  2:625609/1‑149 (MQ=255)
                                                                 gCCATGGTTCCACCATGTTATACGACGCGAAGTTTCGACTTACGCAGGGTCAAATAAATATAGCTGCAGGACAAAAGTGCTCCCTCAGGAGGGCCGAATGGATTAGCGACCGTACAGATCACGTTTTGTCGCTTTAAACGGCTGGGaaa                                                                             <  2:130484/149‑1 (MQ=255)
                                                                 gCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaaa                                                                             <  1:640396/149‑1 (MQ=255)
                                                                 gCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaaa                                                                             <  1:104783/149‑1 (MQ=255)
                                                                    aTGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                                                          <  1:184365/149‑1 (MQ=255)
                                                                          tCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAccacca                                                                    <  2:85966/149‑1 (MQ=255)
                                                                                 gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCGCCAGCACTGc                                                             >  2:417704/1‑149 (MQ=255)
                                                                                   tataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGa                                                           <  1:467616/149‑1 (MQ=255)
                                                                                                   cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGATAAGCGAt                                           >  1:177023/1‑149 (MQ=255)
                                                                                                        ttACTGAGGGTCAAATAAATATCGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaaaa                                      <  1:29244/149‑1 (MQ=255)
                                                                                                            ttaGGGCCAAATAAATATAGCGGTAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTCCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                  <  1:740701/147‑1 (MQ=255)
                                                                                                            tGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                  <  1:139240/149‑1 (MQ=255)
                                                                                                            tGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCGTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                  <  1:532109/149‑1 (MQ=255)
                                                                                                              aGGGTCAAATAAATATAGCGGCAGGAACAAAGAGTTCCGTCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGGTAAAAATACCGa                                >  1:390287/1‑149 (MQ=255)
                                                                                                                                cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCAttt              >  2:212126/1‑149 (MQ=255)
                                                                                                                                            ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGCCCTaaaa  >  1:593745/2‑149 (MQ=255)
                                                                                                                                            ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaaa  >  1:181444/2‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AAGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAA  >  NC_000913/1232333‑1232621

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