Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,484 A→G 17.3% intergenic (+30/+16) 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,4840AG17.3% 51.2 / 6.2 29intergenic (+30/+16)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 A (11/13);  new base G (3/2);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

GGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAA  >  NC_000913/1232338‑1232621
                                                                                                                                                  |                                                                                                                                         
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)
                                                                                                                                                  |                                                                                                                                         
GGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAA  >  NC_000913/1232338‑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: