Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 950,280 T→G 100% intergenic (+133/‑2) lrp → / → ftsK DNA‑binding transcriptional dual regulator, leucine‑binding/DNA‑binding membrane protein required for chromosome resolution and partitioning

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606950,2800TG95.8% 344.9 / 3.4 94intergenic (+133/‑2)lrp/ftsKDNA‑binding transcriptional dual regulator, leucine‑binding/DNA‑binding membrane protein required for chromosome resolution and partitioning
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base G (36/54);  total (40/54)
Fisher's exact test for biased strand distribution p-value = 3.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
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 102 total aligned reads displayed.

GTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTA  >  REL606/950233‑950329
                                               |                                                 
gtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCtt                                                <  1:8211732/51‑1 (MQ=255)
 tCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCcttgcttg                                               <  1:7494940/51‑1 (MQ=255)
 tCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCcttgcttg                                               <  1:1532421/51‑1 (MQ=255)
 tCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCcttgcttg                                               <  1:8203977/51‑1 (MQ=255)
  ccGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCcttgcttg                                               >  1:8904872/1‑50 (MQ=255)
  ccGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGa                                              >  1:7223454/1‑51 (MQ=255)
   cGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGa                                              <  1:7016021/50‑1 (MQ=255)
   cGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAg                                             <  1:4311424/51‑1 (MQ=255)
   cGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAg                                             <  1:8275063/51‑1 (MQ=255)
    gTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGc                                            >  1:6345755/1‑51 (MQ=255)
     ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGcc                                           >  1:7359329/1‑51 (MQ=255)
     ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGcc                                           <  1:8197012/51‑1 (MQ=255)
     ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGcc                                           >  1:5147063/1‑51 (MQ=255)
      tttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCa                                          >  1:5693527/1‑51 (MQ=255)
      tttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCa                                          >  1:2261468/1‑51 (MQ=255)
      tttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCa                                          <  1:8326852/51‑1 (MQ=255)
      tttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCa                                          >  1:6233808/1‑51 (MQ=255)
       ttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAg                                         >  1:9217598/1‑51 (MQ=255)
       ttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAg                                         >  1:1362611/1‑51 (MQ=255)
       ttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAg                                         <  1:2430001/51‑1 (MQ=255)
        tAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAgg                                        >  1:2900950/1‑51 (MQ=255)
         aGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGa                                       >  1:331241/1‑51 (MQ=255)
          gCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGaa                                      >  1:4886809/1‑51 (MQ=255)
           cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt                                     <  1:2746836/51‑1 (MQ=255)
           cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt                                     >  1:6407312/1‑51 (MQ=255)
           cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt                                     >  1:7117742/1‑51 (MQ=255)
           cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt                                     <  1:2563388/51‑1 (MQ=255)
            atggggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATa                                    >  1:9607473/4‑51 (MQ=39)
            aTCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATa                                    >  1:969801/1‑51 (MQ=255)
            aTCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATa                                    >  1:5708003/1‑51 (MQ=255)
             tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATac                                   <  1:4331214/51‑1 (MQ=255)
             tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATac                                   <  1:6938219/51‑1 (MQ=255)
              cGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt                                     <  1:2743553/48‑1 (MQ=255)
              cGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATaca                                  >  1:8752147/1‑51 (MQ=255)
               gggCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATacac                                 >  1:6984412/1‑51 (MQ=255)
               gggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATacac                                 >  1:5091955/1‑51 (MQ=255)
               gggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATacac                                 <  1:9087648/51‑1 (MQ=255)
               gggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATacac                                 <  1:785812/51‑1 (MQ=255)
               gggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATacac                                 <  1:2530290/51‑1 (MQ=255)
                ggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACt                                <  1:3557081/51‑1 (MQ=255)
                ggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACt                                >  1:3230470/1‑51 (MQ=255)
                ggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACt                                <  1:8263227/51‑1 (MQ=255)
                 gCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTg                               >  1:1270137/1‑51 (MQ=255)
                 gCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTg                               <  1:8222326/51‑1 (MQ=255)
                  cAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGa                              <  1:894068/51‑1 (MQ=255)
                  cAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGa                              <  1:2917695/51‑1 (MQ=255)
                  cAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGa                              >  1:2669661/1‑51 (MQ=255)
                  cAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGa                              <  1:10008854/51‑1 (MQ=255)
                   aGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGaa                             >  1:1768534/1‑51 (MQ=255)
                   aGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGaa                             <  1:5074066/51‑1 (MQ=255)
                   aGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGaa                             >  1:6105426/1‑51 (MQ=255)
                   aGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGaa                             <  1:7395130/51‑1 (MQ=255)
                   aGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGaa                             <  1:6621145/51‑1 (MQ=255)
                    ggAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAg                            <  1:566521/51‑1 (MQ=255)
                    ggAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAg                            <  1:7441407/51‑1 (MQ=255)
                    ggAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAg                            <  1:66377/51‑1 (MQ=255)
                       aaaGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACa                         <  1:2590768/51‑1 (MQ=255)
                       aaaGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACa                         <  1:10265491/51‑1 (MQ=255)
                        aaGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACaa                        <  1:5675704/51‑1 (MQ=255)
                         aGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGGAGACaaa                       <  1:10731549/51‑1 (MQ=39)
                           ccTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaaga                     >  1:9354280/1‑51 (MQ=255)
                            cTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAaagca                    >  1:10262275/1‑49 (MQ=255)
                            cTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaa                    <  1:7822231/51‑1 (MQ=255)
                            cTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaa                    <  1:1073189/51‑1 (MQ=255)
                            cTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaa                    <  1:3546612/51‑1 (MQ=255)
                            cTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaa                    >  1:3788703/1‑51 (MQ=255)
                             tGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATTCACTGAAGACAaagaag                   <  1:6000948/51‑1 (MQ=39)
                             tGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaag                   <  1:6263659/51‑1 (MQ=255)
                             tGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaag                   >  1:2137648/1‑51 (MQ=255)
                              gTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGt                  >  1:6981551/1‑51 (MQ=255)
                              gTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGt                  >  1:7574236/1‑51 (MQ=255)
                               tAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTc                 <  1:9868507/51‑1 (MQ=255)
                                aaCCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTca                <  1:8050986/51‑1 (MQ=255)
                                 aCCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcac               <  1:3702726/51‑1 (MQ=255)
                                 aCCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcac               >  1:2438379/1‑51 (MQ=255)
                                 aCCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcac               <  1:3682939/51‑1 (MQ=255)
                                  ccTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcaca              <  1:2886511/51‑1 (MQ=255)
                                  ccTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcaca              >  1:1473235/1‑51 (MQ=255)
                                  ccTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcaca              <  1:10052629/51‑1 (MQ=255)
                                   cTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGCAGTCACAt             >  1:4601571/1‑51 (MQ=255)
                                   cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt             >  1:4825188/1‑51 (MQ=255)
                                   cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt             <  1:7238482/51‑1 (MQ=255)
                                   cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt             >  1:1575407/1‑51 (MQ=255)
                                   cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt             <  1:8183424/51‑1 (MQ=255)
                                   cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt             >  1:4436778/1‑51 (MQ=255)
                                    tGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAtt            >  1:10709043/1‑51 (MQ=255)
                                    tGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAtt            <  1:4830650/51‑1 (MQ=255)
                                     ggAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTg           <  1:1510737/51‑1 (MQ=255)
                                      gagaGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGa          <  1:6051914/51‑1 (MQ=255)
                                        gagCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACa        <  1:9557948/51‑1 (MQ=255)
                                        gagCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACa        <  1:6208542/51‑1 (MQ=255)
                                         agCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACaa       >  1:8583597/1‑51 (MQ=255)
                                         agCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACaa       <  1:10111154/51‑1 (MQ=255)
                                          gCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACaaa      <  1:10447993/51‑1 (MQ=255)
                                          gCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACaaa      <  1:10390166/51‑1 (MQ=255)
                                           ccTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAg     <  1:4061096/51‑1 (MQ=255)
                                            cttgcttgAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGt    >  1:2772943/1‑51 (MQ=255)
                                            cttgcttgAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGt    >  1:8224298/1‑51 (MQ=255)
                                            cttgcttgAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGt    >  1:9969051/1‑51 (MQ=255)
                                              ttCTTGAGCCAGGAATACACTGAAGACAAAGCAGTCACATTGACAAAGTTa  >  1:3004462/1‑51 (MQ=255)
                                               |                                                 
GTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTA  >  REL606/950233‑950329

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