Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,069,802 C→G 19.4% A82G (GCG→GGG)  U712_05265 → Putative quinone oxidoreductase YhfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,069,8020CG19.4% 105.5 / 3.9 53A82G (GCG→GGG) U712_05265Putative quinone oxidoreductase YhfP
Reads supporting (aligned to +/- strand):  ref base C (12/30);  new base G (8/2);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 4.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTGCTTATTCCGGCATTAATTACAAAGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGATGGTGTA  >  CP006881/1069667‑1069940
                                                                                                                                        |                                                                                                                                           
aaGTTGCTTATTCCGGCATTAATTACAAAGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑a                                                                                                                                         <  1:152368/139‑1 (MQ=255)
aaGTTGCTTATTCCGGCATTAATTACAAAGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑a                                                                                                                                         <  4:167589/139‑1 (MQ=255)
  gTTGCTTATTCCGGCATTAATTACAAAGATGGTCTGGCCTGAAAATCTGGAGGATATATCGTCAGAGAGTATCCGATTATTTTTTGCTTTGATGCTGCGGGTTCGGTCGTCACTTCCAATGAT‑CCGCGTTTTGcttagg                                                                                                                                        <  8:183355/139‑6 (MQ=255)
                cATTAATTACAAAGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACTGTCGTGAATTGCAAAGAT‑GCGCGATTTTGGG‑CAGGGGTTGAGGTgg                                                                                                                           >  8:70191/1‑136 (MQ=255)
                     attaataaGATGTTCTGGCCTGTAATGCTGGTGGCAATATGGTCAGAGAGTATCCGCTTTTTTTTGGTTTTGTTGATGCGGGTAAGGTCGTCTATTACAATGAT‑CCGTGTTTTGCGG‑AGGGGGATGAGGTGATCGCGAc                                                                                                                    <  6:7500/133‑1 (MQ=255)
                     attaCAAAGTTGGTCTGGCCTGAAAAGCTGGAGGCAATATCGTCAGAGAGTATCGGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGTATGAGGTTATCGCGAc                                                                                                                    <  6:327929/139‑1 (MQ=255)
                     attaCAAAGATGGTCTGTCCGGAAATTCAGGAGGCAATATCGTCAGAGTGTATCCGCTTATTTTAGGCTTTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGAc                                                                                                                    <  6:180588/139‑1 (MQ=255)
                     attaCAAAGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTTCGG‑TGGGGGATGTGGTGCTCGCGa                                                                                                                     >  7:23315/1‑138 (MQ=255)
                     attaCAAAGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGAc                                                                                                                    <  5:61106/139‑1 (MQ=255)
                     attaCAAAGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGAc                                                                                                                    <  3:244226/139‑1 (MQ=255)
                     attaCAAAGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCGGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGAc                                                                                                                    <  1:151655/139‑1 (MQ=255)
                      ttaCAACGATGGGCTGGCCGGAAAAGCAGGAGGCAATAGCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGGCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGGGAGCGCGACa                                                                                                                   <  7:101769/139‑1 (MQ=255)
                          aaaGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTGTTACAATTAT‑CCGCCTTTTGCTG‑AGGGGGATGAGGTGTTCGCGACAAga                                                                                                                >  7:286892/1‑137 (MQ=255)
                            aGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTAt                                                                                                             <  1:27301/139‑1 (MQ=255)
                              aTGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAAAGAT‑AAGCCTGTTTCGGTGGGGGGGTGTGG‑GATCTCGACAAGCAAt                                                                                                             >  5:204016/1‑137 (MQ=255)
                                ggTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTGGCTTACAATGAT‑AACCGGTTTGGGG‑TGGGGGATGTGGGGATAGCGCCAAGCAATTAGc                                                                                                         >  2:118926/1‑139 (MQ=255)
                                ggTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATTAT‑CCGGGGTTTTGGG‑AGGGGGGTGAGGTGTTCGCGACAACCTATGTGc                                                                                                         >  7:118459/1‑139 (MQ=255)
                                  tCTGGCCGGAAAAGCAGGTGTCAATATCGTCAGTGAGTATCCGCTTATTTAAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCAAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTc                                                                                                       <  8:120105/139‑1 (MQ=255)
                                   cTGGCCGGTAAAGCCTGAGGAAATATCGTCAGTGAGTATACGATTTTTTTAGGCATTGATGCGGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCg                                                                                                      <  6:204016/139‑1 (MQ=255)
                                   cTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTAGCCGCTTATTTTAGGCAGTGATGCGGCGGGTACGGGCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCg                                                                                                      <  6:6887/139‑1 (MQ=255)
                                      gCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTGCAATGAT‑CCGCGGTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGgtg                                                                                                   >  7:23494/1‑139 (MQ=255)
                                                     ggCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑ACGCGGTTTGCGG‑AGGGGGATGTGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCg                                                                                    >  6:67615/1‑139 (MQ=255)
                                                         atatCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGAc                                                                                                                    <  8:335136/103‑1 (MQ=255)
                                                         atatCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGAc                                                                                                                    >  7:335136/1‑103 (MQ=255)
                                                               tCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGGGG‑GGGGGGGTGAGGTGATCGCGACAAGCTATGCGCTCGGTGTCTCACGTGATGGCGGATTAAGTg                                                                           >  5:129204/1‑138 (MQ=255)
                                                               tCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCGCGTCCAATGAT‑CCGCGGTTTGGGG‑TGGGGGGTGAGGTGATAGGGCCAAGCGATTAGCTCGGTGGCTCACGTGATGGCGGATTTAGTg                                                                           >  2:251224/1‑138 (MQ=255)
                                                                     agTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCGGTTGCAATGAA‑ACGCCTTTTTCGG‑AGGGGGATGAGGGGCTCGCGACAACCAATGAGCTCGGTGTCGAACGGGATGGCGGATTAATTGaaaaac                                                                     >  8:169323/1‑134 (MQ=255)
                                                                       tataCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCGCTTACAATGAT‑ACGCGGTTTTTGG‑AGGGGGAAGAGGGGATAGCGCCAAGCAATTTGCTAGCTGTCTCACGTGATGGGGGATTAATTGAATAAGcat                                                                  >  4:178719/1‑137 (MQ=255)
                                                                         tCCGCTTTTTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTc                                                                 <  1:251224/138‑1 (MQ=255)
                                                                         tCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCg                                                                <  1:118926/139‑1 (MQ=255)
                                                                             cTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTGTCTTACAATGATACCGC‑ATTTGGGG‑ATGGGGATGTGGTGGTCGCGCCAAGCAATTAGATCGGGGTCTCACGTGATGGCGGCTTTATTGAAAACGCATCgggga                                                            >  4:93759/1‑136 (MQ=255)
                                                                                aTTTTTGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCATTGAT‑CGGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGcc                                                           <  7:169323/137‑1 (MQ=255)
                                                                                aTTTTAGGCTTTGATGCTGCGGGTACTGTCGTCTCTTCCAATGAT‑CCGCGGTTTGGGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACCTGATGGCGGATTAAGTGTATTCGCTTCGGTGcc                                                           <  7:70191/137‑1 (MQ=255)
                                                                                 ttttAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTgg                                                        >  3:135072/1‑139 (MQ=255)
                                                                                 ttttAGGCATTGATGCTGCGGGTACGGTCGTCGGGTGCAATGAT‑ACGCGGGTTGGGG‑GGGGGGGTGAGGGGGTCGCGGCCAGCGATGTGCGCGGGGGGTCGCGGGAGTGCGGGTTTAGGGGATAGGCGTCGGTGCCggg                                                        >  6:262594/1‑139 (MQ=255)
                                                                                     aGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGa                                                     <  4:135072/138‑1 (MQ=255)
                                                                                           tGATGCTGCGGGTACGGTCGTGGCTTCCAATGAT‑ACGCGGTTTGGGG‑AGGGGGAGGAGGTGTTGGCGGCAAGCGATGAGCTCGGTGTCTGACGGGATGGCGGGTTTATTGAATACGCTTCGGTGGCTGGTGgctggctg                                              >  8:338653/1‑139 (MQ=255)
                                                                                                                      cAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGtt                    >  5:335041/1‑138 (MQ=255)
                                                                                                                        aTGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATgg                                                                                      <  1:77081/70‑1 (MQ=255)
                                                                                                                        aTGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATgg                                                                                      >  2:77081/1‑70 (MQ=255)
                                                                                                                              ccGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAaagaag            >  2:67975/1‑139 (MQ=255)
                                                                                                                               cggtttttGGGG‑TGGGTGTTGTGGTGTTCGCGACAAGTTATTTGTTCGGTGTGTCAGGTGTTGGCGGTTTTTGTGATTAGGTTTTGTTGCTTGGTTATTGGTTGGTGCTTTTGCCACAGAATCTTTCGTTAAaagtag            <  8:269140/134‑1 (MQ=255)
                                                                                                                               cgcgTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAaagaag            <  8:286892/138‑1 (MQ=255)
                                                                                                                               cgcgTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAaagaag            <  3:93759/138‑1 (MQ=255)
                                                                                                                                 cggtttGCGG‑AGGGGGATGAGGTGAGCGCGACAAGCTATGAGCTCGGGGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGGGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGa         <  2:270864/136‑1 (MQ=255)
                                                                                                                                 cgTTTTGCGT‑AGGGGGATGAGGGGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGa         <  1:300058/139‑1 (MQ=255)
                                                                                                                                 cgTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGa                                                                                  >  7:75386/1‑66 (MQ=255)
                                                                                                                                 cgTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGa                                                                                  <  8:75386/66‑1 (MQ=255)
                                                                                                                                 cgTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCTGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGa         <  3:178719/139‑1 (MQ=255)
                                                                                                                                 cgTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGa         <  4:247694/139‑1 (MQ=255)
                                                                                                                                 cgTTTGGCGG‑AGGGGGAGGAGGCGATCGCGACAAGCTATGAGCTCGGGGTCTCACGTGATGGCGGATTAAGGGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGa         <  2:152350/139‑1 (MQ=255)
                                                                                                                                 cgTGTTGCGG‑AGGGGGATGAGGGGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGa         <  8:342961/139‑1 (MQ=255)
                                                                                                                                   ttttGCGG‑TGGGGGTTGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGATg       <  1:304846/139‑1 (MQ=255)
                                                                                                                                   tgttGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATTGCGGATAAAGTGAATACGCTTCGGTGCCTGGGGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGATg       <  4:307001/137‑1 (MQ=255)
                                                                                                                                        cGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGATGGTGTa  >  7:269140/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                           
AAGTTGCTTATTCCGGCATTAATTACAAAGATGGTCTGGCCGGAAAAGCAGGAGGCAATATCGTCAGAGAGTATCCGCTTATTTTAGGCATTGATGCTGCGGGTACGGTCGTCTCTTCCAATGAT‑CCGCGTTTTGCGG‑AGGGGGATGAGGTGATCGCGACAAGCTATGAGCTCGGTGTCTCACGTGATGGCGGATTAAGTGAATACGCTTCGGTGCCTGGTGACTGGCTGGTGCCTTTGCCACAGAATCTTTCGTTAAAAGAAGCGATGGTGTA  >  CP006881/1069667‑1069940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑