Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,758,314 A→C 100% K374Q (AAA→CAA)  ygbK → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,758,3140AC100.0% 426.8 / NA 111K374Q (AAA→CAA) ygbKhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (60/51);  total (60/51)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 111 total aligned reads displayed.

TCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATAAAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTA  >  REL606/2758267‑2758361
                                               |                                               
tCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAg                                              <  1:5569782/51‑1 (MQ=255)
tCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAg                                              <  1:6634818/51‑1 (MQ=255)
tCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAg                                              >  1:2160738/1‑51 (MQ=255)
  gCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAggg                                            <  1:1792828/51‑1 (MQ=255)
   cAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGt                                           <  1:4845200/51‑1 (MQ=255)
   cAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGt                                           <  1:956583/51‑1 (MQ=255)
   cAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGt                                           <  1:8017295/51‑1 (MQ=255)
    aGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGtt                                          <  1:9339275/51‑1 (MQ=255)
     ggCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGttt                                         <  1:1680516/51‑1 (MQ=255)
     ggCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGttt                                         <  1:6165217/51‑1 (MQ=255)
     ggCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGttt                                         >  1:3236854/1‑51 (MQ=255)
      gCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTc                                        >  1:5670539/1‑51 (MQ=255)
      gCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTc                                        <  1:5884768/51‑1 (MQ=255)
      gCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTc                                        <  1:10027876/51‑1 (MQ=255)
      gCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTc                                        <  1:5461682/51‑1 (MQ=255)
      gCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTc                                        <  1:7131786/51‑1 (MQ=255)
        ggTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCat                                      >  1:3048287/1‑51 (MQ=255)
        ggTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCat                                      <  1:4966706/51‑1 (MQ=255)
        ggTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCat                                      <  1:5422800/51‑1 (MQ=255)
        ggTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCat                                      >  1:10732221/1‑51 (MQ=255)
         gTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCata                                     <  1:10852547/51‑1 (MQ=255)
          tGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCatat                                    >  1:8334254/1‑51 (MQ=255)
          tGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCatat                                    <  1:8719076/51‑1 (MQ=255)
           gagaCCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATAtt                                   <  1:1392024/51‑1 (MQ=255)
           gagaCCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATAtt                                   <  1:8432624/51‑1 (MQ=255)
           gagaCCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATAtt                                   >  1:3600253/1‑51 (MQ=255)
            agaCCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTg                                  <  1:4391067/51‑1 (MQ=255)
            agaCCTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTg                                  <  1:3060422/51‑1 (MQ=255)
               ccTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGcc                               >  1:1903063/1‑51 (MQ=255)
               ccTCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGcc                               >  1:2169529/1‑51 (MQ=255)
                 tCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCa                             >  1:6584493/1‑51 (MQ=255)
                 tCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCa                             >  1:7579422/1‑51 (MQ=255)
                 tCCGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCa                             <  1:2113857/51‑1 (MQ=255)
                  ccGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCaa                            <  1:9751678/51‑1 (MQ=255)
                  ccGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCaa                            <  1:4267753/51‑1 (MQ=255)
                  ccGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCaa                            <  1:8036391/51‑1 (MQ=255)
                  ccGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCaa                            >  1:2646303/1‑51 (MQ=255)
                   cGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAAc                           >  1:10842992/1‑51 (MQ=255)
                   cGGCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAAc                           >  1:8833507/1‑51 (MQ=255)
                    ggCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAAcc                          >  1:6665678/1‑51 (MQ=255)
                     gCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCa                         >  1:6748541/1‑51 (MQ=255)
                     gCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCa                         >  1:6070689/1‑51 (MQ=255)
                     gCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCa                         <  1:9605718/51‑1 (MQ=255)
                     gCGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCa                         <  1:6346896/51‑1 (MQ=255)
                      cGTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAt                        <  1:2983876/51‑1 (MQ=255)
                       gTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAtt                       >  1:7223140/1‑51 (MQ=255)
                       gTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAtt                       >  1:7031951/1‑51 (MQ=255)
                       gTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAtt                       >  1:9485139/1‑51 (MQ=255)
                       gTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAtt                       >  1:8989389/1‑51 (MQ=255)
                       gTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAtt                       >  1:1079736/1‑51 (MQ=255)
                       gTAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAtt                       >  1:2378739/1‑51 (MQ=255)
                        tAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAttt                      <  1:1421995/51‑1 (MQ=255)
                        tAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAttt                      >  1:5479162/1‑51 (MQ=255)
                        tAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAttt                      >  1:7154578/1‑51 (MQ=255)
                        tAGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCAttt                      >  1:6196982/1‑51 (MQ=255)
                         aGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTc                     >  1:6477867/1‑51 (MQ=255)
                         aGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTc                     >  1:9001333/1‑51 (MQ=255)
                         aGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTc                     >  1:10822208/1‑51 (MQ=255)
                         aGTCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTc                     >  1:5360767/1‑51 (MQ=255)
                           tCACACTGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTccc                   <  1:10094958/51‑1 (MQ=255)
                           tCACACAGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTccc                   <  1:2543534/51‑1 (MQ=255)
                            cacacaGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTcccc                  <  1:5984952/51‑1 (MQ=255)
                             acacaGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTccccc                 <  1:4143063/51‑1 (MQ=255)
                              cacaGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCg                <  1:8407704/51‑1 (MQ=255)
                                caGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGc              <  1:10010014/51‑1 (MQ=255)
                                caGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGc              <  1:7286559/51‑1 (MQ=255)
                                caGAGCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGc              >  1:4234674/1‑51 (MQ=255)
                                 agagCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCg             >  1:2867048/1‑51 (MQ=255)
                                 agagCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCg             >  1:946558/1‑51 (MQ=255)
                                 agagCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCg             >  1:4898376/1‑51 (MQ=255)
                                  gagCCTTGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGt            >  1:1482479/1‑51 (MQ=255)
                                  gagCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGt            <  1:1493138/51‑1 (MQ=255)
                                  gagCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGt            >  1:969292/1‑51 (MQ=255)
                                   agCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTg           >  1:6753971/1‑51 (MQ=255)
                                   agCCTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTg           >  1:171749/1‑51 (MQ=255)
                                     ccTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGcc         >  1:3287115/1‑51 (MQ=255)
                                     ccTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGcc         <  1:233898/51‑1 (MQ=255)
                                      cTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCg        >  1:2595973/1‑51 (MQ=255)
                                      cTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCg        >  1:10285790/1‑51 (MQ=255)
                                      cTGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCg        >  1:7601152/1‑51 (MQ=255)
                                      cTGGGAATACAAGGGCTTCATATTGGCACAACCATTTCCCCCGGCGTg           >  1:2875648/1‑48 (MQ=255)
                                       tGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGc              <  1:9966999/44‑1 (MQ=255)
                                       tGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCg        >  1:5216830/1‑50 (MQ=255)
                                       tGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGt       <  1:8219514/51‑1 (MQ=255)
                                       tGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGt       >  1:8377982/1‑51 (MQ=255)
                                       tGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGt       <  1:1211134/51‑1 (MQ=255)
                                       tGGGAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGt       >  1:7738309/1‑51 (MQ=255)
                                        gggAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTg      >  1:5929550/1‑51 (MQ=255)
                                        gggAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTg      <  1:2923121/51‑1 (MQ=255)
                                        gggAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTg      <  1:10482414/51‑1 (MQ=255)
                                         ggAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTgg     >  1:8719931/1‑51 (MQ=255)
                                          gNATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTggg    >  1:454149/1‑51 (MQ=255)
                                          gAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTggg    >  1:7909830/1‑51 (MQ=255)
                                          gAATACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTggg    >  1:136004/1‑51 (MQ=255)
                                           aaTACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGt   >  1:4226799/1‑51 (MQ=255)
                                           aaTACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGt   >  1:4645046/1‑51 (MQ=255)
                                            aTACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTa  >  1:6826935/1‑51 (MQ=255)
                                            aTACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTa  >  1:6354888/1‑51 (MQ=255)
                                            aTACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTa  <  1:9900969/51‑1 (MQ=255)
                                            aTACAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTa  <  1:3453749/51‑1 (MQ=255)
                                               |                                               
TCGCAGGCGGTGAGACCTCCGGCGTAGTCACACAGAGCCTGGGAATAAAAGGGTTTCATATTGGCCCAACCATTTCCCCCGGCGTGCCGTGGGTA  >  REL606/2758267‑2758361

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