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,9560GC99.2% 346.8 / ‑5.4 131G208A (GGC→GCG) avtAalanine‑alpha‑ketoisovalerate transaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (66/64);  minor base T (1/0);  total (67/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01
*NC_0009133,737,9570CG99.2% 351.6 / ‑4.8 131G208A (GGC→GCGavtAalanine‑alpha‑ketoisovalerate transaminase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (66/64);  total (67/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 131 total aligned reads displayed.

GTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGGCAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCA  >  NC_000913/3737886‑3738024
                                                                      ||                                                                   
  ctcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCaa                                                               <  1:210598/76‑1 (MQ=255)
  ctcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCaa                                                               >  2:3739364/1‑76 (MQ=255)
   tcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                              >  1:4097049/1‑76 (MQ=255)
   tcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGAGCTGGCGAATCaaa                                                              >  1:2880339/1‑75 (MQ=255)
   ccccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                              <  1:3324685/75‑1 (MQ=255)
    cccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCaaac                                                             >  2:5385068/1‑74 (MQ=255)
    cccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAaca                                                             <  1:918672/76‑1 (MQ=255)
    cccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAaca                                                             >  1:3705204/1‑76 (MQ=255)
     ccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAacac                                                            <  1:5630784/76‑1 (MQ=255)
     ccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGATGAAGCTTGACGCGCTGGCGAATCAacac                                                            >  2:2995605/1‑76 (MQ=255)
      cGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGGTGCTGAAGCTTGACGCGCTGGCGAATCAACACg                                                           >  1:5150302/1‑76 (MQ=255)
       ggCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACgg                                                          >  2:1855352/1‑76 (MQ=255)
       ggCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACgg                                                          >  1:1672184/1‑76 (MQ=255)
       ggCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACgg                                                          <  1:2247932/76‑1 (MQ=255)
       ggCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACgg                                                          <  2:4683454/76‑1 (MQ=255)
        gCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGc                                                         <  2:2058064/76‑1 (MQ=255)
         ccGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCa                                                        >  1:5273091/1‑76 (MQ=255)
          cgccgaATGCAACAGGCCATCTTATTACTGACGAAGCGTTGCTGAATCTTGTCGCGCCGGCGCATCAACACGGCAt                                                       <  1:5252873/73‑1 (MQ=255)
            acgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTc                                                     >  2:1955825/1‑76 (MQ=255)
            acgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTc                                                     >  2:3326839/1‑76 (MQ=255)
            acgaATCCAAAAGGCAATGCGATTACTGACGAAGAGATGCTGAAGCATGACGCGCTCGCGAATCAACACGGCAATc                                                     >  1:1911764/1‑76 (MQ=255)
             cgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATAAACACCGCATccc                                                    >  2:5733272/1‑76 (MQ=255)
              gaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCg                                                   >  1:5461096/1‑76 (MQ=255)
              gaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCg                                                   >  2:4344783/1‑76 (MQ=255)
              gaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCg                                                   >  2:5812822/1‑76 (MQ=255)
               aacccaccaGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTTACTCGCTGGCGGATCAACACGGCATTCCGc                                                  <  1:4305448/73‑1 (MQ=255)
               aaTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGATGCTTGAGGCGGTGGCGACTCACCACGGCATTCGGc                                                  <  2:4003174/76‑1 (MQ=255)
               aaTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCCTTCCGc                                                  >  2:3977670/1‑76 (MQ=255)
               aaTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCAt                                                       <  2:424392/71‑1 (MQ=255)
               aaTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGc                                                  >  2:4932763/1‑76 (MQ=255)
                aTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCt                                                 <  1:2737409/76‑1 (MQ=255)
                aTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCt                                                 >  2:2085689/1‑76 (MQ=255)
                 tCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTg                                                <  1:640820/76‑1 (MQ=255)
                 tCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTg                                                <  1:231046/76‑1 (MQ=255)
                   cAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGt                                              <  1:1790268/76‑1 (MQ=255)
                   cAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGt                                              <  2:1488446/76‑1 (MQ=255)
                    aaCAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtg                                             >  2:3251987/1‑76 (MQ=255)
                    aaCAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtg                                             >  2:3247134/1‑76 (MQ=255)
                     aCAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtga                                            <  1:5529878/76‑1 (MQ=255)
                       aGGCAATGGGATTACTGACGAAGAGTGGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgatt                                          <  2:263738/76‑1 (MQ=255)
                        ggCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattg                                         <  2:508383/76‑1 (MQ=255)
                         gCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattga                                        <  1:3748329/76‑1 (MQ=255)
                         gCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattga                                        >  2:5730382/1‑76 (MQ=255)
                          cAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattgat                                       <  1:4207692/76‑1 (MQ=255)
                           aaTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATa                                      >  2:3256989/1‑76 (MQ=255)
                           aaTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATa                                      >  2:3697414/1‑76 (MQ=255)
                           aaTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATa                                      >  1:977817/1‑76 (MQ=255)
                            aTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATaa                                     >  1:4030296/1‑76 (MQ=255)
                            aTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATaa                                     >  1:3219483/1‑76 (MQ=255)
                            aTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGCGGTTCATaa                                     >  2:5580591/1‑76 (MQ=255)
                             tgtgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGGGATTGATaaa                                    >  1:627291/1‑75 (MQ=255)
                              gtgATTGCTGACGAAGAGTTGCTGAAGCTTTACGCGCTGGCGAATGAGCACGGCATTCCGCTGGTGATTGATAACg                                   <  2:1378916/76‑1 (MQ=255)
                               tgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGc                                  <  2:2510926/76‑1 (MQ=255)
                                 aTTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCtt                                >  1:1312174/1‑76 (MQ=255)
                                 aTTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCtt                                >  2:1356593/1‑76 (MQ=255)
                                 aTTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCtt                                <  1:3329280/76‑1 (MQ=255)
                                 aTTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCtt                                <  1:1908966/76‑1 (MQ=255)
                                  ttACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCCTa                               >  2:452499/1‑76 (MQ=255)
                                   tACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTAt                              <  1:3392200/76‑1 (MQ=255)
                                   tACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTAt                              <  1:3326839/76‑1 (MQ=255)
                                   tACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTAt                              <  1:4885132/76‑1 (MQ=255)
                                    tcTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATg                             <  2:4633867/75‑1 (MQ=255)
                                    aCTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTCTg                             >  2:5843907/1‑76 (MQ=255)
                                     cTGACGAAGAGTTGCTCGAGCTTGACGAGCTTGTCAAGCGACACGGCATTCCGCTGGGGATTTATAACGCTTAttg                            >  1:1662048/1‑74 (MQ=255)
                                      tGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGc                           <  2:1622967/76‑1 (MQ=255)
                                      tGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGc                           <  1:3543214/76‑1 (MQ=255)
                                       ggcGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCg                          <  2:5004043/74‑1 (MQ=255)
                                        aCGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGt                         >  1:370603/1‑76 (MQ=255)
                                          gAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTcc                       <  1:2467028/76‑1 (MQ=255)
                                          gAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTcc                       >  2:4551717/1‑76 (MQ=255)
                                          gAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTcc                       >  1:1455172/1‑76 (MQ=255)
                                          gAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTcc                       >  2:620849/1‑76 (MQ=255)
                                          gAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTcc                       >  2:1929151/1‑76 (MQ=255)
                                            agagTTGCTTAAGCTTGCCGCGCTGGCGACTCGACACGGCATTCCGCTGGTGCTTGATAACGCTTATGGCGTCCCg                     <  1:2010975/76‑1 (MQ=255)
                                            agagTTGCTGAAGCGTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTccca                     <  2:3191183/76‑2 (MQ=255)
                                             gagTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGt                    >  2:5056179/1‑76 (MQ=255)
                                             gagTTGCTGAACCTTGACGCGCTGTCGACTCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGt                    <  1:3513154/76‑1 (MQ=255)
                                              agTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGtt                   <  1:3626146/76‑1 (MQ=255)
                                              agTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGtt                   >  2:1191065/1‑76 (MQ=255)
                                               gTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTc                  <  1:1245754/76‑1 (MQ=255)
                                                  gCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCg               <  1:3585519/76‑1 (MQ=255)
                                                   cTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCgg              <  2:1312174/76‑1 (MQ=255)
                                                   cTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCgg              <  2:4241114/76‑1 (MQ=255)
                                                    tGAAGCTTGACGCGCTGGCGAATCAACACGGCCTTCCGCTGGTGAGTGATAACTCATATGGCGTCCCGTTACCggg             >  1:3233422/1‑76 (MQ=255)
                                                    tGAAGCTTGACGAGCTGGCGAATCAACACGGCATTCCGCTGGCGAGTGAGAACGCTTATGGCGACACGTTCCCcgg             >  1:279579/1‑76 (MQ=255)
                                                     gAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGt            <  2:2414395/76‑1 (MQ=255)
                                                     gAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGt            <  1:86810/76‑1 (MQ=255)
                                                      aaGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTa           <  2:4097048/76‑1 (MQ=255)
                                                       aGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTat          >  1:4613989/1‑76 (MQ=255)
                                                         cTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatca        <  2:2285789/76‑1 (MQ=255)
                                                          ttGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatcat       >  1:2762781/1‑76 (MQ=255)
                                                          ttGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatcat       >  2:2488018/1‑76 (MQ=255)
                                                          ttGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatcat       <  1:5642364/76‑1 (MQ=255)
                                                          ctGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatcat       <  2:1275734/75‑1 (MQ=255)
                                                            gACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCt     >  1:2677149/1‑76 (MQ=255)
                                                            gACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCt     <  2:4833801/76‑1 (MQ=255)
                                                             aCGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCtt    >  1:1566481/1‑76 (MQ=255)
                                                             aCGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCtt    <  1:2633240/76‑1 (MQ=255)
                                                              cgcgcgGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTc   <  2:101569/76‑1 (MQ=255)
                                                               gcgcTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCa  >  2:3627620/1‑76 (MQ=255)
                                                                      ||                                                                   
GTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGGCAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCA  >  NC_000913/3737886‑3738024

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