Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 226,159 C→A 100% P163T (CCT→ACT)  gmhB → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606226,1590CA100.0% 248.6 / NA 73P163T (CCT→ACT) gmhBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/41);  total (32/41)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGAC  >  REL606/226111‑226206
                                                |                                               
gcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACt                                               <  1:6415888/51‑1 (MQ=255)
gcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACt                                               <  1:1108340/51‑1 (MQ=255)
  gAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGa                                             <  1:328493/51‑1 (MQ=255)
  gAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGa                                             <  1:8832968/51‑1 (MQ=255)
   aaCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGaa                                            >  1:10194698/1‑51 (MQ=255)
   aaCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGaa                                            <  1:5799142/51‑1 (MQ=255)
   aaCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGaa                                            <  1:9626131/51‑1 (MQ=255)
    aCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAg                                           <  1:9067389/51‑1 (MQ=255)
    aCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAg                                           >  1:7282673/1‑51 (MQ=255)
     cGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGc                                          <  1:3625115/51‑1 (MQ=255)
     cGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGc                                          <  1:4920655/51‑1 (MQ=255)
      gTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCa                                         <  1:7355063/51‑1 (MQ=255)
      gTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCa                                         >  1:10220757/1‑51 (MQ=255)
      gTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCa                                         <  1:4365542/51‑1 (MQ=255)
       tGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAg                                        >  1:3853927/1‑51 (MQ=255)
       tGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAg                                        >  1:8963850/1‑51 (MQ=255)
        gggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGa                                       <  1:4092680/51‑1 (MQ=255)
         ggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaa                                      >  1:4201095/1‑51 (MQ=255)
         ggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaa                                      <  1:7344122/51‑1 (MQ=255)
         ggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaa                                      <  1:561351/51‑1 (MQ=255)
           aaCAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaaaa                                    <  1:7254487/51‑1 (MQ=255)
            aCAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAc                                   <  1:6748187/51‑1 (MQ=255)
            aCAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAc                                   <  1:1939852/51‑1 (MQ=255)
              aaaaGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAcgc                                 >  1:896699/1‑51 (MQ=255)
               aaaGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAcgcg                                <  1:6773983/51‑1 (MQ=255)
               aaaGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAcgcg                                >  1:5235902/1‑51 (MQ=255)
               aaaGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAAcgcg                                >  1:7058187/1‑51 (MQ=255)
                 aGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggc                              >  1:9552423/1‑51 (MQ=255)
                 aGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggc                              <  1:10771587/51‑1 (MQ=255)
                  gTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggcg                             <  1:5027922/51‑1 (MQ=255)
                  gTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggcg                             <  1:3616034/51‑1 (MQ=255)
                   tGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGcggcgg                            <  1:2209481/51‑1 (MQ=255)
                   tGCTGGTGCGTACGGGCAAACCTATTACGACTGAAGCAGAAAACGcggcgg                            >  1:4016315/1‑51 (MQ=255)
                    gCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGa                           <  1:8583005/51‑1 (MQ=255)
                     cTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGAt                          >  1:10929900/1‑51 (MQ=255)
                      tGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGAtt                         <  1:2381747/51‑1 (MQ=255)
                      tGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGAtt                         >  1:2175915/1‑51 (MQ=255)
                       ggTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTg                        <  1:3690026/51‑1 (MQ=255)
                        gTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTgg                       >  1:3043958/1‑51 (MQ=255)
                         tGCGTACGGGAAAACCTATTACGACTGAAGCAGAAAACGcggcgg                            >  1:796270/1‑45 (MQ=38)
                          gCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTggg                      >  1:8037512/1‑50 (MQ=255)
                            gTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTggg                      >  1:1469603/1‑48 (MQ=255)
                            gTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGgtgt                   <  1:9677372/51‑1 (MQ=255)
                             tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                  >  1:2680817/1‑51 (MQ=255)
                             tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                  <  1:1893633/51‑1 (MQ=255)
                             tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                  >  1:7171823/1‑51 (MQ=255)
                             tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                  >  1:3006609/1‑51 (MQ=255)
                             tACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGtt                  <  1:3621892/51‑1 (MQ=255)
                              aCGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTggg                      >  1:5022179/1‑46 (MQ=255)
                              aCGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTa                 >  1:6507423/1‑51 (MQ=255)
                              aCGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTa                 >  1:2837205/1‑51 (MQ=255)
                                gggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTaaa               <  1:8709497/51‑1 (MQ=255)
                                gggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTaaa               <  1:7033026/51‑1 (MQ=255)
                                gggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTaaa               >  1:5723656/1‑51 (MQ=255)
                                  gTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATa             <  1:4993723/51‑1 (MQ=255)
                                   tAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAg            <  1:119403/51‑1 (MQ=255)
                                     aaCCTATTACGACTGAAGCAGAAAACgcgg                               >  1:1301711/1‑30 (MQ=38)
                                     aaCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAAt              >  1:6829306/1‑47 (MQ=255)
                                      aCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCt         <  1:5624080/51‑1 (MQ=255)
                                      aCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCt         >  1:10305882/1‑51 (MQ=255)
                                      aCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCt         >  1:8304853/1‑51 (MQ=255)
                                      aCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCt         <  1:671791/51‑1 (MQ=255)
                                       ccTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTg        <  1:5883/51‑1 (MQ=255)
                                        cTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTgg       <  1:10665549/51‑1 (MQ=255)
                                          aTTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCa     <  1:6060314/51‑1 (MQ=255)
                                          aTTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCa     <  1:5191097/51‑1 (MQ=255)
                                           ttACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAg    <  1:6316799/51‑1 (MQ=255)
                                           ttACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAg    >  1:2906367/1‑51 (MQ=255)
                                           ttACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAg    >  1:3349489/1‑51 (MQ=255)
                                           ttACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAg    <  1:4316866/51‑1 (MQ=255)
                                            tACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGa   >  1:1951912/1‑51 (MQ=255)
                                            tACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGa   <  1:6690845/51‑1 (MQ=255)
                                             aCGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGAc  >  1:7936821/1‑51 (MQ=255)
                                                |                                               
GCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGAC  >  REL606/226111‑226206

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