Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,356,738 G→A 12.2% intergenic (+27/‑166) glpC → / → rpnE anaerobic glycerol‑3‑phosphate dehydrogenase subunit C/inactive recombination‑promoting nuclease‑like protein RpnE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,356,7380GA12.2% 71.4 / 2.2 31intergenic (+27/‑166)glpC/rpnEanaerobic glycerol‑3‑phosphate dehydrogenase subunit C/inactive recombination‑promoting nuclease‑like protein RpnE
Reads supporting (aligned to +/- strand):  ref base G (13/14);  new base A (2/2);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGT  >  NC_000913/2356607‑2356875
                                                                                                                                   |                                                                                                                                         
tgtccggtggtCACCGTCTGCGACACCGTGAAGTGCCAGACTGAGATGGCCACATGTCTTCGCTGCTATCTTCCGATTTCGCTACTGGCCCAGGCGTTGGCTTAAACTCCTTTCTGATGCCCGGTACCCACATGCTTACCggg                                                                                                                                <  2:534374/138‑1 (MQ=255)
gATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCggg                                                                                                                                <  1:436846/143‑1 (MQ=255)
    tggtggTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGGTACCGGGCAtt                                                                                                                            >  2:173916/1‑143 (MQ=255)
      gtggtCGCCGGCTGCGAGGCCGGTGAGTGGCAGATTGGGATGTCGACAGGTCTTCGCTGCGAACATCCGATGACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAtttt                                                                                                                          <  2:225204/143‑1 (MQ=255)
       tggtCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACACGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAttttt                                                                                                                         <  1:254661/143‑1 (MQ=255)
       tggtCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAttttt                                                                                                                         <  1:572279/143‑1 (MQ=255)
        ggtCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTg                                                                                                                        <  2:391677/143‑1 (MQ=255)
         gtCAGCGACCGCGAACCCTGTGAATGGCAGAGTGAGAGGTCCACAAGTCTGCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTAAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGc                                                                                                                       <  1:93506/143‑1 (MQ=255)
         gtCACCGACTGCGAAACCTGTAAATGGCAGAGTGAGATGTCCACAAGTCTTCGCTGCGCACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTCTGc                                                                                                                       >  1:664917/1‑143 (MQ=255)
          tCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCg                                                                                                                      >  2:462758/1‑143 (MQ=255)
               gACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTacac                                                                                                                 <  2:621843/143‑1 (MQ=255)
               gACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTacac                                                                                                                 >  2:198602/1‑143 (MQ=255)
                           tGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATATGGTTACCGGGCATGTTTGCGGACAAGAGTCCATGccc                                                                                                     >  1:392924/1‑143 (MQ=255)
                                                   acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCCTTTTTGCGTACACGCTTCCGTGCCCAATGTATGCGTTGCAACGCAGTGa                                                                             >  1:580412/1‑143 (MQ=255)
                                                   acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGa                                                                             >  1:51255/1‑143 (MQ=255)
                                                        tCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTACGCATGTGGTTACCGGGCATTTTTGCGTACACGCTTCCGTGCCCAATGTTTCCATTGCAACGCAGTGAACAtt                                                                        >  2:467369/1‑143 (MQ=255)
                                                         cTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCGGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGGAAAGGAATCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTc                                                                       <  1:128954/143‑1 (MQ=255)
                                                           tCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAACATTCct                                                                     >  1:629691/1‑143 (MQ=255)
                                                            cGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCctc                                                                    >  1:318287/1‑143 (MQ=255)
                                                                          gATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTCCACGCTTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGCAAACGTCTCGc                                                      >  2:295280/1‑143 (MQ=255)
                                                                           aTTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCa                                                     >  1:371924/1‑143 (MQ=255)
                                                                                     tGGCCCAGGGGCTGTCGTAACCGCTGTTCTGAGGCCGGGTAATCATGTGGTTACCGGGCATTTTTGCGTCCACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGaaa                                           <  2:207516/143‑1 (MQ=255)
                                                                                      ggCCCAGGCGCTGGCGTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAAc                                          <  1:98469/143‑1 (MQ=255)
                                                                                          cAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGc                                      <  1:198602/143‑1 (MQ=255)
                                                                                             gcgtTGGCTTAAACTCCTTTCTTATGCCCGGTAAGTATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGa                                   <  1:83753/143‑1 (MQ=255)
                                                                                                      tAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCACAgg                          >  1:82931/1‑143 (MQ=255)
                                                                                                          ctctTTTCTGATGCCCGGTAAGCATGTGGTTCCCGGGCATTTTTACGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAAGTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTg                      <  2:28637/143‑1 (MQ=255)
                                                                                                           tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                     >  2:646722/1‑143 (MQ=255)
                                                                                                              tttCTGGTGCCCGGTAGGTATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCGATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGAt                  <  2:363002/143‑1 (MQ=255)
                                                                                                                      gCCCGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTCTGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGCTGTCAAAGGCCCGGGAACACCGtca          >  2:421429/1‑141 (MQ=255)
                                                                                                                              cgCATGTGGTTACCGGGCAGTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAAGTCCTCTGAAAACGGCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGt  <  2:207466/142‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                         
GATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGT  >  NC_000913/2356607‑2356875

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