Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,737,956 2 bp→CG coding (623‑624/1254 nt) avtA → alanine‑alpha‑ketoisovalerate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,737,9560GC100.0% 286.0 / NA 108G208A (GGC→GCG) avtAalanine‑alpha‑ketoisovalerate transaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (57/51);  total (57/51)
*NC_0009133,737,9570CG99.1% 285.9 / NA 108G208A (GGC→GCGavtAalanine‑alpha‑ketoisovalerate transaminase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (56/51);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGGCAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAG  >  NC_000913/3737888‑3738025
                                                                    ||                                                                    
ctcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGACGAGTTGCTGAAGCTTGACGCGCTGGCGAATCaa                                                                >  2:3563621/1‑76 (MQ=255)
ctcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCaa                                                                >  1:425040/1‑76 (MQ=255)
ctcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCaa                                                                >  1:4172929/1‑76 (MQ=255)
 tcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                               >  2:4266951/1‑76 (MQ=255)
 tcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                               <  2:3407938/76‑1 (MQ=255)
 tcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                               <  2:3404354/76‑1 (MQ=255)
 tcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                               <  2:2706793/76‑1 (MQ=255)
 tcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                               <  2:2696997/76‑1 (MQ=255)
  cccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAaca                                                              >  2:4414528/1‑76 (MQ=255)
   ccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAacac                                                             <  1:4533403/76‑1 (MQ=255)
   ccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAacac                                                             >  1:4817555/1‑76 (MQ=255)
   ccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAacac                                                             <  1:1907016/76‑1 (MQ=255)
   ccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAacac                                                             <  2:4050565/76‑1 (MQ=255)
    cGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACg                                                            >  2:2155064/1‑76 (MQ=255)
    cGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACg                                                            <  1:3980403/76‑1 (MQ=255)
      gCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGc                                                          <  2:4812497/76‑1 (MQ=255)
      gCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGc                                                          <  1:4440374/76‑1 (MQ=255)
      gCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGc                                                          <  2:3574822/76‑1 (MQ=255)
       ccGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCa                                                         <  1:3620751/76‑1 (MQ=255)
        cgacgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTTCTGAAGCTTGACGCCCCTGCCAACCCACACCGCAt                                                        >  2:2388758/1‑76 (MQ=255)
        cgacgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGGTGAAGCTTGACGCGCTGGCGAATCAACACGGCAt                                                        >  1:3340659/1‑76 (MQ=255)
        cgacgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCAt                                                        >  2:3683928/1‑76 (MQ=255)
         gacgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCAtt                                                       >  1:1180075/1‑76 (MQ=255)
         gacgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCAtt                                                       >  1:349862/1‑76 (MQ=255)
          acgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATTAACACGGCATTc                                                      >  2:563244/1‑76 (MQ=255)
          acgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTc                                                      <  2:4481266/76‑1 (MQ=255)
          acgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTc                                                      <  1:55066/76‑1 (MQ=255)
           cgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTcc                                                     >  1:2410654/1‑76 (MQ=255)
           cgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTcc                                                     >  1:2456827/1‑76 (MQ=255)
            gaATCCAACAGTGAATGTGATTACTGACGAAGAGTTGCTGCAGCTTGACGCGCTGGCGAATCCACACGGCATTCCg                                                    <  2:2251878/76‑1 (MQ=255)
            gaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCg                                                    >  2:2832678/1‑76 (MQ=255)
             aaTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTTCGc                                                   <  1:3110223/76‑1 (MQ=255)
             aaTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGc                                                   >  2:3641723/1‑76 (MQ=255)
             aaTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCGTGACGCGCTGGCGTATCAACACGGCGTTTCGc                                                   >  1:4843719/1‑76 (MQ=255)
              aTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCt                                                  >  1:4400828/1‑76 (MQ=255)
               accAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTg                                                 <  2:161956/75‑1 (MQ=255)
                 cAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGt                                               <  2:402074/76‑1 (MQ=255)
                  aaCAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtg                                              <  2:4752813/76‑1 (MQ=255)
                  aaCAGGCAAGGAGATTACTGACGAAGGGTTGCTGAAGCTTGACGCGCTCGCGTATCAACACGGCAGTCCGCTgggg                                              >  1:3662392/1‑74 (MQ=255)
                   aCAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtga                                             >  2:1733170/1‑76 (MQ=255)
                     aGGCAATGTGATTTCTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgatt                                           <  2:1503094/76‑1 (MQ=255)
                     aGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgatt                                           >  1:2727449/1‑76 (MQ=255)
                     aGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgatt                                           <  2:474668/76‑1 (MQ=255)
                     aGGCAAGGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgatt                                           <  2:776630/76‑1 (MQ=255)
                       gCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattga                                         >  1:1852234/1‑76 (MQ=255)
                       gCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattga                                         <  2:4521778/76‑1 (MQ=255)
                        cAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattgat                                        >  2:935741/1‑76 (MQ=255)
                        cAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattgat                                        <  1:3852041/76‑1 (MQ=255)
                         aaTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATa                                       >  1:3999752/1‑76 (MQ=255)
                         aaTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATa                                       >  1:3174051/1‑76 (MQ=255)
                          aTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATaa                                      >  1:217296/1‑76 (MQ=255)
                           tgtgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAAc                                     <  1:4335032/76‑1 (MQ=255)
                           tgtgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAAc                                     >  2:1317490/1‑76 (MQ=255)
                            gtgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACg                                    >  1:2225284/1‑76 (MQ=255)
                             tgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGc                                   <  2:2803859/76‑1 (MQ=255)
                              gATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCt                                  >  2:3325662/1‑76 (MQ=255)
                              gATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCt                                  >  1:4480463/1‑76 (MQ=255)
                              gATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGCGATTGATAACGCt                                  >  2:119782/1‑76 (MQ=255)
                               aTTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCtt                                 <  1:1518895/76‑1 (MQ=255)
                               aTTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATCGATAACGCtt                                 <  2:425040/76‑1 (MQ=255)
                                ttACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTa                                <  1:2261566/76‑1 (MQ=255)
                                 tACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTag                               >  2:4666534/1‑75 (MQ=255)
                                  aCTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATg                              >  2:1820091/1‑76 (MQ=255)
                                  aCTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATg                              <  1:2652358/76‑1 (MQ=255)
                                   cTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATgg                             <  2:1602392/76‑1 (MQ=255)
                                     gACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCg                           >  2:4371347/1‑76 (MQ=255)
                                     gACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCg                           >  1:62389/1‑76 (MQ=255)
                                      gcGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGt                          <  2:845655/75‑1 (MQ=255)
                                      aCGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGt                          >  2:1870765/1‑76 (MQ=255)
                                      aCGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGt                          >  2:4795568/1‑76 (MQ=255)
                                      aCGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGt                          >  2:1527495/1‑76 (MQ=255)
                                       cGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTc                         <  1:789550/76‑1 (MQ=255)
                                       cGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTc                         >  2:3758826/1‑76 (MQ=255)
                                        gAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTcc                        <  1:4794512/76‑1 (MQ=255)
                                        gAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTcc                        >  2:4150801/1‑76 (MQ=255)
                                         aaGAGTTGCTGAACCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGAGTATGCTCATGGCATccc                       >  1:315670/1‑76 (MQ=255)
                                          agagTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCg                      <  1:479416/76‑1 (MQ=255)
                                           gagTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGt                     <  1:2388759/76‑1 (MQ=255)
                                           gagTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCCTcccct                     >  2:10497/1‑74 (MQ=255)
                                            agTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGtt                    >  1:3276121/1‑76 (MQ=255)
                                             gTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTc                   <  1:4084754/76‑1 (MQ=255)
                                             gTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTc                   <  2:392764/76‑1 (MQ=255)
                                             gTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTc                   >  1:703734/1‑76 (MQ=255)
                                              ttGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTcc                  >  1:4702414/1‑76 (MQ=255)
                                              ttGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTcc                  >  1:3770764/1‑76 (MQ=255)
                                               tGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTccc                 <  1:3497194/76‑1 (MQ=255)
                                                gCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCg                <  2:3055524/76‑1 (MQ=255)
                                                gCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCg                <  2:1371149/76‑1 (MQ=255)
                                                gCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCg                >  1:545900/1‑76 (MQ=255)
                                                 cTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCgg               <  2:3229760/76‑1 (MQ=255)
                                                    aaGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGCTCCCGGCTa            >  2:673415/1‑76 (MQ=255)
                                                     aGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCCGGTGATTGATAACGCTTATGGCGTCCCGTTCCAGGGTat           >  1:502411/1‑76 (MQ=255)
                                                      gCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatc          <  1:2289120/76‑1 (MQ=255)
                                                       cTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatca         >  2:342698/1‑76 (MQ=255)
                                                        ttGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatcat        >  1:3522415/1‑76 (MQ=255)
                                                        gtGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatcat        <  1:161575/75‑1 (MQ=255)
                                                           gcgcgCTGGCGAATCAACCCGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCtt     <  1:4349804/75‑1 (MQ=255)
                                                            cgcgCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTc    >  2:28110/1‑76 (MQ=255)
                                                            cgcgCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTc    <  1:1840832/76‑1 (MQ=255)
                                                             gcgcTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCa   >  2:2677267/1‑76 (MQ=255)
                                                                    ||                                                                    
CTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGGCAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAG  >  NC_000913/3737888‑3738025

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: