Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,134 A→G 11.8% H103R (CAC→CGC)  yhcG → DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,1340AG11.8% 144.7 / 10.2 51H103R (CAC→CGC) yhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base A (22/23);  new base G (3/3);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

ATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGT  >  NC_000913/3368011‑3368223
                                                                                                                           |                                                                                         
aTGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTgg                                                                            >  3:958836/1‑139 (MQ=255)
aTGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTgg                                                                            <  4:958836/139‑1 (MQ=255)
              aTCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACaa                                                                                   >  3:866184/1‑118 (MQ=255)
              aTCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACaa                                                                                   <  4:866184/118‑1 (MQ=255)
                             aaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGacac                                            <  3:681612/142‑1 (MQ=255)
                             aaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGacac                                            >  4:681612/1‑142 (MQ=255)
                              aGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAaca                                                                  <  3:1073255/119‑1 (MQ=255)
                              aGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAaca                                                                  >  4:1073255/1‑119 (MQ=255)
                              aGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTcc                                                    >  2:821619/1‑133 (MQ=255)
                              aGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTcc                                                    <  1:821619/133‑1 (MQ=255)
                                           aaGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTc                                                     >  1:1021037/1‑119 (MQ=255)
                                           aaGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTc                                                     <  2:1021037/119‑1 (MQ=255)
                                             ggTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTgg                                                                            >  2:1004757/1‑94 (MQ=255)
                                             ggTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTgg                                                                            <  1:1004757/94‑1 (MQ=255)
                                               tataAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTGGTGaa                                                                        >  3:56477/1‑96 (MQ=39)
                                               tataAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCGCCAGACATTGTCTGGTGaa                                                                        <  4:56477/96‑1 (MQ=39)
                                                ataAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCaa                    <  1:518549/147‑1 (MQ=255)
                                                ataAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCaa                    >  2:518549/1‑147 (MQ=255)
                                                    gCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCa                                                                                           <  1:147822/72‑1 (MQ=255)
                                                    gCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCa                                                                                           >  2:147822/1‑72 (MQ=255)
                                                         ggTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGa                                                                         >  4:446673/1‑85 (MQ=255)
                                                         ggTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGa                                                                         <  3:446673/85‑1 (MQ=255)
                                                         ggTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCAtt                                                             >  4:111345/1‑97 (MQ=255)
                                                         ggTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCAtt                                                             <  3:111345/97‑1 (MQ=255)
                                                            tctcTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTc                                                     <  2:583574/102‑1 (MQ=255)
                                                            tctcTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTc                                                     >  1:583574/1‑102 (MQ=255)
                                                               ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCTGCAGACGGTGTCTGCGGAATTCCCaatggtg                     <  1:874906/131‑7 (MQ=37)
                                                               ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCTGCAGACGGTGTCTGCGGAATTCCCaatggtg                     >  2:874906/1‑125 (MQ=37)
                                                                 gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGa                                                                                      >  2:221650/1‑64 (MQ=255)
                                                                 gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGa                                                                                      <  1:221650/64‑1 (MQ=255)
                                                                       aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTc                                        >  3:132803/1‑104 (MQ=255)
                                                                       aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTc                                        >  3:489103/1‑104 (MQ=255)
                                                                       aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTc                                        <  4:489103/104‑1 (MQ=255)
                                                                       aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTc                                        <  4:132803/104‑1 (MQ=255)
                                                                                  aaTGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCg                                          <  3:981524/91‑1 (MQ=255)
                                                                                  aaTGAGGCGTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCg                                          >  4:981524/1‑91 (MQ=255)
                                                                                    tGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTg                                      <  1:754148/93‑1 (MQ=255)
                                                                                    tGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTg                                      >  2:754148/1‑93 (MQ=255)
                                                                                    tGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAAtt                               <  2:355158/100‑1 (MQ=255)
                                                                                    tGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAAtt                               >  1:355158/1‑100 (MQ=255)
                                                                                           tttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGc                                                                           <  2:636293/49‑1 (MQ=255)
                                                                                           tttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGc                                                                           >  1:636293/1‑49 (MQ=255)
                                                                                                   tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTc                              >  4:827266/1‑86 (MQ=255)
                                                                                                   tcatttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTc                              <  3:827266/83‑1 (MQ=255)
                                                                                                    ctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTg             >  4:329944/1‑102 (MQ=255)
                                                                                                    ctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTg             <  3:329944/102‑1 (MQ=255)
                                                                                                             aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAac                                                                   <  2:174525/39‑1 (MQ=255)
                                                                                                                   tGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGAtt               >  4:921126/1‑85 (MQ=255)
                                                                                                                   tGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGAtt               <  3:921126/85‑1 (MQ=255)
                                                                                                                     taaTTCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGc                                     >  2:571227/2‑61 (MQ=25)
                                                                                                                     taaTTCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGc                                     <  1:571227/60‑1 (MQ=25)
                                                                                                                     aaaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt  <  2:996242/96‑1 (MQ=255)
                                                                                                                     aaaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTGGGCAGACGCTGt  >  1:996242/1‑96 (MQ=255)
                                                                                                                           |                                                                                         
ATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGT  >  NC_000913/3368011‑3368223

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