Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,659,541 A→T Y42N (TAT→AAT)  glcK ← glucokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,659,5410AT100.0% 128.5 / NA 39Y42N (TAT→AAT) glcKglucokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (21/18);  total (21/18)

GCTCCATTTACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATATCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTCTGTTCA  >  CP000730/1659405‑1659673
                                                                                                                                        |                                                                                                                                    
gCTCCATTTACTGTACATTTTGCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAATTTCCTTTCAAAAGTGTATTTCCt                                                                                                                                  >  1:921701/1‑141 (MQ=255)
gCTCCATTTACTGTAACTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCt                                                                                                                                  >  1:921704/1‑141 (MQ=255)
   ccATTTACTGTACCATTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTa                                                                                                                               >  1:495763/1‑141 (MQ=255)
      tttACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTa                                                                                                                            >  1:373924/1‑141 (MQ=255)
       ttACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTAt                                                                                                                           <  1:163765/141‑1 (MQ=255)
         aCTGGACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCa                                                                                                                         <  1:215453/141‑1 (MQ=255)
          cTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCACCAAAGGAAACATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAg                                                                                                                        >  1:865347/1‑141 (MQ=255)
                cTTTTTCAAAGTCAACAGGAGCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCTTGTACTATCAGATGTAt                                                                                                                  >  1:651747/1‑141 (MQ=255)
                cTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTAt                                                                                                                  <  1:648509/141‑1 (MQ=255)
                    tttaaaTTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAgt                                                                                                              <  1:165170/138‑1 (MQ=255)
                    ttCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTGCTATCAGATGTATCCgt                                                                                                              >  1:473205/1‑141 (MQ=255)
                        aaGTCAACGGGACCTGGTACACCAATACCTACGCCAAATACATTTGAAAAATTATAATTTTTTTCATTTACTTTTTCAACAAACGAATCATAACTTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGa                                                                                                          <  1:55008/141‑1 (MQ=255)
                        aaGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGa                                                                                                          >  1:579926/1‑141 (MQ=255)
                         aGTCAACAGGACCCGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGaa                                                                                                         <  1:111017/141‑1 (MQ=255)
                          gTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCCGTGTGAAt                                                                                                        >  1:476041/1‑141 (MQ=255)
                          gTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTACAAGTGTGAAt                                                                                                        >  1:168170/1‑141 (MQ=255)
                               cAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCTTCGTTTGCTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGAcc                                                                                                   <  1:440695/141‑1 (MQ=255)
                               cAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGAcc                                                                                                   <  1:142168/141‑1 (MQ=255)
                               cAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGAcc                                                                                                   <  1:470925/141‑1 (MQ=255)
                               cAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGAcc                                                                                                   <  1:417009/141‑1 (MQ=255)
                                       ggTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAAAAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTAACGGTGTGAATAAACCATTTATGt                                                                                           >  1:154801/1‑141 (MQ=255)
                                        gTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTTTACTATCACGTGTATGCGTGTAAATGGACCATTTATCTa                                                                                          >  1:197687/1‑141 (MQ=255)
                                        gTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTCTCAGTGTGAATTGACCAtttattta                                                                                          >  1:197547/1‑141 (MQ=255)
                                          acacCAATTCCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAAt                                                                                        <  1:328659/141‑1 (MQ=255)
                                            acCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTg                                                                                      <  1:655010/141‑1 (MQ=255)
                                                          cAAGTACATTTGAAAAATTATCATTATTTTCTTTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGgtg                                                                        >  1:448291/1‑141 (MQ=255)
                                                          cAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGgtg                                                                        >  1:51946/1‑141 (MQ=255)
                                                            aGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCCGGtgtg                                                                      >  1:433406/1‑141 (MQ=255)
                                                                  tttGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATa                                                                <  1:581954/141‑1 (MQ=255)
                                                                           ttataattatTTTCATTTACTTTTTTAACATACTAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAttttt                                                       >  1:316012/1‑140 (MQ=255)
                                                                           ttataattatTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCAGTTATGTAATGGTTCTATTGCAGGTGTGAAAATTCGTATTTTa                                                       >  1:858679/1‑141 (MQ=255)
                                                                                     tttCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTcc                                             <  1:150759/141‑1 (MQ=255)
                                                                                              cttcttCAACAACCGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGCCCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACAt                                    <  1:249447/141‑1 (MQ=255)
                                                                                                             aaTCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCTGGAGAGAGAATACCTAAACTAAAAGTAGTCCCGCCTACTTCAGCTGCTAAaatat                     >  1:907244/1‑140 (MQ=255)
                                                                                                             aaTCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGCAAATACCTAATTTACAAGACGTCCCGCCTACATCAGCTGCTAAaataa                     >  1:114190/1‑141 (MQ=255)
                                                                                                              aTCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAaataat                    >  1:24571/1‑141 (MQ=255)
                                                                                                                          tttCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTc        <  1:457431/141‑1 (MQ=255)
                                                                                                                           ttCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTCt       <  1:243277/141‑1 (MQ=255)
                                                                                                                                aaGTGTATTTCCCGTACTATCAGATGTATCAGTGTTAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTCTGTTCa  <  1:418508/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
GCTCCATTTACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATATCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTCTGTTCA  >  CP000730/1659405‑1659673

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