Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,557 2 bp→GC coding (451‑452/930 nt) ilvE → branched‑chain amino‑acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,5570CG99.3% 360.2 / NA 138R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (68/69);  total (68/70)
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
*NC_0009133,950,5580GC100.0% 368.7 / NA 139R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (68/71);  total (68/71)

Only 100 of 139 total aligned reads displayed.

GAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCG  >  NC_000913/3950489‑3950626
                                                                    ||                                                                    
 aGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAac                                                               >  1:2875141/1‑76 (MQ=255)
  gCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAaca                                                              >  1:1070187/1‑76 (MQ=255)
  gCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAaca                                                              >  1:4186154/1‑76 (MQ=255)
   cGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAacac                                                             <  2:2622121/76‑1 (MQ=255)
    gTATCTGGGCGCAGCAGCGATGGAGCAGGGGATCGCTGCGAAGGCTCCCTCCTGGAACCGCGCTGCACCAAACAcc                                                            >  1:137469/1‑76 (MQ=255)
    gTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAcc                                                            >  1:4136676/1‑76 (MQ=255)
      aTCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAt                                                          >  2:4708520/1‑76 (MQ=255)
        cTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                        <  2:1699091/76‑1 (MQ=255)
         tGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATccc                                                       <  2:3376113/76‑1 (MQ=255)
         ggggCGCACAAGCGCTGGAGCAGGGGATCGATGCGATAGTTTCCGCCTGGAACCGCGCAGCACCAAACACCATccc                                                       >  2:4916669/2‑76 (MQ=255)
          gggCGTAGAGGTGCTGGAGCAGGGGCGCGAGGCGAGGGATTCCGCCTGGAAGCGCGCAGCACCACACACCATCCCg                                                      >  1:5340750/1‑76 (MQ=255)
          gggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCg                                                      >  1:1187968/1‑76 (MQ=255)
          gggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCg                                                      >  1:1676601/1‑76 (MQ=255)
            gcgcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                    <  2:1070187/76‑1 (MQ=255)
            gcgcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                    <  1:5577654/76‑1 (MQ=255)
            gcgcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                    <  1:4893701/76‑1 (MQ=255)
               cAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggc                                                 <  1:1549442/76‑1 (MQ=255)
               cAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggc                                                 >  2:1427335/1‑76 (MQ=255)
               cAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAccgc                                                 >  2:668394/1‑76 (MQ=255)
               cAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAACATCCCGAcggc                                                 >  1:4056430/1‑76 (MQ=255)
                aGAAGGGCGGGAGCGGGGGAGCGATACGATGGTTTCCTCCGGGTACCGCGCCGCACTAATCACCATCCCGAcggcg                                                <  2:255335/76‑1 (MQ=255)
                aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggcg                                                <  1:701215/76‑1 (MQ=255)
                aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggcg                                                <  1:758166/76‑1 (MQ=255)
                  aaGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACggcggc                                              <  1:3219869/76‑1 (MQ=255)
                  aaGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACggcggc                                              <  1:905840/76‑1 (MQ=255)
                   aGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACGAAACACCATCCCGACGGCGGCa                                             <  2:4200960/76‑1 (MQ=255)
                   aGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCa                                             <  2:1249025/76‑1 (MQ=255)
                   aGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCa                                             >  1:471356/1‑76 (MQ=255)
                   aGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCa                                             <  2:2413464/76‑1 (MQ=255)
                   aGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCa                                             <  1:4105175/76‑1 (MQ=255)
                   aGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCa                                             <  2:5595107/76‑1 (MQ=255)
                     cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGTAACCTCGCAGCACCAACCGCCCTCCCGACGGCGGCaaa                                           <  1:3344091/76‑1 (MQ=255)
                     cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaa                                           <  2:1279719/76‑1 (MQ=255)
                     cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaa                                           <  1:4326275/76‑1 (MQ=255)
                     cgcTGGAGCAGGGGATCCATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaa                                           >  2:4165759/1‑76 (MQ=255)
                      gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGATCCGCGCAGCACCATACACCATCCCGACGGCGGCaaaa                                          >  1:70334/1‑76 (MQ=255)
                      gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGCCGGCGGCcaaa                                          >  2:434670/1‑76 (MQ=255)
                      gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                          <  1:4557089/76‑1 (MQ=255)
                      gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                          >  2:4977293/1‑76 (MQ=255)
                      gatGGAGCAGGGGATCGATGCGACGGTTTCCTCACGGAAACGAGCAGCACCAAACAACATCACGACGGCGGCgaaa                                          >  1:3672151/3‑76 (MQ=255)
                       gtGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAg                                         <  2:3249373/75‑1 (MQ=255)
                       cTGGGGCAGGGGATCGCTGCGAGGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAg                                         <  2:4579805/76‑1 (MQ=255)
                        tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                        >  1:2933464/1‑76 (MQ=255)
                         ggAGCGGGGGATCGATGCGATGCTTTCCTCCTGGAACGGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGcc                                       <  1:2864763/76‑1 (MQ=255)
                          gAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                      >  2:3503144/1‑76 (MQ=255)
                          gAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                      <  2:1532813/76‑1 (MQ=255)
                           aGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCCAACGCCgg                                     >  2:1678607/1‑76 (MQ=255)
                           aGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                     >  2:5267456/1‑76 (MQ=255)
                           aGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                     <  1:5667518/76‑1 (MQ=255)
                           aGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                     >  1:1191576/1‑76 (MQ=255)
                            gCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCCTCCCGACGGCGGCACAAGcccga                                    >  2:5268375/1‑73 (MQ=255)
                            gCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggg                                    >  1:488333/1‑75 (MQ=255)
                              aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtgg                                  <  1:2960973/76‑1 (MQ=255)
                               ggggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                 <  2:1476762/76‑1 (MQ=255)
                                 ggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGAGGTaa                               >  1:205707/1‑76 (MQ=255)
                                  ggtCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                              >  1:2499857/3‑76 (MQ=255)
                                  gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                              >  2:146696/1‑76 (MQ=255)
                                   aTCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTACCt                             >  2:5339766/1‑76 (MQ=255)
                                    gcgatgcgaggGTTTCCGCCGGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTa                            <  2:5825976/75‑1 (MQ=255)
                                      gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAcc                          >  2:2161819/1‑76 (MQ=255)
                                      gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAcc                          >  2:3447301/1‑76 (MQ=255)
                                         gcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                       >  1:5692688/1‑76 (MQ=255)
                                          cgatgGTTTCCTCCTGGACCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                      <  1:1188328/76‑1 (MQ=255)
                                          cgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                      >  1:5554904/1‑76 (MQ=255)
                                           gatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctct                     >  1:1030149/1‑76 (MQ=255)
                                           gatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctct                     <  2:4657115/76‑1 (MQ=255)
                                            atgGTTTCCTCCTGGACCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCtt                    <  1:5713023/76‑1 (MQ=255)
                                             tgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTc                   >  1:1106509/1‑76 (MQ=255)
                                             tgGTGTCCTCCGGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTc                   <  2:2411207/76‑1 (MQ=255)
                                              ggTTTCCTCCTGGAACGGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTcc                  <  2:4944518/76‑1 (MQ=255)
                                                 ttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctg               <  1:4921149/76‑1 (MQ=255)
                                                 ttCCTCCTGCAACGGCGAAGCACCAAAGACCAGCCCAACGGCGGCAAAAGCCGGTCGCAACGACCTCTCATCCctt               >  2:2685861/1‑75 (MQ=255)
                                                  tcctccTGGACCCGCGCAGCACCAAACCCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgc              <  1:722089/76‑1 (MQ=255)
                                                  tcctccTGGAACCGCGCAGCCCCCCACACCATCCCGACGGCGTCAAACGCCGGGGGTCACTACCTCTCTCCCctgc              >  2:3829987/1‑76 (MQ=255)
                                                  tcctccTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgc              >  2:3796664/1‑76 (MQ=255)
                                                  tcctccTGAACCCGCGCAGCACCAACCACCACCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgc              <  1:2385450/76‑1 (MQ=255)
                                                   cctcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct             >  1:1296077/1‑76 (MQ=255)
                                                   cctcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct             >  1:3013382/1‑76 (MQ=255)
                                                   cctccgGGACCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct             <  1:2631637/76‑1 (MQ=255)
                                                   ccccctGGACCCACGCGCCGCCAAACCCCACCCCGACGGCGGCAAACGCCGGTGGTAACTACCTCTCTTCCctgct             <  1:805057/73‑1 (MQ=255)
                                                    ctcctGGACCCGCGCAGCACCAAACCCCACCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgctg            <  1:232277/76‑1 (MQ=255)
                                                    ctcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGGGGTACCTACCTCTCTTCCccgctg            >  2:115976/1‑76 (MQ=255)
                                                     tcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtgg           >  1:1339387/1‑76 (MQ=255)
                                                      ccgggACCCGCGCAGCACCAAACACCATCCCGCCGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt          <  2:3277353/73‑1 (MQ=255)
                                                       ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                              >  1:5324257/1‑55 (MQ=255)
                                                       ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg         >  1:4224142/1‑76 (MQ=255)
                                                       ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg         >  1:175388/1‑76 (MQ=255)
                                                       ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg         >  1:4460358/1‑76 (MQ=255)
                                                       ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg         >  2:4981491/1‑76 (MQ=255)
                                                        tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg        <  1:5094298/76‑1 (MQ=255)
                                                        tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg        <  1:4977294/76‑1 (MQ=255)
                                                        tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg        >  1:5319174/1‑76 (MQ=255)
                                                        tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg        >  1:2061122/1‑76 (MQ=255)
                                                         ggAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTggg       >  1:3326412/1‑76 (MQ=255)
                                                          gacccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt      <  1:1795667/73‑1 (MQ=255)
                                                          gAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt      >  2:3973329/1‑76 (MQ=255)
                                                           acccGCGCAGCACCAAACACCACCCCGCCGGCGGAAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTa     <  1:2710749/74‑1 (MQ=255)
                                                           aaCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTa     >  1:2449257/1‑76 (MQ=255)
                                                            cccGCGCAGCACCAAACCCCACCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAg    <  2:5424066/75‑1 (MQ=255)
                                                              cgcgCAGCACCAAACACCACCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCg  <  2:1522638/76‑1 (MQ=255)
                                                                    ||                                                                    
GAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCG  >  NC_000913/3950489‑3950626

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