Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,098,045 G→C 9.8% G54R (GGG→CGG)  ghrA → glyoxylate/hydroxypyruvate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,098,0450GC 9.8% 120.1 / 5.1 49G54R (GGG→CGG) ghrAglyoxylate/hydroxypyruvate reductase A
Reads supporting (aligned to +/- strand):  ref base G (27/17);  new base C (3/2);  total (30/19)
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.37e-01

ATACCCAATGGTGGATTGAGGCACTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGG  >  NC_000913/1097917‑1098186
                                                                                                                                |                                                                                                                                             
aTACCCAATGGTGGATTGAGGCACTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGTCATCCTCCTGTTGAAATTCTGGCAGGGCGCGATATaaaa                                                                                                                                 >  1:160475/1‑143 (MQ=255)
    ccAATGGTGGATTGAGGCACTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCAGGCAGGGCGCGATGTTAAGGCgg                                                                                                                             >  2:554834/1‑143 (MQ=255)
    ccAATGGTGGATTGAGGCACTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCAGTCATTCACAGGTTGAAATGCTAGCAGGGCGCGATCTTAAAGcag                                                                                                                             >  2:747849/1‑141 (MQ=255)
      aaTGGTGGATTGAGGCACTGCGCAAAGCTATTCCTCAGGCAAGCGTCAGCGCATGGAAACGCGGAGATAATGACTCTGATGATTATACTTTATTCTGGCATCCTCATGTTGAAATGCTGGCAGGGCGGGATCTTAAAGCGgtg                                                                                                                           >  1:371659/1‑143 (MQ=255)
        tggtggATTGAGGCACTGCGCAAAGCTATTCCGCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGtt                                                                                                                         <  2:33949/143‑1 (MQ=255)
         ggtggATTGAGGCACTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTc                                                                                                                        <  1:372343/143‑1 (MQ=255)
               ttGAGGCACTGCGCAAAGCTATTCCTCAGGCAAGAGACAGAGCAGGGAAAAGCGGAGATAATGAATCTGAGGACTATGCGATAGTCTGGCATCCGCCTGTGGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTCCGCACta                                                                                                                  >  2:126939/1‑142 (MQ=255)
                tGAGGCACAGCGGAGAGGTAGTCGTCAGGCAAGAGGCAGAGCAGGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCg                                                                                                                 <  1:430758/143‑1 (MQ=255)
                  aGGCACTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTCTGCTTTAGTCTGGCATAATCATGTTAAAATGCTGGCAGAGCGCGATCTTAAAGCGGTGTTCGCACTCggg                                                                                                               >  1:18578/1‑143 (MQ=255)
                  aGGCACTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCggg                                                                                                               >  1:414950/1‑143 (MQ=255)
                     cACTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCCGGCATCCTCCTGTTGAAATGCTGCCACGTCGCGATCTCAAAGCGGCGCTCGCACTCGGGGcc                                                                                                            >  2:744375/1‑143 (MQ=255)
                      aCTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATACTGGCAGGGAGCAATATTAAAGCGGTGTTCGCACTGGGGGcag                                                                                                           >  1:200270/1‑141 (MQ=255)
                       cTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTCTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCATCAAAAGATAATCTCACTCGGGGACgg                                                                                                          >  2:570548/1‑143 (MQ=255)
                          cgcAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCACGGCGCGATTTTAAAGCGGTGTTCGCACTCGGGGCCGgtgt                                                                                                       >  1:547099/1‑143 (MQ=255)
                          cgcAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCACGGCGCGATTTTAAAGCGGTGTTCGCACTCGGGGCCGgtgt                                                                                                       >  1:547098/1‑143 (MQ=255)
                             aaaGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGa                                                                                                    <  2:720401/143‑1 (MQ=255)
                                gCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCCCCCGGGGCCGGCGATGACTc                                                                                                 >  2:256332/1‑143 (MQ=255)
                                  tATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCg                                                                                                                       <  1:321887/119‑1 (MQ=255)
                                  tATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCg                                                                                                                       >  2:321887/1‑119 (MQ=255)
                                     tcttCAGGCAAGATTCAGAGCATGTAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGGTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTAttt                                                                                            <  1:543537/140‑1 (MQ=255)
                                      ccgcAGGCCAGAGTCAGAGCCTGGAGGAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCAGCCTCCTGTTTAAATGCTGGCAGGGCGCGATCTTTAAGCGGTGTTCGCTCTCGGGGCCGGTGTTGATTCTAtttt                                                                                           <  2:488531/140‑1 (MQ=255)
                                                   tgagagCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAg                                                                              <  2:176546/141‑1 (MQ=255)
                                                     agagCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAACGCGGTATTCGCAACCGGGGGAGGTGTTACTACTATATAGAGCAAGCTACACGc                                                                            >  2:312017/1‑143 (MQ=255)
                                                           tGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACAccc                                                                      >  1:715404/1‑143 (MQ=255)
                                                           tGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACAccc                                                                      >  2:580395/1‑143 (MQ=255)
                                                                 aacGGAGATAATGATTCTGCTGATTATGCTTTAGTCTGGCATTCTCGTGTTGAAATGCTGGTAGGGCGCGATATTAAAGCGGTGTTCGCACTAGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAAt                                                                <  1:357026/141‑1 (MQ=255)
                                                                  gCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGCTTCTATTTTGAGCAAGCTACAGGCACACCCTGaacgg                                                               >  2:549522/1‑140 (MQ=255)
                                                                   aagagaTAATGACTATGCTGATTAGGCTGGAGCCTGGCATCCTCCAGTTGAAATGCCGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGc                                                              <  1:31210/141‑1 (MQ=255)
                                                                                  tgctgAGGATACTTGAGTGTGGGATCCTCCAGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTGTGCAGAGGGGGCCGGGGGGGATTCGATATTGAGCAAGCTACAGAGACACGCTGAAATGGTGAACGCCTCAGtta                                               >  2:366328/1‑142 (MQ=255)
                                                                                    ctgATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGCTTCTATTTTGAGCAAGCTACAGGCACACACTGAAATGCTGAACCCTTCTGTTcaa                                             >  1:8925/1‑141 (MQ=255)
                                                                                    ctgATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACCCAGAGGTGGTTTAAAGTCTGTGTTGAGCAAGCTACAGGCACACACTGAAATACTGAACACTACTGttaaa                                             >  1:172845/1‑140 (MQ=255)
                                                                                       aTTATGCTTTAGTCTGGCATCGTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGCGGGCGGTGTTGATTCTATTTTGAGCAAGCTACAGTCATACTCTGAAATGCTGAACGCTTCGGTTTCACtt                                          >  2:643914/1‑143 (MQ=255)
                                                                                             cgttAGTCTGGCCTCCTCCTGTTTAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCATTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGc                                    <  2:222196/141‑1 (MQ=255)
                                                                                                  gTCTGTCATCCTCGTGTTGATATGCTGGCCGGGCTCGATTTTAAAGCGGTGTTTGCACTCGGGGCCGGTGTTGATTCTATTGTGAGCATGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGa                               <  2:143175/143‑1 (MQ=255)
                                                                                                   tCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGaa                              >  1:295269/1‑143 (MQ=255)
                                                                                                    cTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTggccg                             >  1:770089/1‑140 (MQ=255)
                                                                                                    cTGGCATCCTCCTGTTGAAATGCTGGCAGCGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGTGCAAGCTACAGCCACACACTGAAACGCTAAGACCTTCTGTTCCACTATTACGCCAGAAAg                             >  2:153408/1‑143 (MQ=255)
                                                                                                       gCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATa                          <  1:580395/143‑1 (MQ=255)
                                                                                                       gCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATa                          <  2:715404/143‑1 (MQ=255)
                                                                                                            ctcctGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGACTGGAAGATACCGGt                     >  1:73567/1‑143 (MQ=255)
                                                                                                             tcctGTTGAGATGCTGGGCGGGCGGGATGCTAGAGCGGTGGTCGCACGGGGGGGCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTCCCACTTTTGCGCCTGGAAGATACCGGTa                    <  1:652597/143‑1 (MQ=255)
                                                                                                             tcctGTTGAAATGCTGGCAGAGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGCGTTGAATCTCTATTGAGCAAGCTACAGGCACACCCTGAACAGCTCCAGCCCTTTTTCCGCCTTTTTCGCCTGGAAGATACCgcga                    >  2:538766/1‑140 (MQ=255)
                                                                                                              cctGTTGAAATGCTGGCACGGCGCGATTTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTAt                   <  2:547098/143‑1 (MQ=255)
                                                                                                              cctGTTGAAATGCTGGCACGGCGCGATTTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTAt                   <  2:547099/143‑1 (MQ=255)
                                                                                                                     aaaTGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAGCCCTTCTGTTCTACTTTTTCGCCTGGTTTCTACCGCTATGGGCggt            >  1:35278/1‑141 (MQ=255)
                                                                                                                       aTGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGCAGAACCCGATATGGGCGAGcc          >  2:88087/1‑142 (MQ=255)
                                                                                                                            ggCAGGGCGCGATCTTAAAGCGGGGTTCGCACTCGGGGCCGGGGTTGAGTCTATTTTCAGCAAGCTACAGGTGCATTCTGACATGCTGAACCATTCTGCTCCACTTTTTCGCCTGGACAATACCGGTATAGGCGAACACATGc     >  2:10204/1‑143 (MQ=255)
                                                                                                                             gCAGGGCGCGATCTTAAAGCGGTTTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCGGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCa    <  1:261247/143‑1 (MQ=255)
                                                                                                                               aGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAgg  <  1:693396/143‑1 (MQ=255)
                                                                                                                                |                                                                                                                                             
ATACCCAATGGTGGATTGAGGCACTGCGCAAAGCTATTCCTCAGGCAAGAGTCAGAGCATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGG  >  NC_000913/1097917‑1098186

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