Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,122,270 (A)6→5 coding (9/924 nt) yceM → Gfo/Idh/MocA family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,122,2650A.96.5% 111.6 / ‑1.2 29coding (4/924 nt)yceMGfo/Idh/MocA family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (10/18);  minor base G (1/0);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGAC  >  NZ_CP009273/1122179‑1122350
                                                                                      |                                                                                     
ggCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaa                                                                                    >  1:1219347/1‑90 (MQ=255)
  cccGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaat                                                                                  >  2:760530/1‑89 (MQ=255)
   ccGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaatt                                                                                 <  2:557774/90‑3 (MQ=255)
     gcgTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaattac                                                                               >  2:1028886/1‑86 (MQ=255)
     gcgTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaattac                                                                               >  2:914872/1‑86 (MQ=255)
     gcgTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaattac                                                                               >  2:1410695/1‑86 (MQ=255)
       gTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGt                                                                            <  1:1230938/90‑1 (MQ=255)
             gCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGc                                                                      <  2:1219347/90‑1 (MQ=255)
                cTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTa                                                                   <  1:152580/90‑1 (MQ=255)
                       tCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGAt                                                            <  1:168235/90‑1 (MQ=255)
                       tCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGAt                                                            <  2:724320/90‑1 (MQ=255)
                       tCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGAt                                                            <  2:933975/90‑1 (MQ=255)
                                       gAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCa                                            >  2:1254613/1‑90 (MQ=255)
                                         aaCAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCaaa                                          >  1:145868/1‑90 (MQ=255)
                                             gCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGc                                      >  2:1415479/1‑90 (MQ=255)
                                             gCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGc                                      >  1:892268/1‑90 (MQ=255)
                                              cGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCg                                     <  1:104683/90‑1 (MQ=255)
                                                  ttGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGt                                 >  2:952460/1‑90 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  1:182449/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  1:1028886/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  1:1410695/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  2:1471630/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  2:145868/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  1:861434/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTCAcc                             <  1:341321/90‑1 (MQ=255)
                                                           tGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGc                        >  1:937825/1‑90 (MQ=255)
                                                                    aTCTGGGAGATTAAAGGGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGGGGTTACCGGTGCTGGCGGCAg               >  1:1393920/1‑90 (MQ=255)
                                                                            gaTTAAAGTGCAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGAc       <  1:1518565/90‑1 (MQ=255)
                                                                              ttAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTg     >  2:672924/1‑90 (MQ=255)
                                                                              ttAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTg     >  2:1234223/1‑90 (MQ=255)
                                                                              ttAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTg     >  2:886172/1‑90 (MQ=255)
                                                                                aaaGTGAAAAATTACGTATTGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGa   <  1:914872/90‑1 (MQ=255)
                                                                                aaaGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGa   <  1:952460/90‑1 (MQ=255)
                                                                                aaaGTGAAAAATTACGTACCGGCGTAGTGGGATTAGGTGGCATTGCGCTAAAAGCGGGGTTATTGGTGCTGGCGGCAGCGTCTGACGGGa   <  1:1345517/90‑1 (MQ=255)
                                                                                 aaGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGAc  >  2:648208/1‑90 (MQ=255)
                                                                                      |                                                                                     
GGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGAC  >  NZ_CP009273/1122179‑1122350

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