Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,012,261 T→C G267G (GGT→GGC yqeA → putative amino acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,012,2610TC97.4% 139.2 / ‑3.9 39G267G (GGT→GGCyqeAputative amino acid kinase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (15/23);  total (15/24)
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.77e-01

AAATGATTGATGCCGATATGCTGATCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGTTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTG  >  NC_000913/3012116‑3012408
                                                                                                                                                 |                                                                                                                                                   
aaaTGATTGATGCCGATATGCTGATCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTAt                                                                                                                                                <  2:333058/151‑1 (MQ=255)
 aaTGATTGATGCCGATATGCTGATCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTAtg                                                                                                                                               <  1:532567/151‑1 (MQ=255)
        gatGCCGATATGCTGATCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGcc                                                                                                                                          <  2:198479/149‑1 (MQ=255)
           gCCTATATGCTGATCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAg                                                                                                                                     <  1:145511/151‑1 (MQ=255)
              gATATGCTGATCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTgg                                                                                                                                  <  1:157569/151‑1 (MQ=255)
               atatGCTGATCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTgg                                                                                                                                  >  2:305908/1‑150 (MQ=255)
                   gCTGATCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGc                                                                                                                              <  1:142647/150‑1 (MQ=255)
                       aTCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTgcc                                                                                                                          <  2:317623/150‑1 (MQ=255)
                                          aaaaaGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCtctc                                                                                                       >  2:623980/1‑150 (MQ=255)
                                                   ccATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCg                                                                                             >  1:131657/1‑151 (MQ=255)
                                                         aCTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGTTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGc                                                                                       <  2:486303/151‑1 (MQ=255)
                                                          cTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGc                                                                                       >  2:438140/1‑150 (MQ=255)
                                                             cGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACt                                                                                    <  1:122479/150‑1 (MQ=255)
                                                                 aaaGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTa                                                                               <  2:356514/151‑1 (MQ=255)
                                                                     agaATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGt                                                                           <  1:104284/151‑1 (MQ=255)
                                                                     agaATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGt                                                                           <  1:491938/151‑1 (MQ=255)
                                                                      gaATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGt                                                                           >  2:587971/1‑150 (MQ=255)
                                                                        aTGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTAtt                                                                        >  1:198026/1‑151 (MQ=255)
                                                                          gAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTg                                                                       >  2:451683/1‑150 (MQ=255)
                                                                            aCAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGc                                                                    <  1:93258/151‑1 (MQ=255)
                                                                              aGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTCTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCaa                                                                  <  2:350296/151‑1 (MQ=255)
                                                                                     cGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGa                                                             <  2:598734/149‑1 (MQ=255)
                                                                                           ccTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAggg                                                     >  2:19095/1‑151 (MQ=255)
                                                                                           ccTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAggg                                                     >  1:308824/1‑151 (MQ=255)
                                                                                                cATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTg                                                 <  1:145080/150‑1 (MQ=255)
                                                                                                cATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTg                                                 <  1:32625/150‑1 (MQ=255)
                                                                                                             tGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAgg                                    >  1:286985/1‑150 (MQ=255)
                                                                                                             tGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAgg                                    >  1:399510/1‑150 (MQ=255)
                                                                                                              gAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGaa                                  <  1:404890/151‑1 (MQ=255)
                                                                                                                  gCTTTATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGt                              <  2:198021/151‑1 (MQ=255)
                                                                                                                     ttATCGAAGAAGGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGa                            <  2:5127/150‑1 (MQ=255)
                                                                                                                            agaagGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAg                    <  2:353640/151‑1 (MQ=255)
                                                                                                                             gaagGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGt                   <  1:179234/151‑1 (MQ=255)
                                                                                                                              aagGTCATTTTGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTa                  <  1:380880/151‑1 (MQ=255)
                                                                                                                                     ttttGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTaa            >  2:435214/1‑150 (MQ=255)
                                                                                                                                     ttttGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTaa            >  2:11216/1‑150 (MQ=255)
                                                                                                                                       ttGCGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAg          <  2:290627/150‑1 (MQ=255)
                                                                                                                                          cGAAAGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCg      >  1:525223/1‑151 (MQ=255)
                                                                                                                                              aGGCTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCgggtg  >  2:178013/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAATGATTGATGCCGATATGCTGATCATTCTGACTGCGGTAGAAAAAGTGGCCATTAACTTCGGTAAAGAGAATGAACAGTGGCTCGATCGCCTGTCATTAAGTGATGCTGAACGCTTTATCGAAGAAGGTCATTTTGCGAAAGGTTCTATGCTGCCTAAAGTGGAAGCTGCCGCCTCGTTTGCTCGCTCTCGCGCAGGTCGTGAAGCACTGATTACCGTATTGAGCAAAGCGAAAGAAGGGATTGAAGGCAAAACAGGAACGGTGATTTGTCAGTAATTTAAAGGCCGGGTG  >  NC_000913/3012116‑3012408

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